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Liquid Crystal Large Screen Projection Displays 


@ Akihiro Mochizuki @ Yasuhiro Yoneda @ Toshiaki Narusawa 


A liquid crystal large screen projection display for presentations 
is reviewed. Two types of addressing mode, an optically addressed 
and an electrically addressed projection display are introduced. 
Several types of electrically addressed liquid crystals are possible 
for the presentation device. Those currently available include the 
supertwisted nematic (STN) and nematic-cholesteric phase tran- 
sition (NCPT) liquid crystal projectors. The applications of these 
for a still image display device are compared. 
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Recent Advances in the Application of Transmission 
Electron Microscopy to the Evaluation of IIIl-V Compound 
Semiconductor Thin Films 


@ Osamu Ueda 


This paper reviews recent advances in the application of trans- 
mission electron microscopy to the evaluation of structures, 
interfaces and defects in compound semiconductor thin films. It 
has been shown microscopically for the first time that the growth 
of (GaAs)/(GaP), strained layer superlattices by atomic layer 
epitaxy proceeds in a precisely layer-by-layer manner. Detailed 
characterizations of interfaces and defects in AlAs, GaAs and GaP 
on Si, also GaAsSb and Fe-doped InP on InP are described; later, 
the mechanisms which generate such defects and techniques for 
their elimination are discussed. Studies of atomic ordering in 
InGaP and InGaAs are presented, with reference to the mecha- 
nisms which generate ordered structures and their influence on the 
Properties of the material. 
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Pulsed-Jet Epitaxy of IIl-V Compounds 


@ Masashi Ozeki @ Nobuyuki Ohtsuka @ Yoshiki Sakuma 


Pulsed-jet epitaxy (PJE) for IIl-V compounds (GaAs, InP, GaP, 
and InAs) and its possibilities for process technology are described. 
A model involving selective reactions has been proposed to explain 
the self-limiting process in PJE. A reduction of carbon contamina- 
tion has been achieved for GaAs growth using PJE and carrier 
concentrations ranging from 10'* cm™ to 107° cm™® for n-type 
GaAs, and from 10!5 cm”? to 107! cm™® for p-type GaAs. These 
can be obtained by controlling growth conditions and doping 
levels. The growth of (GaAs), (GaP), superlattices has demon- 
strated the high potential of PJE in superlattice growth. PJE offers 
unique possibilities for low temperature growth, selective growth, 
and uniform thickness growth. 
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Metalorganic Vapor Phase Epitaxy Using Organic Group-V 
Sources 


@ Toshihide Kikkawa @ Hitoshi Tanaka @ Junji Komeno 


Safer MOVPE growth of GaAs, AlGaAs, and InGaAs on 3-inch 
GaAs substrates was investigated using organic group-V sources. 

Tertiarybutylarsine (TBAs) and monoethylarsine (EtAs) were 
used with a horizontal atmospheric pressure reactor. GaAs and 
AlGaAs grown using TBAs were of high purity and extremely 
uniform, equivalent to those using arsine. With either silane or 
disilane doping, Si incorporation efficiency using TBAs or EtAs 
was always higher than that using arsine. The group-V pyrolysis 
simulation implies that H,AsSiH, dominates in the gas phase 
when using TBAs, corresponding to high doping efficiency. High 
quality, high electron-mobility transistors were achieved with 
TBAs. These results verify that epitaxial layers grown using 
TBAs are of sufficiently high quality for device applications. 
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Molecular Dynamics Simulation on an AP1000 Dis- 
tributed Memory Parallel Computer 


@ Hiroyuki Sato @ Yasumasa Tanaka @ Toru Yao 


Parallel processing techniques have been applied to the molecular 
dynamics calculation module of AMBER developed by the Uni- 
versity of California. The target machine is the AP1000 distributed 
memory parallel computer developed at Fujitsu Laboratories 
Limited. AP1000 consists of up to 1024 processor elements 
connected with three different networks. To obtain a higher 
degree of parallelism and better load balance between processors, a 
Particle division method was developed to randomly allocate 
Particles to processors. Experiments showed that a problem with 
41095 atoms can be processed 226 times faster with a 512-pro- 
cessor AP1000 than by one processor. 
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Automated Visual Inspection System for Bonded IC Wires 


@ Hiroyuki Tsukahara @ Masato Nakashima 
@ Takehisa Sugawara 


This paper discusses an automated visual inspection system for 
bonded IC wires. The system uses high-contrast image capture, an 
algorithm for accurately measuring the bonding balls, and a wire 
inspection algorithm. The high-contrast image capture system 
consists of a wide-area camera and two types of lighting optics, 
one for ball sensing and the other for wire sensing. The algorithm 
for measuring the bonding balls is based on morphological tech- 
niques, and the wire inspection algorithm is based on a border 
following algorithm. The automated inspection system measures 
ball diameters to within +5 um accuracy, which corresponds to 
half a pixel of the captured picture. The system takes 0.2 seconds 
to inspect a wire and ball. Combining the inspection system with 
an automatic wire bonder enables a fully automatic bonding 
system. 
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Low Drive Voltage and Low Chirp Modulator Integrated 
DFB Laser Light Source for Multiple-gigabit Systems 


@ Toru Watanabe @ Keiji Sato @ Haruhisa Soda 


This paper discusses a low drive voltage and low chirp modulator 
integrated DFB laser (MI-DFB laser). In this research, the composi- 
tion of the absorption layer was adjusted experimentally to 
obtain a low drive voltage, and the absorption layer thickness 
was analytically optimized. The drive voltage for an extinction 
ratio of -13dB is found to be -3 V, which is small enough for 
high-speed IC drivers. A SiN film on the front facet is optimized 
for an antireflective (AR) coat, and gives a low chirp width of 
0.016 nm at -15dB (0.004nm at -3 dB) under 10 Gbit/s NRZ 
pseudo random modulation. The low drive voltage and low chirp 
performance is promising for use in the next generation of multi- 
ple-gigabit long-haul optical fiber transmission systems. 


UDC 548-144:681.327 


Liquid Crystal Large Screen Projection 
Displays 


@ Akihiro Mochizuki @ Yasuhiro Yoneda @ Toshiaki Narusawa 
(Manuscript received September 9, 1991) 


A liquid crystal large screen projection display for presentations is reviewed. 

Two types of addressing mode, an optically addressed and an electrically addressed projec- 
tion display are introduced. Several types of electrically addressed liquid crystals are possible 
for the presentation device. Those currently available include the supertwisted nematic 


(STN) and nematic-cholesteric phase transition (NCPT) liquid crystal projectors. 
The applications of these for a still image display device are compared. 


1. Introduction 

Large projection displays are currently popu- 
lar for presentations, public announcements, and 
showing moving pictures. For presentations, a 
liquid crystal display (LCD) pad is now com- 
monly used together with a personal computer 
and an overhead projector. The LCD pad can 
conveniently be used with a conventional over- 
head projector and a personal computer, this 
kind of system may be called an electronic over- 
head projector (EOHP). Since the LCD pad is 
light and thin, it can easily be carried and by 
using the simple optical system of an overhead 
projector, a low cost LCD projector is possible. 
However, such a display requires high informa- 
tion content, high resolution, a large display 
area, a high luminance screen and good compat- 
ibility with various office automation equip- 
ment. To meet these demands, we proposed a 
projection display using nematic-cholesteric 
phase transition (NCPT) type liquid crystals!). 
Consequently, several types of liquid crystal 
projection display for EOHPs have been pro- 
posed and realized*)*). The most popular 
liquid crystal projection display is the so-called 
electronic overhead projector using a _ super- 
twisted nematic (STN) LCD, and this projection 
display is quite portable. A high resolution pro- 
jection display using laser beam addressed 
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smectic liquid crystals is now also available. 

Despite such progress, the liquid crystal 
projection display is still in the development 
stage due to such drawbacks as poor screen 
luminance and information content. We have 
proposed and achieved an NCPT projection 
display which overcomes these drawbacks. We 
now compare and discuss the characteristics of 
the several types of liquid crystal large screen 
projection displays, particularly those used for 
meetings, education, conferences and public 
announcements. 


2. LCD modes for projection displays 

Several types of liquid crystal drive mode are 
first introduced, then the advantages and disad- 
vantages of the LCDs for EOHP applications are 
discussed. 

Liquid crystal projection displays can be 
classified into two types, an optically addressed 
LCD projector, and an electrically addressed 
LCD projector. Generally, the laser beam ad- 
dressed type offers higher information content 
and a higher luminance screen image than the 
electrically addressed type. However, the elec- 
trically addressed type is more compact and 
less expensive. 
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Projection lens 


Fig. 1—Structure of laser beam addressed liquid crystal 
projection display system. 


Orientation films 


Light absorption layer ITO 
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Read out 


Glass substrate Glass substrate 


A 
ITO 
Dielectric mirror 


Liquid crystal layer 
Fig. 2—Structure of liquid crystal panel for laser beam 
addressed display. 


2.1 Optically addressed projector 

2.1.1 Laser beam addressed projector 

In 1972, Kahn, F.J. proposed a laser beam 
addressed projection display using cholesteric 
liquid crystals®””’. Since then, improvements 
have been made for this type. Currently, smectic 
A and C phase liquid crystals are used for this 
type of projection display. Figure | shows the 
structure of the projection display system and 
Fig. 2 shows the structure of the liquid crystal 
display. The laser beam addressed liquid crystal 
projector uses phase transition induced by 
thermal transfer from the laser beam to the 
liquid crystals. The projector is initially trans- 
parent since the liquid crystal alignment in the 
panel is homeotropic as shown in Fig. 3. 

The laser beam is scanned to a particular 


Electric field 


Fig. 3—Liquid crystal molecular alignment of the laser 
beam addressed projection system: 
Transparent state. 


Fig. 4—Liquid crystal molecular alignment of the laser 
beam addressed projection system: 
Scattering state. 


point in the panel, and the liquid crystal at that 
point is heated to the isotropic phase. When an 
electric field is applied to the panel and the laser 
beam removed, the heated liquid crystal cools 
with a parallel molecular alignment, giving a 
transparent state. When an electric field is not 
applied, the liquid crystal alignment becomes 
irregular with no alignment force during the 
cooling process, resulting in the light scattering 
as shown in Fig. 4. 

The transparent and light scattering states 
are projected on to a large screen by the reflec- 
tive optical system shown in Fig. 1. In this sys- 
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tem, the transparent state in the liquid crystal 
panel provides a bright image on the screen and 
the scattering state provides a dark image. Color 
images can be provided by using color filters. 

The advantages of this display are a bright 
screen, high resolution and high information 
content. The laser beam addressed system uses 
light scattering induced by thermal transfer, thus 
there is no need to use polarizers which reduce 
the screen luminance. Also, since the laser beam 
addressed system uses the liquid crystal’s memo- 
ry effect, which preserves both light transparent 
and scattering states, there is no limitation to 
the information content. A system with nine 
million pixels was proposed by Usubuchi in 
1987®). However, the speed of writing is slow, 
and the system is heavy and bulky, and thus 
expensive. Laser beam addressing works one 
point at a time, thus the writing time is propor- 
tional to the number of pixels. To write nine 
million pixels, for instance, it takes 10 us times 
nine million, or 90s for a whole screen. Fur- 
thermore, heat transfer from the laser beam to 
the liquid crystal requires a high power laser; 
rapid heat transfer requires a several hundred 
mW class laser. The laser beam scanning requires 
beam transducers and beam modulators, and 
these make the liquid crystal projection system 
bulky and expensive. 

The laser beam addressed liquid crystal pro- 
jector has an extremely high information 
content and resolution with a bright screen, and 
it produces shades of gray by dithering. Such a 
high resolution image with gray shades is consid- 
ered suitable for computer aided design (CAD) 
displays. To improve the writing speed, a laser 
beam addressed system using a chiralsmetic C 
phase liquid crystal has been proposed by 
Bone”). The fast response of the chiralsmetic C 
phase liquid crystal and the limited information 
content combine to achieve a video rate liquid 
crystal projection display. Moreover, the con- 
trast ratio of this display is as high as 100:1 due 
to the higher order liquid crystal molecular 
alignment of the chiralsmetjc C phase. 

2.1.2 CRT addressed projector 

A CRT addressed liquid crystal projector has 
been proposed for moving picture images by 


FUJITSU Sci. Tech. J., 28, 1, (March 1992) 


Photoconduction layer —_ Orientation films 
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Dielectric mirror Liquid crystal layer 


Fig. 5—Structure of CRT addressed liquid crystal 
projector. 


Koda!®. The structure of the liquid crystal 
panel is the same as the laser beam addressed 
projector, but the writing method is different. 
This system uses the rapid scanning of a CRT 
beam to give a moving picture. While CRT ad- 
dressing can produce video rate projection, the 
resolution is not so good because of divergence 
of the CRT beam. To suppress beam divergence, 
an optical fiber guide was used between the CRT 
and LCD, as shown in Fig. 5. Using the optical 
fiber guide, Koda realized a 640 x 480 pixel 
video rate liquid crystal projector. The liquid 
crystal panel used tunable birefringence with a 
dielectrically positive anisotropy nematic liquid 
crystal material. This mode requires two polar- 
izers. Thus the screen luminance is as low as a 
conventional twisted nematic liquid crystal 
display and the CRT writing equipment makes it 
a rather heavy and bulky system. 

This display is being used for a spatial light 
modulator (SLM) with optical processability in 
conjunction with an optical fiber guide. SLMs 
have great potential optically addressed liquid 
crystal devices for optical switching systems. 


2.2 Electrically addressed LC projector 

There are several types of electrically ad- 
dressed liquid crystal projector. Generally, the 
electrically addressed projector is smaller and 
lighter than the optically addressed projector, 
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Table 1. Electrically addressed liquid crystal drive modes for a projection display 


STN DAP SH NCPT Ferroelectric AM-TN 
Screen luminance Poor Poor Poor Excellent Fair Poor 
Aperture ratio Good Good Good Good Good Fair 
Information 640 1 200 2640 22240 =2000 1 200 
content (dot) x 480 x 1000 x 480 X 2240 x 2000 x 1000 
Drive voltage (V) 20 60 60 25 25 10 
Production cost Good Fair Fair Good Fair Poor 
eae) Good Good Good Good Good Excellent 
Contrast ratio Good Good Good Good Fair Excellent 
Pacer carne Poor Poor Fair Poor Fair Good 


and is thus suitable for EOHPs. 

In principle every liquid crystal drive mode 
can be used in the liquid crystal projector. 
Several electrically addressed liquid crystal drive 
modes which are suitable for projection display 
are introduced below. The advantages and disad- 
vantages of each liquid crystal drive mode are 
summarized in Table 1. 

2.2.1 Tunable birefringence mode 

In this mode, the liquid crystal molecular 
alignment is controlled by an electric field. The 
liquid crystal’s birefringence changes as its align- 
ment changes, and this property is used for 
polarized light switching. A dielectrically posi- 
tive nematic liquid crystal material is used with 
homogeneous alignment, a dielectrically negative 
liquid crystal material is used with homeotropic 
alignment. A hybrid molecular alignment panel, 
whose upper substrate has homeotropic align- 
ment and lower substrate has homogeneous 
alignment, also demonstrates the tunable 
birefringence effect’ a 
1) Nematic LC with homogeneous alignment 

In this liquid crystal molecular alignment, a 
dielectrically positive liquid crystal material is 
used. The most popular device of this type is 
supertwisted birefringence effect (SBE) or 
supertwisted nematic (STN) LCDs. The liquid 
crystal’s birefringence is controlled by an exter- 
nal applied voltage. The transmission of incident 
light is controlled by the retardation And of the 
liquid crystal which is sandwiched between two 


polarizers. The STN type of liquid crystal is now 
used for projection displays, and its relatively 
high information content makes it attractive for 
presentation displays. Some of the STN projec- 
tion systems are well designed and equipped 
with sophisticated devices, a detailed discussion 
is given in the next chapter. 
2) Nematic LC with homeotropic alignment 

This type of LCD makes use of the deforma- 
tion of the vertical aligned phase (DAP)!”?. It 
uses a dielectrically negative liquid crystal mate- 
rial with homeotropic alignment. When there 
is no applied voltage, linearly polarized incident 
light passes through the liquid crystal medium 
with no polarization change, hence it cannot 
pass through the polarizer on the viewer’s side, 
and this makes the dark state. When a drive volt- 
age higher than the threshold is applied, the 
molecular optical axis inclines in proportion to 
the drive voltage. The polarized axis of the 
incident light and the molecular optical axis are 
then no longer coincident, resulting in 
birefringence of the liquid crystal. This 
birefringence rotates the polarization axis of the 
incident light, and the rotational axis allows 
transmission through the polarizers on the 
viewer's side. The intensity of the transmitted 
light is dependent on the external applied volt- 
age. 

Tunable birefringence is considered suitable 
for projection displays. However, the bire- 
fringence of liquid crystal material is usually 
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Linear polarizer 
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Liquid crystal panel ([|\ 


Optical compensator 


Linear polarizer i 


Fig. 6—Display configuration of super homeotropic (SH) 
mode with optical compensator. 


quite sensitive to the ambient temperature, so 
this mode may need strict temperature control 
when used for OHP display. However, heat 
absorption by the two polarizers makes temper- 
ature control difficult. 

2.2.2 Electrically controlled birefringence 

mode 

Electrically controlled birefringence (ECB) 
mode with homeotropic alignment has good 
multiplexing capability, achromaticity and gray 
scale performance!) particularly when used for 
projection. However, the viewing angle of an 
ECB mode LCD is narrow, thus it is quite dif- 
ficult to use it for a flat panel display which is 
viewed directly. This type of LCD may therefore 
be suitable for projection display'*). Recently, 
an ECB mode with an optical compensator has 
been demonstrated for a direct viewing dis- 
play!*15), The display configuration is shown 
in Fig. 6. The liquid crystal panel and the optical 
compensator are sandwiched between crossed 
polarizers, and the LCD uses a dielectrically 
negative liquid crystal with homeotropic align- 
ment. This type of LCD is called super homeo- 
tropic (SH). Generally when the ECB with 
homeotropic alignment is perpendicular to the 
substrate, and crossed polarizers are at 45° to 
the molecular tilt direction, no birefringence 
occurs, resulting in the dark state. When an 
external voltage is applied to the LCD, the liquid 
crystal molecules tilt and light is transmitted. 

However, the molecular alignment of the 
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liquid crystal sometimes deviates from perpen- 
dicular to the substrate when an external voltage 
is applied, due to reverse tilt of the liquid crys- 
tal. To avoid this molecular alignment oblique to 
homeotropic alignment is needed, but this 
induces both light leakage and a narrow viewing 
angle, which reduces the legibility. However, a 
narrow viewing angle is not a disadvantage for 
projection displays. The fixed optical system of 
an OHP does not need a wide viewing angle, but 
light leakage due to deviations reduces the 
contrast ratio of the display. Light leakage is an 
intrinsic problem of the ECB types of projec- 
tion display. 

The response of the super homeotropic LCD 
is rapid compared to tunable birefringence 
LCDs. The SH LCD also gives a high information 
content because of the steep threshold in the 
electrooptical relation. However, this LCD has a 
high drive voltage, reaching 60V for a 
640 x 400 pixel panel. The high drive voltage 
raises the cost, particularly for a high informa- 
tion content LCD. However, the ECB mode 
LCD has a great potential for high contrast 
projection displays. Further improvements in 
the materials to overcome both the deviations 
and the high drive voltage will be quite impor- 
tant to achieve an ECB mode projection display 
in practice. 

2.2.3 Nematic bistable LCDs 

Heffner!® suggested that when a nematic or 
cholesteric liquid crystal material is placed 
between two surfaces, the directors can be 
constrained throughout the cell to certain equi- 
librium configurations. The configuration of the 
liquid crystal directors throughout the cell is 
determined by the surface boundary conditions, 
the material constants such as helical pitch, and 
the cell gap. For certain geometries and material 
parameters, two or more configurations are 
taken as the solutions to the Oseen-Frank equa- 
tions!”). In these configurations, a first-order 
phase change may occur from the state of higher 
free energy to the configuration with the lower 
free energy. This first-order phase change pro- 
vides the bistable cholesteric twist cell'®?. 

In this cell, bistability can occur between 
two topologically equivalent states with differ- 
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Fig. 7—Nematic bistable configuration. Orientation of 
molecules in three states of holding voltage cell. 
The up and down orientations are the two 
bistable states. 
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Fig. 8—Electrooptical property of the NCPT liquid 
crystal display. 


ent twist behavior. A dielectrically positive 
nematic liquid cystal material with added 
cholesteric liquid crystal is used with pretilted 
homogeneous alignment. The two states, which 
are described as Up or Down, may exist either 
with No field or in the presence of a holding 
field as shown in Fig. 7. Barberi!®? also reported 
the nematic bistable configuration. He suggests 
that the bistability is caused by surface anchor- 
ing and flexoelectricity. Sandwiched between 
two polarizers, these configurations can be used 
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a) Focal conic b) Homeotropic c) Quasi-homeo- 
state tropic state 
Fig. 9—Liquid crystal molecular alignment of the NCPT 

display. 


for LCDs. However, reliable bistable behavior is 
required for a practical LCD, and particularly 
for projection displays, bistability will govern 
the performance of both the LCD and the liquid 
crystal projection display. 

2.2.4 Nematic-cholesteric phase transition 

(NCPT) mode 

A nematic-cholesteric phase transition type 
liquid crystal mixture with positive dielectric 
anisotropy changes state when the applied volt- 
age is increased as shown in Fig. 8. It moves 
from a cloudy cholesteric state, cholesteric 
phase F, to a transparent state, nematic phase H. 
When the voltage is decreased, this mixture 
moves from phase H to a metastable transparent 
state, nematic phase H’, then to phase F. In 
phase F, the liquid crystal in the panel has the 
helical structure shown in Fig. 9a). Because the 
axis of the helix is parallel to the glass sub- 
strates, incident light is scattered, making phase 
F cloudy. In phase H, the liquid crystal does 
not have a helical structure but has a home- 
otropic structure instead, as shown in Fig. 9b). 
Incident light passes through the homeotropic 
structure, making phase H transparent. In phase 
H'’, the structure of the liquid crystal layer is 
homeotropic. However, in the center, the liquid 
crystal molecules are thought to be slightly 
inclined as shown in Fig. 9c). Phase H’ is a 
metastable state, a supersaturated state. Super- 
saturation might stabilize hysteresis in the 
cholesteric-nematic phase transition'”?®. Thus 
the stability of the hysteresis effect is though to 
be affected by the specific balance between 
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Fig. 10—Bistability of NCPT mode. 


Fig. 11—Bistable projection configurations of the 
ferroelectric liquid crystal display. 


surface anchoring and the thermal properties of 
the helical structure. 

The width of the hysteresis effect is A, as 
shown in Fig. 10. This is the difference between 
the voltage V,?° which gives 20 percent trans- 
mittance in the cholesteric to nematic phase 
transition, and the voltage V4°° which gives 
90 percent in the nematic to cholesteric phase 
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Fig. 12—High resolution mechanism of ferroelectric 
liquid crystal display. 


Fig. 13—Nematic liquid crystal molecular alignment 
along with lateral electric field. 


transition. The drive voltage Vg is set between 
V,7° and Vq°° as shown in Fig. 10. 

Addressing of the NCPT panel is done using 
electrooptical hysteresis. After initializing all the 
pixels, a writing task is done using the electro- 
optical behavior. This is done one line at a time, 
thus the writing time of the whole screen is 
proportional to the number of scanning lines. 
Since the average writing time for each scanning 
line is 2.2 ms, it takes 2.2 s (2.2 ms x | 000 lines) 
to write 1 000 scanning lines. 

2.2.5 Ferroelectric LC mode 

Unlike nematic liquid crystal materials, 
ferroelectric liquid crystal shows spontaneous 
polarization coupled with externally applied 
voltage. Strong coupling with the applied voltage 
gives a sub-millisecond response time for the 
ferroelectric liquid crystal display. Furthermore, 
the ferroelectric liquid crystal molecule has 
bistable polarization configurations as shown in 
Fig. 11, and this bistability makes a memory 
type drive possible. Thus the ferroelectric liquid 
crystal display has no limitation of information 
contents with a direct drive panel. The liquid 
crystal molecular movement is also parallel to 
the substrates, as shown Fig. 12, and parallel 
movement and bistability combine to give a high 
resolution of the display. A fine, transparent 
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electrode pattern provides the lateral electric 
field. 

Since nematic liquid crystal materials which 
exhibit induced polarization are polarized longi- 
tudinally, the liquid crystal molecules at the 
outer edge of the electrode align with the lateral 
electric field as shown in Fig. 13. The optical 
axis of a nematic liquid crystal is the long 
molecular axis, hence the light transmission 
changes continuously with the lateral electric 
field. 

In contrast, the spontaneous polarization of 
a ferroelectric liquid crystal is perpendicular to 
the molecular long axis, thus even though there 
is a lateral electric field, the molecular alignment 
stays almost perpendicular to the substrate. 
Therefore, the ferroelectric liquid crystal display 
has a high resolution and high contrast ratio?!. 

The ferroelectric liquid crystal display needs 
two polarizers to reduce screen luminance for 
projection display. However, high resolution is 
necessary for projection displays. Since the fer- 
roelectric liquid crystal display has both high 
information content and high resolution, a good 
large projection screen image can be produced. 
Iwai?” reported a four million pixel full color 
projection display using a ferroelectric liquid 
crystal. This projection display consisting of 
three ferroelectric liquid crystal panels, and each 
panel switched light. Red, green and blue lights 
were switched by each liquid crystal panel, and 
full color was achieved by dithering. 

2.2.6 Other liquid crystal modes 

Several other types of liquid crystal drives 
may be used for projection displays. A polymer 
dispersed liquid crystal (PDLC)?*) may be an 
excellent candidate due to its light scattering 
ability. It can be made into sheet film for a large 
liquid crystal display viewed directly, such as for 
may also be used with public announcement 
boards’. A large screen projection display may 
also be used with a polysilicone thin-film transis- 
tor (TFT) panel*®), 

The z-cell?® is another candidate for a large 
liquid crystal projection display, it offers fast 
optical switching at 1 ms. The m-cell is being 
used for a time sequential optical shutter in 
conjunction with a CRT display?”?. The z-cell 


may require TFTs for a high information con- 
tent display, and the need to use polarizers 
makes it hard to use it for large projection dis- 
plays. 


3. Comparison of currently available drive 
modes 

There are three types of liquid crystal drive 
mode currently in use: a nematic-cholesteric 
phase transition (NCPT) type liquid crystal, a 
supertwisted nematic (STN) liquid crystal, and 
an active matrix (AM) drive twisted nematic 
(TN) liquid crystal. Both the NCPT and the STN 
type are driven by a multiplexed system in a 
simple matrix panel. The active matrix TN is 
driven by thin-film transistors (TFTs). 

The characteristic of a simple matrix and an 
active matrix for a liquid crystal projector are 
first compared, then application in the electron- 
ic overhead projector is discussed. 

The active matrix driven TN liquid crystal 
projector can produce images at video speed, so 
it is suitable for a video projector. However, the 
aperture ratio of the active matrix panel is lower 
than that of the simple matrix panel, and for a 
high resolution display such as the one million 
pixel projector, the aperture ratio is lower than 
50 percent. The low aperture ratio and two 
polarizers significantly reduce screen luminance 
which makes the legibility of the image poor. 

In contrast, a simple matrix panel has a large 
aperture ratio. The current LCD production 
technology gives a much higher productivity of 
an A4-size LCD panel with 640 x 400 pixels 
than for an active matrix panel. An A4-size LCD 
panel is required for an electronic overhead 
projector for compatibility with conventional 
OHP transparencies. 

Since the electronic overhead projector is 
used mainly for meetings, conferences, public 
announcements and so on, a character image is 
mainly required. For the first stage, a still black 
and white image is necessary. Then, a still multi- 
color, and a still full color image should be 
provided. Next, a video rate multicolor and full 
color image will be required, and finally, a real- 
time three dimensional full color projector can 
be expected. 
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In the next section, the characteristics of the 
two types of currently available EOHPs using a 
simple matrix liquid crystal panel for presenta- 
tions are compared. 


3.1 NCPT projection display 

The NCPT liquid crystal projector was 
reported by Mochizuki in 1985!). This projec- 
tion display needs no polarizer, so screen lumi- 
nance is high. Without polarizers, not only is the 
screen brighter, but also is absorbed the liquid 
crystal panel heats up less since less projected 
light is absorbed. The optical system of the 
conventional overhead projector gives poor 
lighting uniformity. Sometimes, the light inten- 
sity in the center of the image is twice as bright 
as the outer part. Therefore, a liquid crystal 
projector equipped with polarizers should not be 
used for many hours continuously. Continuous 
usage may cause part of the image to deteriorate 
due to high temperatures in the liquid crystal 
panel. 

The total transmittance of the NCPT panel is 
more than 80 percent, which depends on the 
aperture ratio and also thus affects the resolu- 
tion of the panel. The productivity of the NCPT 
panel is also much higher than the STN panel. 
The NCPT needs no rubbing treatment for liquid 
crystal molecular alignment, whereas, the STN 
panel requires careful rubbing for correct align- 
ment. The molecular orientation of the liquid 
crystal is one of the most sensitive thus most 
difficult processes in the production of LCDs. 
Therefore, the NCPT LCD panel production 
process, which does not require rubbing, is more 
efficient. Particularly for the production of large 
panels, such as A4-size, non-rubbing is preferable 
to produce a uniform molecular orientation. 

The other advantage of the NCPT projector 
is that it has no limitation in information con- 
tent. Yabe?®? reported a five million pixel NCPT 
projection display in 1991 (see Fig. 14). The 
bistability of the NCPT type gives an extremely 
high information content, as high as the laser 
beam addressed liquid crystal projector system is 
quite simple, as with conventional STN LCDs, 
and is therefore inexpensive. Moreover, since the 
NCPT needs no polarizers, the screen can be as 
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Fig. 14—The 5-million pixel NCPT liquid crystal 
projection display. 
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Fig. 15—Structure of double-layer NCPT liquid crystal 
projection display. 


a laser beam addressed system. To produce five 
million pixels high resolution projection uniform 
molecular alignment over a large area is done 
without rubbing and no heat up without polar- 
izers. Although this display uses the bistability 
of the NCPT liquid crystal, the extremely high 
resolution enables a fine gray scale image by 
dithering as shown in Fig. 14. 

The NCPT projector can also produce color 
images without using color filters, since it 
controls the direction of the incident light. 
The direction of scattered light is governed by 
the helical pitch of the liquid crystal and the 
retardation And of the material, where An is 
the birefringence of the liquid crystal and d 
is the panel gap. Due to birefringence the 
direction of the incident light is dependent on 
the wavelength of the light. Thus when only the 
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green component of the incident light is scat- 
tered, the red component passes through the 
liquid crystal panel, and a red image can be 
obtained on the screen. 

The most effective scattering wavelength 
dependends on the liquid crystal material and 
the panel gap, therefore one NCPT panel can 
produce a single color image. A double layer 
NCPT projector can produce a multicolor 
image”). The structure of the double layer 
NCPT liquid crystal projection display is shown 
in Fig. 15. When the liquid crystal is in the 
cholesteric phase, the green light scattering 
panel (layer R: upper panel) allows only the 
red component of the incident light to pass 
through, and scatters the other color compo- 
nents. The red light scattering panel (layer G: 
lower panel) allows only the green component 
to pass through, and scatters the other color 
components. With this system, red, green, white 
and black can be displayed. The scattered 
components are dependent on both the selective 
Bragg reflection caused by the liquid crystal’s 
helical pitch, and the retardation of the liquid 
crystal. This system can produce four colors, 
However, by adding a blue light scattering panel, 
an eight-color multicolor projector is possible, 
and dithering technology can produce a full 
color image. 


3.2 STN projection displays 

The STN is the most popular type of LCD. 
Many office automation devices use STN LCDs 
with an A4 panel size. Thus it is quite natural 
to use STN LCDs for electronic overhead 
projectors. Several types of STN electronic 
overhead projectors are currently available from 
several companies. These STN projectors offer 
640 x 480 pixels or 640 x 400 pixels. However, 
the size in panel temperature governs the charac- 
teristics of these projectors and it is quite 
difficult to produce 1200 x 800 pixel STN 
projectors. 

Conner?” proposed a_ high resolution, 
1200 x 1000 pixel multicolor electronic over- 
head projector using STN technology. Figure 16 
shows the structure of the high resolution STN 
projection system. Two STN panels are stacked 
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Fig. 16—Double-layer STN fine resolution panel for 
projection display. 
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Fig. 17—Structure of multicolor STN projection display. 


with each panel’s pixels being compensated as 
shown in Fig. 16. Each panel is driven by a 
1/256 duty ratio, which has enough drive 
window with a wide temperature range, giving a 
total resolution of 1 200 x 1 000 pixels. A color 
image is produced by birefringence of the 
liquid crystal — in principle, STN LCDs have a 
birefringence effect — and incident white light 
turns colored by this effect. Using three STN 
layers, Conner realized a multicolor projection 
display as shown in Fig. 17. The first layer is a 
red panel, which turns incident white light to 
red. The second layer is a cyan panel, which 
turns incident light to cyan. The third layer is a 
yellow panel, which turns the light to yellow. 
This subtracting system achieves a multicolor 
projection image. 

Although these clever techniques have 
been used, heat absorption of the STN panel 


FUJITSU Sci. Tech. J., 28, 1, (March 1992) 


A, Mochizuki et al.: Liquid Crystal Large Screen Projection Displays 


limits the information content. The limitations 
of the STN projection system are thus low 
information content, image deterioration due 
to higher temperatures and poor luminance of 
the projection screen. 


4. Conclusion 

The liquid crystal projection display is an 
effective system for producing a large image. 
The optically addressed projection system 
offers high information content and high resolu- 
tion, but is so complicated that it is bulky, 
heavy and expensive. This system may be 
suitable for special applications such as a huge 
control panel for a command control. 

The electrically addressed liquid crystal 
projector is generally compact and inexpensive. 
The active matrix TN type is particularly useful 
for video rate screen images. However, the 
active matrix projector requires its own optical 
system, so for electronic overhead projectors, 
the simple matrix LCD is more suitable. The 
STN projector is currently the most popular. 
However, this display needs two _ polarizers 
which reduce screen luminance absorb projec- 
tion light, and raise the temperature of the LCD. 

In contrast, the NCPT projector needs no 
polarizers, giving high screen luminance, and its 
bistability gives an extremely high information 
content, reaching five million pixels. This non- 
polarized projection display is suitable for 
displays requiring good screen luminance and 
high image quality. Although the NCPT pro- 
jector can display only still images, the high 
resolution can produce gray scales. 

The NCPT projector offers the same high 
luminance, resolution and information content 
image as a laser beam addressed projector, 
yet is as simple, compact and inexpensive as a 
conventional STN projector. It can also produce 
full color images. Therefore, the NCPT projector 
is an effective presentation device. 
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Semiconductor Thin Films 
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This paper reviews recent advances in the application of transmission electron microscopy to 


the evaluation of structures, interfaces and defects in compound semiconductor thin films. 
It has been shown microscopically for the first time that the growth of (GaAs),,/(GaP),, 
strained layer superlattices by atomic layer epitaxy proceeds in a precisely layer-by-layer 
manner. Detailed characterizations of interfaces and defects in AlAs, GaAs and GaP on Si, 
also GaAsSb and Fe-doped InP on InP are described; later, the mechanisms which generate 
such defects and techniques for their elimination are discussed. Studies of atomic order- 
ing in InGaP and InGaAs are presented, with reference to the mechanisms which generate 
ordered structures and their influence on the properties of the material. 


1. Introduction 

IlIl-V compound semiconductor thin films 
such as GaAlAs/GaAs and InGaAs(P)/InP grown 
by the advanced growth techniques of metalor- 
ganic vapor phase epitaxy (MOVPE), molecular 
beam epitaxy (MBE) and atomic layer epitaxy 
(ALE), are widely used as materials for 
optical devices, electron devices, and their inte- 
grated circuits. In order to ensure the best 
possible device characteristics and high device 
performance, one must achieve precise control 
of the layered structures, i.e. control of layer 
thicknesses and uniformity, and abruptness of 
the hetero-interface, at the atomic level. To 
achieve high device reliability, it is also necessary 
to solve certain materials issues in the thin films, 
i.e. generation of defects in the crystal, and 
thermal instability of the crystal. 

To confirm controllability in the growth of 
layered structures and to solve the materials 
issues described above, one must perform 
structural evaluation of these materials micro- 
scopically or even at the atomic level. Transmis- 
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sion electron microscopy (TEM) |i.e. both 
conventional TEM and _ high-resolution TEM 
(HRTEM)}, is a powerful structural evaluation 
technique which can provide a considerable 
degree of information regarding structures, 
interfaces and defects in the thin films. 

In this paper, the current state of our appli- 
cation of TEM to the evaluation of structures, 
interfaces and defects (including those due to 
thermal instability) in a variety of III-V semi- 
conductor thin films is systematically reviewed. 


2. TEM 

In this chapter, three major techniques in- 
volved in TEM are described. They are specimen 
preparation techniques, observation techniques, 
and computer simulation techniques of HRTEM 
images. 
2.1 Specimen preparation techniques”? 

For TEM observation, preparation of speci- 
mens is a crucial factor. A good TEM specimen 
must: 
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1) be thin enough; 

2) have a flat surface; 

3) be undamaged; 

4) have no surface contamination; 

5) have no bowing. 

In particular, the thickness of the sample is 
most important since the incident electron beam 
should pass through the specimen. Whether or 
not the sample is transparent to the electron 
beam also depends on the voltage used to 
accelerate the electrons, and the material. When 
200 kV TEM is used, the maximum thicknesses 
for Si and GaAs are | wm and 0.5 um, respec- 
tively. However, the optimum thickness for con- 
ventional TEM is 0.1-0.2 ym, and for HRTEM, 
the thickness should be less than 50 nm. 

Specimen preparation techniques can be 
classified into two types, chemical etching and 
ion milling. In chemical etching, Br,/CH,OH 
solution is mostly used for II-V compounds, 
after polishing the sample down to about 
100 pm. In the case of layered structures such as 
InGaAsP/InP double-heterostructures, preferen- 
tial etching techniques can be used, in which the 
substrates are removed preferentially leaving 
only the epitaxial layers. This technique is 
useful for the preparation of plan-view TEM 
specimens. To prepare cross-sectional TEM 
specimens, we must thin the sample by the other 
technique, ion milling. 

Preparation consists of the following steps: 
1) Sticking the sample epitaxial wafers face-to- 

face with epoxy glue; 

2) Cutting the sample processed in step 1) into 
cross-section slices with a thickness of about 
0.5 mm; 

3) Polishing the specimen from both sides 
down to about 50 ym; 

4) Mounting the specimen onto a single-hole 
TEM mesh; 

5) Reducing the specimen to final thickness by 
ion milling from both sides. 


2.2 Observation techniques®? 

In this study, we used a Topcon EM-002B 
ultra-high resolution analytical electron micro- 
scope equipped with an energy dispersive X-ray 
spectrometer (EDX) and electron’ energy 


loss spectrometer (EELS). We applied both 
conventional TEM techniques and HRTEM 
techniques. 

2.2.1 Conventional TEM techniques 

In conventional TEM, there are two main 
methods used: bright-field and  dark-field 
imaging. A bright-field image is obtained from 
the transmitted electron beam, whereas a dark- 
field image is from one of the Bragg reflected 
beams. 

Thus, these images are complimentary. In 
these observations, several diffraction conditions 
are used: 

1) Two-beam condition 
One particular reflection is excited. 
2) Multi-beam condition 

More than two reflections are excited. 
3) Quasi-kinematic condition 

No reflections are excited. 

2.2.2 HRTEM 

The HRTEM technique provides us with 
precise atomic level contrast; the contrast is 
“phase-contrast” which is produced as a result 
of multiple-scattering of transmitted electrons in 
the specimen”). In this situation, a zone-axis 
diffraction condition is used where the direction 
of the incideni electron beam is set exactly 
normal to the crystallographic plane of the spec- 
imen (in this condition, all equivalent reflected 
beams are equally excited), The actual image is 
made by allowing many beams through a circu- 
lar objective aperture at high magnification 
(typically x 300000). The HRTEM image 
consists of lattice-fringes corresponding to low- 
index spacings of the materials (200, 111, 220, 
400,...); thus this is called a “lattice image’’ or 
“structure image’’. 


2.3 Computer simulation of HRTEM images 

The HRTEM images described in subsection 
2.2.2 strongly depend on the defocus value of 
the objective lens, the thickness of the specimen, 
the optical parameters of the electron micro- 
scope, etc. Therefore, computer simulation is 
required for a complete understanding of the 
HRTEM images. The most important method 
for computing such images is the multi-slice 
method®). In_ this study, we have also applied 
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this method to simulate the HRTEM images of 
various types of III-V superlattices. 


3. Evaluation of layered structures 

In order to ensure the best possible device 
characteristics and high device performance, one 
must achieve precise control of the layered 
structure, i.e. atomic-level control of layer 
thicknesses and uniformity, and abruptness of 
the hetero-interface. In this chapter, recent 
results on structural evaluation of GaAs/GaP 
strained layer superlattices grown by atomic 
layer epitaxy, which allows layer-by-layer 
growth of compound semiconductors, is de- 
scribed. 


3.1 GaAs/GaP strained layer superlattices grown 

by ALE 

ALE is a very attractive technique in which 
crystals can be grown layer-by-layer through a 
self-limiting mechanism!). It can offer great 
benefits in fabricating various III-V heterostruc- 
tures with extremely high uniformity in thick- 
ness, excellent abruptness of the hetero-interface 
and novel artificial structures such as monolayer 
superlattices, quantum wires and boxes. We have 
previously reported the development of a new 
ALE technique called pulsed-jet epitaxy (PJE) 
which can achieve layer-by-layer growth of III-V 
compounds with extremely good _ control- 
lability®"1!). Moreover, we have described 
successful growth of (GaAs), /(GaP), (m, n 
=1, 2, 3,...) strained layer superlattices by 
PJE'”). Structural analyses using X-ray diffrac- 
tion'?)? and Raman scattering spectroscopy!» 
have indicated that growth of the superlattices 
was completely controllable at an atomic level 
even for the case of a (GaAs), /(GaP); mono- 
layer superlattice. Furthermore, luminescence 
and reflectance measurements strongly suggested 
that the (GaAs),/(GaP), superlattice has a 
direct bandgap'?”!*. However, no one has ever 
reported any direct evidence for layer-by-layer 
growth of III-V compound semiconductors by 
ALE. 

In this section, we describe for the first time 
a detailed TEM evaluation of various 
(GaAs) /(GaP), (m,n = 1,2,3,...) superlattices 
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grown on (001) GaAs substrates by ALE. In this 
work, (GaAs), /(GaP),, (GaAs)3/(GaP),; and 
(GaAs)4/(GaP), strained layer superlattices were 
examined by TEM. The (GaAs),,, /(GaP),, super- 
lattices were grown on (001) GaAs (just) sub- 
strates by ALE (PJE). Both GaAs and GaP were 
grown at 500°C, where both materials can be 
grown layer-by-layer. Trimethylgallium (TMGa), 
arsine (AsH; ) and phosphine (PH; ) were used as 
source gases. The alternating pulses of TMGa 
and arsine (phosphine) were separated by a 
hydrogen purge pulse to prevent the reactants 
from mixing. A vapor pressure of 20 Torr was 
maintained during the growth {growth proce- 
dures are described in detail elsewhere!?)}. The 
thicknesses of these superlattices are approxi- 
mately 50 nm. In this work, we have performed 
cross-sectional TEM observations. 

3.1.1 Structural evaluation of (GaAs), / 

(GaP), monolayer superlattices 

Figure 1 shows a_ typical transmission 
electron diffraction (TED) pattern obtained 
from the (110) cross-section of such a (GaAs), / 
(GaP), monolayer superlattice. Since the 
superlattice is 50nm thick and the selected 
area aperture used in this experiment was 
150nm in dia, this TED pattern is mostly 
from the GaAs substrate. However, super- 
structure spots are observed at positions indexed 


Fig. 1—A TED pattern from a (110) cross-section of a 
(GaAs), /(GaP),; monolayer superlattice. 
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Fig. 2—A (110) cross-sectional HRTEM image of a 
(GaAs), /(GaP),; monolayer superlattices. 


as 001, 110, 112,... (see diffraction spots 


indicated by arrows.). This particular set of 


superstructure spots is associated with the 
(GaAs), /(GaP), monolayer superlattice. In 
order to clarify if such superlattices are gener- 
ated in real space, HRTEM observation was 
carried out. Figure 2 is a (110) cross-section 
HRTEM image of the superlattice. In this image, 
one bright spot corresponds to two adjacent 
atoms, i.e. Ga-As or Ga-P along the <110> 
direction. The interface between the epitaxial 
layer and the GaAs substrate is indicated by an 
arrow. Computer simulation of (GaAs), /(GaP), 
confirms that in Fig. 2, the bright spots with 
stronger and weaker intensities correspond to 
Ga-As and Ga-P (specimen thickness = 2 nm, 
amount of defocus =38 nm), respectively (see 
also arrows in the epitaxial layer). Note that the 
hetero-interface is extremely abrupt with no 
imperfections and no alloyed region, and that 
the (GaAs), /(GaP), monolayer superlattice is 
formed uniformly in the entire epitaxial layer 
region. This is the first direct evidence for layer- 
by-layer growth of III]-V compound semiconduc- 
tors by ALE. 

3.1.2 Structural evaluation of (GaAs); / 

(GaP); and (GaAs), /(GaP ), 

We have also evaluated multiple-layer 
superlattices of (GaAs);/(GaP),; and (GaAs), / 
(GaP),. A TED pattern from the (110) cross- 
section of (GaAs)3/(GaP), is shown in Fig. 3. 
A sharp satellite-spot pair is present around the 
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Fig. 3—A TED pattern from a (110) cross-section of the 
(GaAs)3/GaP)3 superlattice. 


Fig. 4—A (110) cross-sectional HRTEM image of a 
(GaAs)3/(GaP)3 superlattice. 


000 (direct) spot (see spots around the 000 
indicated by arrows), which corresponds to 
the periodicity of (GaAs)3;/(GaP),. The spacing 
between the satellite and 000 spots is approxi- 
mately one sixth of that between the 002 and 
000 spots. This result indicates that the super- 
lattice has periodicity of 6 atomic layers. In 
the higher order spot region, the fundamental 
spot splits into two individual spots as shown in 
the upper part of the TED pattern. In this 
region, spots for 004 reflections of GaAs and 
GaP are observed with very narrow spacing. 
Since the midpoint between two satellite spots 


FUJITSU Sci. Tech. J., 28, 1, (March 1992) 


O. Ueda: Recent Advances in the Application of TEM to the Evaluation of Ill-V Compound... 


Eg 


ee heeeees 
ANNA GaAs 


mw GaAs 


2nm 


Fig. 5—A (110) cross-sectional HRTEM image of a 
(GaAs)q/(GaP) 2 superlattice. 


nearly coincides with that between the 004 
spots, the periodicity of the superlattice can be 
treated as that of a GaAsP mixed crystal. 

Figures 4 and 5 show (110) cross-section 
HRTEM images of (GaAs)3;/(GaP),; and 
(GaAs)4/(GaP), superlattices. These images also 
confirm that bright spots with stronger and 
weaker intensities correspond to Ga-As and 
Ga-P, respectively. The locations of the GaAs 
and GaP layers are indicated in each figure. 
GaAs and GaP layers of exactly three mono- 
layers have grown uniformly throughout the 
epitaxial layer. In both cases, neither alloyed 
regions nor oxidized regions were found at the 
hetero-interface. Also note that at each hetero- 
interface in the superlattice, {111} planes 
deform to an ‘‘S’’-shape. This is probably due 
to lattice distortion of the crystal at each inter- 
face due to lattice-mismatch. Computer simula- 
tion is required to confirm this. 


4. Interface structure and defects 

Defect generation is one of the important 
materials issues in II-V semiconductor thin 
film growth. Defects introduced during growth 
can be classified into two types: interface 
defects and bulk defects. 

Defects belonging to the first type are stack- 
ing faults, V-shaped dislocations, dislocation 
clusters, microtwins, inclusions, and misfit 
dislocations. Interface defects due to contamina- 
tion, thermal damage of the substrate, local 
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b) V-shaped dislocations 


a) Stacking fault : 
) Steyits and dislocation clusters 


d) Inclusions 


c) Microtwins 


Fig. 6—Schematic diagrams of interface defects due to 
contamination, thermal damage of the substrate, 
or local segregation. 


a) A misfit dislocation b) A misfit dislocation 
caused by a threading caused by a threading 
dislocation gliding out dislocation bending 
from the substrate out 


c) Misfit dislocations in 
the form of half-loops 


NY; 


Lj 


d) A microtwin 


Fig. 7—Schematic diagrams of interface defects due to 
lattice mismatch. 


segregation at the interface, and those due to 
lattice-mismatch are schematically shown in 
Figs. 6 and 7. 

The stacking faults shown in Fig. 6a) are 
generated from the hetero-interface, extending 
along the four equivalent {111} planes, forming 
a stacking fault tetrahedron. In MBE-grown 
crystals, they often correspond to surface 
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defects (see section 4.3). Pit-pairs are often 
observed on the etched surface of thin films. 
In the TEM image of the area which includes a 
pit-pair, small dislocation segments are observed. 
From the shape of the dislocations, this dis- 
location pair is assumed to comprise a V-shaped 
dislocation {see Fig. 6b)}. Occasionally, pref- 
erential etching takes place at such V-shaped 
dislocations! *). Microtwins are generated at the 
hetero-interface through relaxation of stress, 
extending along the {111} planes, as shown in 
Fig. 6c). They occasionally terminate with {211} 
type boundaries in mid-layer. Inclusions are 
accidentally introduced during liquid phase 
epitaxy (LPE) growth by incorporation of 
remaining droplets into the solid solution when 
the Ga- or In-rich melt is wiped off from the 
substrate'® {see Fig. 6b)}. Since the first three 
kinds of defect described above are eliminated 
by clearing the substrate surface prior to growth 
and/or the protection of the substrate against 
evaporation of substrate atoms, they are pre- 
sumably generated by contamination of, and/or 
thermal damage to the substrate. 

Thus, control of the hetero-interface is 
extremely important for the elimination of 
defects. 

If there is a difference in lattice constant 
between the epitaxial layer and the substrate, 
misfit dislocations are generated when the film 
thickness exceeds the critical thiciness f, 
which depends on the structure and elastic 
properties of the material, and the type of dis- 
locations. There are several possible mechanisms 
for the generation of misfit dislocations as 
shown in Figs. 7a) to c). First, threading dis- 
locations propagating from the substrate may 
glide out {see Fig. 7a)} or bend out{see Fig. 7b)} 
to the edge of the wafer. It is also expected that 
two dislocations with the same Burgers vectors 
in the substrate can terminate as one misfit 
dislocation at the hetero-interface, forming a 
half-loop {see right figure of Fig. 7c)}. 

If the dislocation density in the substrate 
is low, other mechanisms which may occur are: 
1) nucleation of micro-half-loops at the sur- 

face; 

2) expansion of loops down to the interface; 


a) Precipitates b) Precipitates accompanying 


punched-out dislocation 


c) Faulted and/or unfaulted d) Faulted and/or unfaulted 
loops due to local loops due to modulated 
segregation of dopant structure generation 
atoms 


Fig. 8—Schematic diagrams of bulk defects due to local 
segregation of dopant atoms or native point 
defects. 


3) gliding out of the threading segments to the 
edge {left figure of Fig. 7c)}'7. 

In highly mismatched heteroepitaxy such 
as GaAs/Si, microtwins {see Fig. 7d)} and/or 
stacking faults are very often generated, partic- 
ularly during the initial stage of growth!®). 
They propagate along the four equivalent {111} 
planes. However, they often annihilate in mid- 
layer, i.e. terminating in {211} type boundaries 
for microtwins!?*®) and coalescence of two 
stacking faults on different planes!®?. 

Figures 8a) to d) show typical features of 
bulk defects. In most cases, they are generated 
in heavily-doped crystals as shown in Figs. 8a) 
to c). When supersaturated dopants are present 
in the crystal, they tend to segregate locally, 
either during growth or during the cooling 
process after growth, forming precipitates 
{see Fig. 8a)}. When there is large lattice- 
mismatch between the precipitates and the 
matrix, dislocation loops are often punched out 
by the relaxation of stress {see Fig. 8b)}. These 
precipitates are often observed in heavily doped 
LPE-grown InGaAsP crystals on InP?!. In the 
formation of precipitates, excess native point 
defects can be introduced. They may condense 
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at some nucleation centers to form dislocation 
loops (either faulted or unfaulted) as shown in 
Fig. 8c). Interstitial type dislocation loops are 
often observed in heavily doped LPE-grown 
GaAlAs crystals with Ge??). On the other hand, 
dislocation loops are present even in undoped 
crystals. We have found both vacancy type and 
interstitial type loops in undoped LPE-InGaP 
crystals*>), Since modulated structures (see 
Chap. 5) are generated in the crystal, these loops 
are presumably caused by the condensation of 
excess point defects induced as byproducts of 
modulated structures. 

In this chapter, to illustrate some of the 
topics described above, we present recent 
results on the TEM evaluation of interfaces 
and defects in highly mismatched AlAs, GaAs 
and GaP grown on (001)Si. Analyses of interface 
defects in GaAsSb on (001)InP heterostructures 
and bulk defects in heavily Fe-doped InP 
crystals are also described. 


4.1 Interfaces and defects in AlAs on Si(001) 

grown by ALE 

Highly mismatched heterostructures such as 
GaAs/Si?*?®, InP/Si?®, and SiGe/Si?!) are 
very attractive materials for use in opto-electronic 
integrated circuits, solar cells and heterobipolar 
transistors. The main problems in these systems 
are the differences in lattice constants and in 
thermal expansion coefficients between the 
substrate and epitaxial layers. At the initial 
stage of crystal growth, although the lattice 
mismatch may be relaxed by the presence of 
misfit dislocations at the hetero-interface, three- 
dimensional island growth often takes place, 
generating in the crystal a high density of de- 
fects such as threading dislocations, stacking 
faults and microtwins. In GaAs/Si the defect 
density can be reduced to about 10°/cm? by 
using the standard two-step growth method?”, 
applying in situ cyclic annealing?®»?®, and 
introducing strained layer superlattices? ’»?®), 
However, during the cooling process after 
growth the difference in the thermal expansion 
coefficients causes additional stress in the 
system, thus enhancing the multiplication of 
dislocations via glide motion?®’. Therefore, in 
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order to achieve a drastic reduction of the defect 

density, the following approaches are required: 

1) Operating a two-dimensional growth mode 
from the initial stage; 

2) Fabricating the heterostructure above the 
buffer layer at lower temperatures. 

We have previously reported that GaAs and 
AlAs single crystals can be directly grown on Si 
substrates at 500°C using ALE?”)35), From 
Auger electron spectroscopy and Raman scatter- 
ing spectroscopy, it has been shown that growth 
by ALE on Si substrates is initially three-dimen- 
sional but that it changes to layer-by-layer 
growth at a relatively early stage. We have also 
found that coverage improves remarkably when 
initiating the growth with AlAs instead of GaAs. 
However, the atomic structure of the hetero- 
interface, the critical thickness for two- 
dimensional growth and the defect structures in 
the GaAs and AIAs layers have not yet been 
clarified. In this section, we describe for the first 
time a detailed structural evaluation of ALE- 
grown AlAs layers on (001)Si substrates by 
TEM. 

Thin AlAs layers were grown on a (001)Si 
substrate (tilted 2° toward the <110 > direc- 
tion) at 500 °C. Prior to growth, the substrates 
were heated to 1000°C for 20 min in H, to 
evaporate the oxidized surface layer. Since 
exposure to air heavily damages the surface of 
AlAs, a 10-atom layer of GaAs was deposited as 
protection. We also grew GaAs (3 wm)/AlAs/Si 
heterostructures to evaluate defect generation 
in thick GaAs layers. In this structure, GaAs 
layers were grown by MOVPE. Arsine (AsH3) 
and trimethylgallium (TMGa) or trimethylalumi- 
num (TMAI) were used as source gases and 
sprayed on the substrates as pulsed-jets®®. 

4.1.1 Atomic structure of the hetero-inter- 

face 

First, the hetero-interface between the AlAs 
layer and the Si substrate was evaluated at the 
atomic level. As a reference, we also examined 
a GaAs/Si heterostructure grown by standard 
two-step MOVPE. Figure 9 shows a typical 
(110) cross-sectional HRTEM image of the 
reference sample. The hetero-interface is not 
abrupt, but highly uneven. The typical interface 
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Fig. 9—A (110) cross-sectional HRTEM image of a 
hetero-interface of the GaAs/(001)Si 
heterostructure grown by standard two-step 
MOVPE. 


roughness is in the range of 5-10 atomic layers. 

An amorphous-like contrast is also observed 

locally. These features are presumably caused by 

three-dimensional growth and the reaction 
between the GaAs buffer layer and the substrate 
during the annealing process (2nd step) of the 
buffer layer growth. On the other hand, the 

AlAs/(001)Si interface is found to be extremely 

sharp |see Fig. 10a)! compared to that of 

GaAs/Si. Also note that atomic steps are clearly 

observed {see arrows in Fig. 10a)} and that 

no amorphous or alloyed regions are found. 

Figure 10b) is a high magnification image of part 

of Fig. 10a). 

At the hetero-interface, two distinct features 
are found: 

1) The (111) fringes in both Si and AlAs 
regions do not coincide with each other as 
indicated by A,-A3 and B,-B;; 

2) An extra layer with anomalous fringe 
contrast, which is not observed at the GaAs/ 
Si interface, appears (see the layer indicated 
by X). In this layer, the spacing of the 
brighter spots along the <110> direction 
(from left to right) is half of that in the Si 
and GaAs regions. 

As previously reported for ALE of GaAs on 

a GaAs substrate, one atomic layer of GaAs is 

grown in one growth cycle®®’. However, for 

AlAs on GaAs, one cycle allows the growth 
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b) A high magnification image of the part of 
the HRTEM image a) 


Fig. 10—(110) cross-sectional HRTEM images of a 
hetero-interface of the AlAs/(001)Si hetero- 
structures grown by ALE. 


of two atomic layers of AlAs. The interpretation 
of the difference in growth rate between GaAs 
and AlAs is as follows?®): during GaAs deposi- 
tion, one Ga atom for every As atom is adsorbed 
on the (001 )GaAs surface, whereas, during AlAs 
deposition, the (001)Al face touches the GaAs, 
resulting in adsorption of two Al atoms per As 
atom?”). This is due to a very small lattice 
mismatch (1.3 percent) between the Al and 
GaAs. Therefore, the extra layer is presumably 
one atomic layer of metallic (001 )Al. We assume 
that a similar mechanism operates during the 
growth of AlAs on Si. To confirm this, 
computer simulation is required. 
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4.1.2 Initial stage of growth 

Figure 11 shows a _ bright-field plan-view 
TEM image of a thin AlAs layer on (001)Si 
taken under zone-axis illumination. In _ this 
sample, a 6-atom layer of AlAs is covered by a 
10-atom layer of GaAs to protect the AlAs 
surface. Since the plan-view specimen includes 
Si, AlAs and GaAs, Moiré fringes are observed 
as a result of the interference between the trans- 
mitted wave and the doubly reflected waves at 
the AlAs/Si interface. As shown in Fig. 11, 
numerous patchy regions with Moiré fringes 
in two directions of nearly <110> and 
<110> are formed. These regions are thought 
to be three-dimensional islands of AlAs on the 
substrate. They are 5-20 nm in dia. Rotational 
misalignment of the Moiré fringes is observed, 
indicating that the islands rotate slightly at the 
initial stage of growth to minimize the strain 
energy. Cross-sectional observation was also 
carried out on the GaAs/AIAs/(001)Si hetero- 
structures. Figures 12a) and b) show cross- 
sectional HRTEM images of AlAs/(001)Si at an 
initial stage of growth. In Figs. 12a) and b), 
6- and 26-atom layers of AlAs were grown, 
respectively (in both cases, 10-atom layers of 
GaAs were grown to protect the surface of the 
AlAs). In Fig. 12a), islands of 5-20 nm in dia are 


—— 
100 nm 


Fig. 11—A bright-field plan-view TEM image of 
AlAs/(001)Si heterostructures taken under 
zone-axis illumination. 
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observed uniformly in the entire region. Thus, 
one can conclude that a_ three-dimensional 
growth mode is dominant at the initial stage of 
growth of AlAs. However, as the AlAs layer 
grows thicker, the growth mode becomes more 
two-dimensional, forming a rather smoother 
surface {see Fig. 12b)}. 

Figures 13a) and b) show cross-sectional 
dark-field TEM images of GaAs/AlAs/(001)Si 
heterostructures with 36- and 80-atom layers of 
AlAs, respectively, capped by a 3 um thick 
GaAs layer, obtained under the 002 reflection. 
Since the intensity of the 002 reflection is very 
low in GaAs and Si, the AlAs layer appears 


brighter. In Fig. 13a), the interface between 
AlAs and GaAs is nearly smooth although 


regions with dark contrast are often observed in 
the AlAs layer. When the AlAs was grown 
thicker {see Fig. 13b)}, the hetero-interface was 
found to be just as smooth. These results can 
lead us to conclude that the critical thickness for 
achieving two-dimensional growth in AlAs on 
(001 )Si is approximately 36 atomic layers. 

Next, in order to compare the growth mode 
of AlAs on (001)Si with that of GaAs on 
(001 )Si, we have grown heterostructures consist- 
ing of three stacked buffer layers, the first AlAs 
(a 70-atom layer), the second GaAs (a 70-atom 
layer), and the third AlAs (a 70-atom layer), 
capped by a GaAs (3 um) layer. A cross-sec- 
tional TEM image of this structure is shown in 
Fig. 14. The interface between the first AlAs 
layer and the second (GaAs) layer (the region 
with darker contrast) is smooth, whereas the 
interface between the second (GaAs) and the 
third (AlAs) layer is very uneven. However, the 
interface smoothness between the third layer 
and the thick, upper GaAs layer has improved. 
These results indicate that the critical thickness 
for a change from three-dimensional to two- 
dimensional growth in GaAs on (001)Si is much 
larger than that for AlAs on (001)Si. As 
described in section 3.1, two atomic layers are 
grown per ALE cycle for AlAs, achieving 
coverage twice as good as for GaAs (the bond 
strength of Al-Si, which is greater than that of 
Ga-Si, could also enhance the coverage)”. 
Thus, the difference in the critical thickness 
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Fig. 12—(110) cross-sectional HRTEM images illustrating 
the initial stage of growth of AlAs layers on 
(001)Si substrates. In both cases, a 10-atom 
layer of GaAs was grown to protect the surface 
of the AlAs layer. 


50 nm 
Fig. 14—A (110) cross-sectional dark-field TEM image of 
GaAs with AlAs/GaAs/AIAs buffer layers 
grown on (001)Si substrate taken under the 
002 reflection. 
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b) An 80-atom layer of AlAs 22"™ 


Fig. 13—(110) cross-sectional dark-field TEM images of 
GaAs (2 um)/AIlAs/(001)Si heterostructures 
grown on (001)Si substrates taken under the 
002 reflection. 


between AlAs and GaAs could be explained by 
the improvement in coverage in the case of AlAs 
deposition. 

4.1.3 Defects in the AlAs layer 

Two types of defect are present in the AlAs 
layer, stacking faults and microtwins. They are 
generated at the initial stage of growth as shown 
in Figs 12a) and b). In both micrographs, stack- 
ing faults and microtwins are indicated by S and 
M, respectively. Self-annihilation of these 
defects is observed at regions indicated by A in 
Fig. 12b). As the AlAs grows thicker, the 
density of these defects increases as shown in 
Fig. 13b). Generation of such a high density of 
defects in AlAs may be associated with a low 
stacking fault energy and the softness of the 
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material. However, surprisingly, the defects 
annihilate at the interface between the AlAs and 
the thick upper GaAs layer. From the images 
shown in Figs. 12b) and 13b), possible models 
for the self-annihilation of stacking faults and 
microtwins are as follows: 

1) Stacking faults propagating on different 
{111} planes annihilate where they meet 
each other; 

2) Microtwins terminate on the {211} planes 
normal to the twin boundary at the inter- 
face between the AlAs and the thick upper 
GaAs layer to minimize the strain energy. 
The results described above lead us to the 

conclusion that using an AlAs layer grown by 

atomic layer epitaxy as a buffer layer is a very 
promising way to achieve high-quality GaAs/Si 
heterostructures with a low defect density. 


4.2 Anti-phase boundaries in GaAs and GaP 

crystals grown on (001)Si by MOVPE 

In highly mismatched hetero-epitaxy, apart 
from the problems of lattice mismatch and 
thermal mismatch, the generation of anti-phase 
domains (APDs)?®? in the growth of polar semi- 
conductors such as GaAs on (001) nonpolar 
semiconductors such as Si, has_ received 
considerable attention in the past. It is now 
recognized that they can be controlled in 
practice by using substrates with an off-angle of 
more than 2° off Si®®. However, the exact 
nature and mechanism for their generation (and 
annihilation, if any) is not yet fully understood. 
At present, two unique TEM techniques can be 
employed for surveying APDs. One is the 
convergent beam electron diffraction (CBED) 
approach*®”). The other is the +002 dark-field 
imaging technique proposed by Kuan et al*!). 
The latter method is a very simple one: we have 
also employed this method in our experiments. 
The technique enables us to produce distinct 
dark and bright contrast on both sides of the 
anti-phase boundary (APB), although the 
contrast is strongly related to the foil thickness. 
In this section, we describe for the first time a 
detailed TEM investigation of APDs and present 
direct evidence for their self-annihilation. 

Crystals were grown by MOVPE. Two types 
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of heterostructure are examined in this study. 

Heterostructures I consist of a GaP buffer layer 

(0.1 mm), followed by strained layer super- 

lattices (SLSs) of GaP/GaAsp.5P)., and 

GaAso,.5Po,;/GaAs (5 periods for each SLS and 

20 nm for each layer), and finally a thick GaAs 

layer (4mm). The growth temperatures for the 

GaP, SLSs and GaAs were 900 °C, 750 °C and 

750 °C, respectively. TMGa, AsH, and PH, 

were used as reactants {the growth procedure 

is described in detail elsewhere*?”**)|. Hetero- 
structures I] are grown by standard two-step 

MOVPE. They consist of a GaAs buffer layer 

(12.5 nm) grown first at 400°C, and a thick 

GaAs layer (3 um) grown at 650°C. In both 

cases, the heterostructures were grown on Si 

substrates with the following orientations: 

Types Ia and Ila were grown on (001) +0.5° Si 

substrates and Types Ib and IIb on (001)2° off 

toward the [110] direction. 

4.2.1 Evaluation of Types Ia and Ila 

structure 

First, plan-view and cross-section TEM ob- 
servations were performed for Type la and Ila 
structures. The plan-view observations reveal 
uniformly distributed domain-like regions of the 
type shown in Fig. 15. We conclude that these 
regions correspond to APDs for the following 
reasons: 

1) No extra spots for small-angle tilt boundaries 
are found in the electron diffraction pattern 
from the area containing the domain-like 
regions. 

2) A part of the boundary has fringe contrast, 
but this is not due to a stacking fault since 
the fringe is not asymmetric in the dark- 
field image (stacking faults exhibit asym- 
metric fringe contrast). 

3) In electron diffraction patterns from areas 
which include the domain-like regions in the 
(110) and the (110) cross-sections, no twin- 
related spots are observed. 

4) From cross-sectional observations, the con- 
trast of the domain itself in a dark field 
image is reversed when the reflection vector 
changes from 002 to 002 as illustrated in 
Fig. 18 (Fig. 18 is described in detail later). 


23 


O. Ueda: Recent Advances in the Application of TEM to the Evaluation of II-V Compound... 


Fig. 15—A (001) plan-view TEM image of APDs in the 
GaAs top layer of the Type Ia structure, taken 
under the < 001 > zone-axis illumination. 


OS am 


Fig. 16—A (001) plan-view bright-field TEM image of 
an APD in the GaAs top layer of the Type Ia 
structure, taken under the 220 reflection. 


The domain size is 0.5-2 wm. Each APD has 
a different geometry; they are found to be three- 
dimensional extended polygon-shaped structures 
consisting of sub-boundaries nearly normal to 
the (001) plane and inclined to the (001) 
plane. Figure 16 shows a typical plan-view 
bright-field TEM image of an APD taken under 
the 220 reflection. Since the sub-boundaries 
indicated by A-B and D-E do not display fringe 
contrast, they are presumably normal to the 
(001) plane. In view of their elongations, they 
are possibly (110)- and (100)-oriented bounda- 
ries, respectively. From plan-view observations, 
it is found that most of the APBs normal to the 
(O01) plane are (110) or (100) boundaries. The 
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Fig. 17—A (001) plan-view bright-field TEM image of an 
APD in the GaAs top layer of the Type Ila 
structure, taken under the 220 reflection. 


other boundaries indicated by B-C-D-C’-B and 
A-F'-E-F are curved and exhibit fringe contrast. 
This may be due to the phase-shift of the elec- 
tron beam when it enters the APD from the 
matrix crystal. Thus, they are inclined to the 
(O01) plane and are perhaps a mixture of tilted 
APBs indexed as (11m), such as (111) and (112) 
and (10n), such as (101) and (102). In each 
boundary, a number of dislocations can be seen. 
In this case, these boundary dislocations are 
invisible under the reflection of 400 or 040. 
Since most of the APBs are not lying on the 
{111} plane, these boundary dislocations cannot 
glide. Also in the case when stacking faults and 
microtwins are generated in the APD, they can- 
not be expected to extend to the region outside 
the APD, since there is a clear difference in 
atomic arrangement between the APD and the 
matrix crystal. From. stereoscopic observa- 
tions, it is found that the circuit indicated by 
A-B-C-D-E-F-A in Fig. 16 is lying below the top 
surface and that the circuit indicated by 
A-B-C’-D-E-F’-A is lying on the top surface. This 
finding suggests that the domain size becomes 
smaller during crystal growth with self-annihila- 
tion occurring as a final stage. 

In Type Ila structures, similar APDs to those 
described above were also observed as shown in 
Fig. 17, although not all of the APDs annihilate. 

Surprisingly, we have discovered direct 
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evidence for this speculation in the cross-section 
TEM image shown in Fig. 18. Figures 1 8a) and b) 
are dark-field images taken under weak 002 and 
002 reflections obtained from the (110) cross- 


section specimen of Type la. A high density of 
defects (mostly dislocations) is present in both 
the GaP buffer layer and the SLSs, since the 
crystal is grown on (001) +0.5° Si. In both 
figures, a domain extending from the interface 
between the Si substrate and the GaP buffer 


b) g = 002 
Fig. 18—(110) cross-sectional dark-field TEM images of 
the Type Ia structure. 
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layer into the GaAs is observed (see the domain 
denoted by APD). Note that in Fig. 18a) the 
domain contrast is bright whereas that of the 
matrix region is dark. On the other hand, the 
opposite result is obtained in Fig. 18b). Since 
the crystal structures on both sides of an APD 
are related to each other by 180° rotation about 
the <110> direction, the 002 reflection from 
the crystal on one side of an APB should be 
equivalent to that of the 002 reflection from the 
other side of the APB. We have calculated the 
amplitude of the transmitted beam and the 002 
and 002 beams in an (110) diffraction pattern 
from perfect GaP and GaAs crystals as a func- 
tion of foil thickness (see Fig. 19). It is found 
that due to the lack of twofold axes along the 
<110> direction in the zinc-blende structure, 
the amplitudes of the 002 and the 002 reflec- 
tions are not equal for most crystal thicknesses. 
Thus, one can conclude that the image intensi- 
ties from both sides of the APB must be differ- 
ent except for certain foil thickness. In Figs. 18a) 
and b), the upper part of the micrograph cor- 
responds to the edge of the TEM specimen. The 
thickness is almost constant in the direction 
parallel to the interface and gradually increases 
from the upper part (edge of the foil) to the 
lower part. The estimated thickness of the 
major part of the foil where the APD is clearly 
observed (region denoted by APD) is roughly 
30-50 nm. Based on this consideration, the 
domain observed in Figs. 18a) and b) is identi- 
fied to be an APD similar to that shown in 
Fig. 16. As can be seen in Figs. 18a) and b), 
the APD annihilates in the intermediate layer 
of GaAs by changing the orientation of sub- 
boundaries from normal to (001) to higher- 
index-planes (see also the region denoted by 
arrows). 

On the basis of these results, the following 
scenarios are derived for the self-annihilation 
mechanism of APDs during crystal growth: 

1) At the initial stage of crystal growth, three- 
dimensional island growth proceeds, at the 
step-edges. By the coalescence of islands 
with two different phases (i.e. matrix-phase 
and anti-phase), APDs are structurally intro- 
duced. The APDs may be originally co- 
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Fig. 19—Calculated thickness dependence of the 
amplitudes of the 000, 002, and 002 reflections 
from GaAs and GaP crystals. 


lumnar-shaped consisting predominantly of 
sub-boundaries nearly normal to the (001) 
plane such as (110) and (100). Alternately, 
the sub-boundaries may be inclined to the 
(001) plane from an early stage of crystal 
growth, in which case step 3) takes place 
immediately. 

2) When the thickness of the APD reaches a 
critical value, these sub-boundaries change 
their orientation to a higher-index plane 
inclined to the (O01) plane {such as the 
(112) plane, which is the orientation with 
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lowest energy*”?} , so as to minimize the 

total energy in the system. 

3) By the coalescence of APBs with opposite 
orientation, the APD finally annihilates itself 
at the intermediate layer. At present, the 
reason for the existence of (100) APBs at 
the initial stage of crystal growth is not 
clear. Since not all of the APDs annihilate, 
it is also difficult to estimate the critical 
thickness for step 2). 

4.2.2 Evaluation of Type Ib and IIb 

structures 

Plan-view observations of both Type Ib and 
IIb structures reveal no APDs similar to those 
shown in Figs. 16 and 17 in the surface region of 
the GaAs layer. 

Figures 20a) and b) are 002 and 002 dark- 
field TEM image from the (110) cross-section 
of a Type Ib structure, illustrating the interface 
between the GaP buffer layer and Si substrate. 
Hillock-like micro-domains are observed in the 
GaP buffer layer near the interface (see arrows). 
The variation in the diffraction contrast of these 
micro-domains is very similar to that in the 
APDs shown in Fig. 4. Moreover, in this case, 
reasons | )-4) described previously (see the begin- 
ning of subsection 4.2.1), are satisfied. There- 
fore, it is concluded that they are also APDs, 
although very small, with diameters of 
10-20 nm. Provided that there are a certain 
number of monolayer steps on the surface of a 
(001) 2° off Si substrate (there is still contro- 


versy about the existence of double-layer steps 
on the vicinal surface of Si), the frequent pres- 


ence of APDs is well explained by the average 
spacing between monolayer steps, i.e. approx- 
imately 4 nm. In the Type IIb structures, hill- 
ock-like APDs similar to those in Type Ib were 
also observed. 

Regarding the self-annihilation of APDs 
at the initial stage of growth on tilted substrates, 
there have been two competing models, one by 
Pukite and Cohen**) and the other by Kawabe 
et al*®. Although the sub-boundaries of the 
hillock-like APDs shown in Fig. 20 are lying 
nearly on {211} planes, our findings strongly 
support the latter model*®?. 
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GaP 


b) g = 002 
Fig. 20—(110) cross-sectional dark-field TEM images of 
the GaP/Si interface of a Type Ib structure. 


4.3 Defects in GaAsSb crystals grown on 
(001)InP by molecular beam epitaxy 
GaAsSb is an important material for opto- 

electronic and high-speed devices*”). Recently, 

we have reported on the growth of GaAsSb 
crystals on (001)InP substrates by molecular- 
beam epitaxy (MBE)*®). To produce devices 
with suitable properties and reliability, it is 
important to obtain defect-free and thermally 
stable materials. In this section, we describe the 
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TEM evaluation of defects in GaAsSb crystals 
grown on (001)InP substrates by MBE. 

GaAs -x Sb, (x = 0.49) is grown on (001)InP 
substrates by MBE at temperatures in the range 
460-590 °C*®). During growth, the substrates 
are rotated at 1l5r/min to ensure uniform 
composition, thickness and carrier concentra- 
tion. The growth rate and the V/III partial pres- 
sure ratios are kept at 1.2 mm/h and 7.5, respec- 
tively. The epitaxial layers are 2.5 um thick. 

Three types of defect generated at the 
hetero-interface have been observed. They are 
stacking fault tetrahedra, microtwins and misfit 
dislocations. 

4.3.1 Stacking faults 

Figure 21 illustrates a typical optical micro- 
graph of the surface of a GaAsSb crystal grown 
at 590°C. A high density of surface defects is 
observed. They are hillock-like defects lying 
along the <110> direction, and are very similar 
to the oval defects in MBE-Ga(AI)As crystals*??. 
The defect structure observed by TEM at these 
surface defects is shown in Fig. 22. Stacking 
faults lying on the four equivalent {111} planes 
are seen. Since the thin foil includes only the 
surface region of the epitaxial layer, one can 
expect that the surface defects correspond to 
stacking fault tetrahedra®”. The stacking faults 
often multiplied as denoted by X. In addition, 
threading dislocations and small stacking faults 
are found inside the tetrahedra. By considering 
the length of the stacking faults at the surface, 
their basal {111} planes, and the epitaxial layer 
thickness, it is deduced that they originated at 
the hetero-interface. The density of these 
defects can be reduced by lowering the growth 
temperature and cleaning the substrate with 
appropriate As over-pressure. Therefore, the 
dominant factors affecting their generation 
could be the evaporation of phosphorous atoms 
from the substrate and/or the existence of some 
oxide or contamination on the substrate surface. 

4.3.2 Microtwins 

Microtwins are occasionally observed in 
epitaxial layers grown at higher growth tem- 
peratures. They also originate from the hetero- 
interface, extending to the surface area with 
{111} boundaries. Some of them terminate in 
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20 um 


Fig. 21—An optical micrograph of defects revealed on 
the surface of an MBE-grown GaAsSb crystal 
on (001 )InP substrate. 


——— 
0.5 “am 


Fig, 22—A (001) plan-view bright-field TEM image of a 
stacking fault tetrahedron corresponding to the 
surface defect shown in Fig. 21, taken under 
the 220 reflection. 
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Fig. 23—A low magnification (001) plan-view HRTEM 
image of the GaAsSb/(001 )InP interface. 
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mid-layer with {211} boundaries. A HRTEM 
image of a GaAsSb crystal grown at 590°C 
from a region close to the hetero-interface is 
shown in Fig. 23. Brighter regions along the 
interface are probably due to the existence of 
some oxides. In this micrograph, microtwins 
are generated in two regions, I and II. In both 
regions, amorphous-like areas (labeled A), and 
shallow pit-shaped patterns (labeled SP), are 
observed at the interface just below the amor- 
phous-like area. Microtwins (labeled MT) are 
generated only in the region adjacent to the 
areas. 

Based on these results, the following mech- 
anism is suggested for the generation of micro- 
twins: 

1) Evaporation of phosphorous atoms from 
the substrate surface (formation of In-rich 
droplets on the surface); 

2) Generation of amorphous areas as a conse- 
quence of an In-rich melt; 

3) Growth of crystals with different lattice 
constants around the amorphous areas; 

4) Generation of microtwins by the relaxa- 
tion of the stress resulting from the lattice 
mismatch between the normal GaAsSb 
areas and those formed in step 3). 

Thus, microtwins can be eliminated by 
methods similar to those described for stacking 
faults. 

4.3.3 Misfit dislocations 

Since we have successfully grown GaAsSb 
on (001)InP with nearly lattice-matching con- 
ditions (Aa/a <1 x 1073), generation of misfit 
dislocations is unexpected. However, we have 
occasionally found dislocations at the hetero- 
interface with spacing in the range 30-100 nm. 
All of these dislocations have dissociated into 
two partials. Figure 24 is an HRTEM image 
which displays a stacking fault between the 
two partials X and Y. There is clear anisotropy 
in the widths of the stacking fault segments 
contained in the substrate and epitaxial layer. 
This may be due to the different behavior of 
90°- and 30°-partials under different strain 
fields'”). 

A possible model for the formation of 
these dislocations is as follows: 
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Fig. 24—A (110) cross-sectional HRTEM image of a 
dissociated misfit dislocation generated at the 
GaAsSb/(001 )InP interface. 


1) At the initial stage of growth, due to poor 
control of the gas flow (particularly over- 
shooting of the Sb-beam), an Sb-rich ternary 
alloy is grown, resulting in a compressive 
and a tensile strain field in the epitaxial- 
side and the substrate-side, respectively; 

2) As the crystal grows thicker, 60°-type mis- 
fit dislocations are generated by the relaxa- 
tion of the resulting stress; 

3) Due to different behaviors of 90°- and 30°- 
partials, the 90°-partial first extends into 
the substrate, then the 30°-partial extends 
into the epitaxial layer. 

To prevent the introduction of such defects, 
more precise control of the gas flow at the 
initial stage is required. 


4.4 Precipitates in heavily Fe-doped InP 

crystals grown by metalorganic chemical 

vapor deposition 

Semi-insulating InP is an important mate- 
rial for opto-electronic and high-speed devices. 
Until now, much effort has been devoted to 
achieving high resistivity®!)5>), It has previous- 
ly been reported that resistivities of 10° Q:cm 
have been successfully obtained for Fe-doped 
InP grown by MOVPE*”?. However, iron- 
phosphorous precipitates are produced in the 
crystal during growth. Nakahara et al.°5) have 
undertaken a detailed investigation of the pre- 
cipitates by TED analysis and EDX and have 
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concluded that they are composed of FeP. 
However, little is known about the nature and 
distribution of the precipitates, or their genera- 
tion mechanism. In this section, more detailed 
characteristics of the precipitates and the rela- 
tionship between the Fe-doping gas flow rate 
and precipitate density are described based on 
the TEM analysis. 

InP crystals are grown on (001)InP sub- 
strates at 650°C by atmospheric pressure 
MOVPE. Phosphine (PH,) and trimethylindium 
(TMIn) are used as reactants. The V/III partial 
pressure ratio is kept at 20. The electron con- 
centration of the non-doped layer is about 
2x 10'* cm”. Ferrocine Fe(C;H;), vapor is 
used as the Fe-doping gas and its flow rate is 
varied in the range 0-100 ml/min. 

4.4.1 Microstructure of precipitates 

First of all, an Fe-doped InP layer grown 
with the maximum gas flow rate (100 ml/min) 
was examined. As_ previously reported by 
Nakahara et al.°*), a large number of circular 
precipitates are observed in a dark-field image, 
particularly in the area close to the edge of the 
thin specimen. They may be expected to stay 
on the surface of the specimen after etching 
since they are thought not to be dissolved by 
the Br,/CH,OH etchant®*’. However, in the 
region slightly away from the edge of the 
specimen, the density of precipitates is found 
to be nearly constant as shown in Fig. 25. 
Figures 25a) and b) are bright-field and weak- 
beam dark-field TEM images of precipitates 
taken under the 220 and the 660 reflections 
with positive s, which is a deviation parameter 
from an exact Bragg condition. 

From these images, one can differentiate 
between precipitates buried in the matrix and 
those lying on the surface as follows: 

1) Precipitates exhibiting parallel Moiré fringes 
with constant spacing are buried in the 
matrix (we call these Type A precipitates); 

2) Precipitates which do not exhibit Moiré 
fringes are lying on the surface (Type B). 

In the latter situation, the probability of 
direct double reflection between precipitates 
and matrix is considered to be much lower as 
a result of possible oxide presence and the small 
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a) Bright-field image 


b) Weak-beam dark-field image 


Fig. 25—Plan-view TEM images of Fe-precipitates in sample | (the gas flow rate is 100 ml/min) after chemical etching. 


contact area between the two materials. Type B 
precipitates are generally out of contrast in the 
weak-beam image since there is no strain associ- 
ated with them. However, they occasionally do 
exhibit bright contrast presumably as a result 
of sticking force between the precipitate and 


110 


110 


the matrix at the surface. In order to confirm 
this, bright-field images were obtained under 
many beam conditions. 

Typical images of Type A precipitates ob- 
tained from the same specimen are illustrated 
in Fig. 26. Note that all of these precipitates 


10 nm 


Fig. 26—A series of plan-view bright-field TEM images of precipitates taken under many beam diffraction 


conditions [001] zone-axis condition. 
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Note the Moiré fringes in three different orientations. 
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5nm 
Fig. 27— A plan-view HRTEM image of the precipitate 
taken under nearly zone-axis condition. Two 
sets of 200-lattice frings can be seen. A small 
amount of lattice distortion is visible at the 
precipitate-matrix interface (see arrows). 


exhibit three kinds of Moiré fringe with differ- 
ent orientations and periodicities. Therefore, 
there is a certain orientation relationship be- 
tween the precipitates and the matrix crystal. 
Although nano-probe electron diffraction ex- 
periments (2-5 nm probe dia) were performed 
on the precipitates, the exact orientation of the 
precipitates could not be determined due to 
double reflection of electrons at the precipi- 
tate-matrix interface. Figure 27 is a plan-view 
HRTEM image of a precipitate. The precipitate 
also exhibits Moiré fringes similar to those 
shown in Fig. 26. Lattice fringes corresponding 
to dx99 for InP can be seen. Slight misalignments 
of lattice fringes are occasionally found at the 
precipitate-matrix interface (see arrows). This 
result suggests that there is little strain field at 
the interface, which is consistent with the fact 
that the lattice parameter of InP is very close to 
that of FeP. Thus, the precipitates are coherent. 
Although precipitates were determined to be 
FeP by TED analysis®**), compositional analysis 
using EDX was carried out on them. Figures 
28a) and b) are typical EDX spectra obtained 
from the matrix and the precipitate using an 
electron probe 2 nm in dia. Fe is clearly detect- 
ed within the precipitate | see Fig. 28b)}. 

4.4.2 Spatial distribution of precipitates 

A (110) cross-sectional TEM observation 
was performed on Fe-doped InP crystals grown 
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b) Precipitate area 
Fig. 28—EDX spectra from the Fe-doped InP crystal 
(probe dia is nearly 2 nm). 


with the maximum gas flow rate; the resulting 
image is shown in Fig. 29. High magnification 
images of the three areas indicated by a, b and 
c are shown in inserts a’, b’, c’, respectively. 
Precipitates are again observed as circular images 
indicating that they are spherical particles. 
Also note that the precipitates are uniformly 
distributed in the crystal. Since the interface 
between the epitaxial layer and the substrate 
is not clearly observed |the broken line in 
Fig. 29 is assumed to be close to the real inter- 
face from consideration of the thickness of the 
epitaxial layer determined from a scanning elec- 
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tron micrograph (SEM) image of a cleaved and 
stain-etched cross-section of the homo-junc- 
tion}, it seems that there is no preferential 
segregation of Fe-atoms at the interface. More- 
over, precipitates attached to the surface were 
not found at all. 

4.4.3 Dependence of the density and diam- 

eter of precipitates on gas flow rate 

The effect of the Fe-doping gas flow rate 
on the density and diameter of precipitates was 
also studied. Figures 30a) to d) show plan-view 
bright-field TEM images of precipitates in 
Fe-doped InP crystals grown with gas flow rates 
of 100, 30, 20 and 5 ml/min, respectively. 
Precipitates are observed in crystals grown with 
gas flow rates higher than 20 ml/min. Their 
density increases as the gas flow rate increases. 
However, their diameter is almost constant. 
Precipitates are not observed in crystals grown 


Surface 


with gas flow rates lower than 5 ml/min. These 

results are summarized in Table 1. From second- 

ary ion mass spectroscopy (SIMS) measure- 

ments, the total Fe concentration in sample 4 

(gas flow rate 5 ml/min) in Table | is calculated 

to be approximately 1 x 10'7 cm™?, which is 

the solubility limit of Fe in InP at 650 °C. 
4.4.4 Generation mechanism of precipitates 
Based on the results described in subsections 

4.4.1 to 4.4.3, three generation mechanisms for 

FeP precipitates can be proposed as follows: 

1) FeP particles are first formed in the gas 
phase, then they are deposited onto the 
growth surface and are buried in the growing 
InP crystal. 

2) FeP precipitates are formed directly on the 
growth surface through migration and segre- 
gation of Fe-atoms, and the strong reaction 


between Fe-atoms and P-atoms. 
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/ 
/ 


Y 


InP epitaxial layer 
eA 


Uf 
7 


7 InP substrate 


es 
0.2 wm 


Fig. 29— A (110) cross-sectional bright-field TEM image from an Fe-doped InP crystal taken under g=220,s>0. 
Inserts labeled a’, b’ and c’ correspond to high magnification images of the areas denoted by a, b, and c. 
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3) In the growing InP layer, when the doped critical value, they can not grow any more 
Fe-concentration exceeds the solubility due to the strain field accumulated around 
limit, supersaturated Fe-atoms tend to segre- them. 
gate at some nucleation centers via rapid Models 1) and 2) are unlikely since no FeP 
diffusion, and FeP precipitates start to grow. particles are found on the surface of As-grown 
When the diameter of precipitates reaches a InP layers, and no precipitates are observed in 


Table 1. Dependence of Fe-doping gas flow rate on the density and diameter of precipitates 


Sample Gas-flow Epitaxial layer Precipitate Precipitate diameter (nm) 
No. (ml/min) oy m3 Average (max/min) 
4] 100 2~=~———<“<«‘ x10 8.3 (3.8/11.2) 
2 30 2 3x 10° 5.8 (4.1/11.6) 
3 20 2.5 8 x 10)? 12.9 (5.8/18.6) 

4 5 4 <5 x 10/9 = 

5 l 3 <5x10!° = 

6 ; 3 <5«10'° = 

7 0.2 3 <5 X10" = 


Oi (a. 


b) 30 ml/min 


4 = 


« . 
8220 
iia F: 
as 
dl 
————! 
0.1 ~m 

c) 20 ml/min d) 5 ml/min 


Fig. 30—Typical plan-view bright-field TEM images of FeP-precipitates in Fe-doped InP crystals grown 
under different gas flow rates, obtained by g = 220,s>0. 
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crystals grown with gas flow rates lower than 

5 ml/min. 

At present, model 3) is the most probable, 
based on the following evidence: 

1) Precipitates are generated when the doped 
Fe-concentration exceeds the solubility 
limit; 

2) Precipitates are uniformly distributed in the 
crystal and have almost constant diameter; 
and 

3) Precipitates are spherical and coherent. 


5. Evaluation of ordered structures in III-V 
alloy semiconductors 

Apart from the materials issue of defect 
generation, in most III-V alloy semiconductors 
there are two major issues regarding the thermal 
instability of the materials; namely the genera- 
tion of “modulated structures®®)®°” and the 
generation of “ordered structures®!)“®>)”’, 

The modulated structures are generated in 
alloy semiconductors such as InGaAs, InGaP and 
GaAsSb whose compositions are inside the 
spinodal region®®). They are observed as quasi- 
periodic diffraction contrast in bright- or dark- 
field images under two-beam or multi-beam con- 
ditions. They develop in the two equivalent 
directions of <100> and <010>, and they are 
columnar in the <00Ol> growth direction. 
It is also clear that structures with both long 
(50-200 nm) and short (5-20 nm) periodicities 
are formed when the composition is outside the 
spinodal region®’)©®), Thus, one can expect 
that the former structures are formed during 
growth by spinodal decomposition, and the later 
ones during the cooling process after growth. 
Moreover, the compositional fluctuations along 
the structure measured by EDX are in the 
range 2-3 percent®®’, which is far smaller than 
expected. This may be due to the fact that 
in semiconductors, atomic diffusion in the bulk 
or on the surface is very slow compared with 
metal alloys. Their generation does not depend 
on the growth method. The decomposition 
process may be enhanced by atomic diffusion 
at the liquid-solid interface for LPE growth and 
at the growth surface for MOVPE, VPE and 
MBE. 
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CuAu-l CuPt 
a) Unit cells 


Ll Bp 


CuPt Chalcopyrite 
b) Reciprocal lattices 


Chalcopyrite 


Fig. 31—Schematic diagrams of three types of 
substitutionally ordered structure in the fcc 
lattice of AB-type binary compounds. 


On the other hand, recently, ordered struc- 
tures (or natural superlattices), in which atoms 
of two different elements in column III or V 
arrange periodically in lattices, have been 
observed in various III-V alloy semiconduc- 
tors®!)-85) In most cases, the ordered structures 
are generated locally in the crystal, and the 
degree of ordering is strongly affected by the 
growth method and growth conditions. It is 
known that there are three types of substitu- 
tionally ordered structures in face centered 
cubic (fcc) lattices of AB-type binary alloy sys- 
tems. They are CuAu-I (LI,), CuPt (LIl,) and 
Chalcopyrite (El,), whose atomic arrangements 
and corresponding reciprocal lattices are sche- 
matically shown in Fig. 31°4). 

The results for ordered structures in II]-V 
alloy semiconductors obtained up to now are 
summarized as follows: 

1) The most. common structure in crystals 
grown on.a (001) substrate is CuPt-type (see 
Fig. 31), and other structures are only occa- 
sionally observed. 

2) They are very often generated in crystals 
grown by MOVPE, VPE and MBE, whereas 
no ordering takes place in LPE-grown crys- 
tals |except for one report of ordered struc- 
tures in InGaAs by Nakayama et al.””, 
which is not yet confirmed elsewhere} . 
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3) Regarding CuPt-type structures, one can ob- 
serve only two variants in the (110) cross- 
section. In particular, atomic steps on the 
growth surface are believed to play an im- 
portant role in the generation of the ordered 
structures, since one of the two variants is 
preferentially enhanced when substrates 
tilted toward the <110> direction are used. 

4) The ordering is not perfect and the ordered 
regions are plate-like microdomains lying on 
planes nearly parallel to the growth surface. 

5) Defects such as anti-phase boundaries are 
often generated in the ordered region. 

6) The degree of ordering depends on the 
growth temperature, V/II] partial pressure 
ratio (in the case of MOVPE and MBB), and 
rotation velocity of the substrate. 

7) Only strongly ordered InGaP crystals grown 
by MOVPE exhibit abnormal band gap 
energies (up to 50 meV lower than the nor- 
mal value). 

In this chapter, we describe the results of 
detailed TEM evaluations of ordered InGaP crys- 
tals grown on (001)GaAs and ordered InGaAs 
crystals grown on (110)InP. 


5.1 Ordered structures in InGaP crystals grown 
on (001)GaAs by MOVPE 
Regarding atomic ordering, Gomyo et a 
have been the first to report that in InGaP crys- 
tals grown by MOVPE, there is a strong relation- 
ship between the optical properties of the crys- 


1,6) 


tals and the degree of ordering. They have also 
concluded from the electron diffraction analysis 
that they are most likely to be superlattices of 
(110) planes ... Ga/Ga/In/In/Ga/Ga/In/In/ ... 
in the direction of [110] in the column II sub- 
lattice, i.e. (InP),(GaP), super-lattices on the 
(110) plane. However, evidence of ordering is 
available only from a (001) plan-view specimen, 
and the three-dimensional atomic structure of 
the ordered InGaP was not clear. 

In this section, the three-dimensional atomic 
structure of the MOVPE-grown InGaP crystals 
on (001)GaAs is described in more detail based 
on plan-view and cross-section TED analyses, 
and HRTEM. 

Epitaxial growth is carried out by atmos- 
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Fig. 32—A TED pattern from the (001) plan-view of a 
Type I InGaP crystal. 


pheric pressure MOVPE. The InGaP layers are 
grown on (Q01)-oriented GaAs substrates tilted 
2.5° towards <110>, and the growth tempera- 
ture is kept at 630 °C. PH;, TMIn and TMGa are 
used as reactants, and the V/III partial pressure 
ratio is kept at 100. The epitaxial wafers exam- 
ined in this experiment consist of a GaAs buffer 
layer and an InGaP layer of 800-900 nm and 
10-400 nm, respectively. Two types of InGaP 
crystals are examined. One is an InGaP crystal 
grown with rotation of the substrate (20 r/min, 
Type I), and the other is an InGaP crystal grown 
without rotation of the substrate (Type II). In 
both types of InGaP crystal, the photolumines- 
cence peak wavelength of the crystals is 652 nm 
at 77K and the typical lattice-mismatch Aa/a 
between GaAs and InGaP is 2 x 10°3. 

5.1.1 Evaluation of Type I crystals 
1) Plan-view observation 

First of all, plan-view observations were 
carried out. In the dark-field TEM image of the 
InGaP crystal taken under the 220 reflection, 
both coarse and fine modulated structures with 
periodicities of 30-50 nm and 5-10 nm, respec- 
tively, are clearly observed along the two equiv- 
alent directions of [100] and [010]*°”®”. 
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a) (110) cross-section 


b) (110) cross-section 


Fig. 33—TED patterns from the (110) and the (110) cross-section of a Type I InGaP crystal. 


Next, a TED analysis was performed on 
the (001) plan-view specimen as illustrated in 
Fig. 32. Extra-spot pairs are observed on both 
sides of the matrix spots along only one <110> 
direction. 

The extra-spot pairs are indexed as nearly 
[h +1/2, k ¥1/2, 0] for an [hkO] matrix spot®!. 
However, careful measurement of the positions 
of these extra spots confirmed that each one is 
located slightly inside the ideal position de- 
scribed above. 

2) Cross-sectional observation 

Cross-sectional TEM observation has been 
carried out on both the (110) and the (110) 
orientations. Pronounced extra spots (here refer- 
red to as “superstructure spots’’) and streaks are 
present in the (110)TED pattern {see Fig. 33a)}, 
but not in the (110)TED pattern {see Fig. 33b)}. 

The characteristic features of the diffrac- 
tion pattern are as follows: 

i) the presence of superstructure spot 

pairs at positions indexed as 
{h 1/2, k ¥1/2, 1 +1/2] for an [hkl] 
matrix spot {for example, the positions 
indicated by arrows in Fig. 33a)}; 

ii) the appearance of streaks around the 

superstructure spots lying along both the 
[001] and the [001] directions, but 
slightly tilting (a few degrees) toward 
the [110] and the [110] directions, 
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Fig. 34—Schematic diagram of the reciprocal lattice of 
an ordered InGaP crystal (CuPt-type). 
Large and small closed circles represent the 
positions of the fundamental diffraction spots 
and superstructure spots. 


respectively {however, no superstructure 
spots are found at positions indexed as 
pnt /2,. k 21/2 
Fig. 33b), nor are any observed at posi- 
tions indexed as [h 1/2, k ¥1/2, O] as 
shown in Fig. 33a) which are expected 
from the (001)TED pattern shown in 
Fig. 32}; 

ili) in both types of diffraction pattern, 
weak and short streaks are observed 
along the [001] and the [001] direc- 
tions around each matrix spot. 


1 +1/2], as shown in 
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On the basis of the three diffraction patterns 
described above, a reciprocal lattice of the 
InGaP was constructed schematically as shown 
in Fig. 34. This leads us to the conclusion that 
the extra-spot pairs appearing in the (001)TED 
pattern (see Fig. 32) correspond to intersecting 
points of the streaks on the (001) reciprocal 
plane and that they are not real superstructure 
spots®!), The fact that the extra spots appearing 
in the (001)TED pattern are located slightly in- 
side the ideal positions of the [h +1/2, k ¥1/2, 0] 
is well understood by the tilting of the streaks. 
The positions denoted by Pp and Qp are indexed 
as [1/2, 1/2, 0] and [1/2, 1/2, 0], respectively, 
and P and Q show observed locations of pair 
spots on both sides of the 000 spot in Fig. 32. 
Since streaks are present around the superstruc- 
ture spots, it is assumed that the ordering is not 
perfect, and that ordered regions correspond to 
microdomains of a finite size. The reciprocal 
lattice illustrated in Fig. 34 is associated with a 
CuPt-type ordered structure®”?, which is one of 
these previously shown in Fig. 31. The ordering 
occurs on the (111) plane of column II atoms 
with doubling in periodicity, ice. the ordered 
structure corresponds to a monolayer superlat- 
tice of (InP), (GaP), on the (111) plane. 

A (110) cross-sectional HRTEM image of the 
InGaP crystal is illustrated in Fig. 35. Doubling 
in periodicity of the (111) lattice fringe, i.e. the 
presence of arrays of brighter spots on each 
second (111) fringe is clearly observed. This is 
presumably associated with the superstructure 
spots in the (110)TED pattern {see Fig. 33a)}, 
whereas no abnormal contrast is observed in the 
(110) cross-sectional HRTEM image. The 
presence of non-ordered regions is very often 
observed, indicating that the ordering is not 
perfect. These ordered structures are observed in 
all InGaP crystals with different thickness. 

As described above, it should be recognized 
that the ordering occurs on only one of the four 
equivalent (111) planes. In order to clarify the 
cross-sectional shape of the ordered region, a 
dark-field image is taken directly from one of 
the superstructure spots, as shown in Fig. 36. 
In this image, brighter regions correspond to 
ordered regions since the intensity is associated 
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Fig. 35— A (110) cross-sectional HRTEM image of 
ordered regions in a Type I InGaP crystal. 
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Fig. 36—A (110) cross-sectional dark-field TEM image of 
ordered regions in a Type I InGaP crystal, 
obtained from one of the superstructure spots. 


with only the superstructure spot. The ordered 
regions are found to be plate-like microdomains 
of 5-15 nm wide and about 2 nm thick, lying on 
the plane slightly inclined to the (001) plane. 
These results are consistent with the appearance 
of streaks around the superstructure spots 
and their tilting in the TED pattern shown in 
Fig. 33a). 

5.1.2 Evaluation of Type Il crystals 
1) Plan-view observation 

Plan-view observations reveal both coarse 
and fine modulated structures with periodicities 
similar to those observed in Type crystals. 


Si 
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a) Region-1 


5 


b) Region 2 


Fig. 37—TED patterns from a (110) cross-section of a 
Type II InGaP crystal. 


Thus, there is no influence of substrate rotation 
on the generation of modulated structures. The 
electron diffraction patterns obtained from 
Type II crystals also exhibit superstructure spots 
similar to those observed in the Type I crystals. 
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2) Cross-sectional observation 

Quite different intensity profiles of the 
superstructure spots are found in the TED pat- 
terns obtained from two different regions of a 
Type II crystal {see Figs. 37a) and b)}. Although 
one series of superstructure spots is present, 
they are ‘“‘spotty”’ rather than “streaky’’, indicat- 
ing fairly long range ordering. Moreover, the 
intensity of the superstructure spots is much 
stronger than in Typel crystals, and weaker 
streaks can also be seen along the [001] direc- 
tion. Occasionally, two series of superstructure 
spots are observed as illustrated in Fig. 37b) 
(see spots indicated by a, and a,). However, 
there is asymmetry in the intensity of each 
series. Furthermore, very weak superstructure 
spots, which are associated with a CuAu-I-type 


80),81 
structure®”? 2 


are occasionally observed at 
positions indexed as 001, 003, and 110 {see also 
spots indicated by b in Fig. 37b)}. These re- 
sults suggest that CuPt-type structures are pre- 
dominantly generated in the crystals grown on 
(O01 )-oriented substrates. 

Figure 38a) illustrates an HRTEM image ob- 
tained from a (110) cross-section of Type IJ 
crystal. Doubling in (111) lattice fringes is clear- 
ly observed over a considerable region; this indi- 
cates long range ordering. It is also emphasized 
that planar defects are often observed in the 
high resolution image, i.e. the region indicated 
by APBs in Fig.38b), where a phase-shift of 
the (Clits fringes is observed. They might be 
anti-phase boundaries as suggested by Kuan et 
al®®. The HRTEM image from a (110) cross- 
section does not exhibit periodic intensity in 
any lattice fringes. 

A cross-section dark-field image of ordered 
regions obtained from one of the superstructure 
spots in a Type Il crystal is shown in Fig. 39. 
The size and shape of the ordered domains are 
found to be very different. In most cases, the 
domains are much larger than those observed in 
Type I crystals (see Fig. 36). These domains also 
consist of plate-like segments lying nearly on the 
(001) plane, possibly associated with anti-phase 
boundaries. In.this case, approximately 80 per- 
cent of the total area is dominated by ordered 
regions, which is much more than in Type I 
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a) An HRTEM image 5nm 


b) A higher magnification image of a part of the 
ordered region with anti-phase boundaries 
shown in a) 


Fig. 38—(110) cross-sectional HRTEM images of ordered 
structures in a Type II InGaP crystal. 


crystals. Moreover, the boundaries between 
ordered domains and non-ordered regions are 
also slightly tilted. 

5.1.3 On the difference in the generation of 
ordered structures between Type I 
and II crystals 

In order to clarify the difference in the gen- 

eration of ordered structures, one must consider 
their generation mechanism. 

Based on the findings of the ordered struc- 
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Fig. 39—A (110) cross-sectional dark-field TEM image of 
ordered regions in a Type II InGaP crystal, 
obtained from one of the superstructure spots. 


tures in MOVPE-grown InGaP crystal described 
in subsections 5.1.1 and 5.1.2, and those in 
other III-V alloy semiconductors®!)§?), it is 
realized that ordering occurs in III-V alloy semi- 
conductors grown by MOVPE or MBE where 
deposited atoms can migrate and rearrange easily 
on the surface during growth. Thus, the mobility 
of the matrix atoms on the surface is expected 
to be one of the most dominant parameters for 
the generation of ordered structures. One of the 
most direct pieces of evidence for this hypothe- 
sis can be found in the cross-sectional feature of 
ordered domains shown in Fig. 38, i.e. the bound- 
aries between the ordered and non-ordered 
regions are tilted from the [001] direction, 
indicating that the boundaries are moving to- 
wards the [110] direction during crystal growth. 
In view of these considerations, one can 
suggest the following as factors which might 
affect the formation of ordered structures 
during crystal growth with or without substrate 
rotation: 
1) Mobility of atoms on the growth surface; 
2) Distribution of the gas flow near the growth 
surface; 
3) Effective temperature at the growth surface. 
Since there is no major difference between 
the average thickness in Type I crystals and that 
in Type II crystals, factor 3) seems to have little 
contribution to the ordering. 
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Thus, one can assume that a combination of 
factors 1) and 2) affect the formation of ordered 
structures on the growth surface as follows: 

i) When the InGaP crystal is grown on a 
rotating substrate, local fluctuation is 
present in the gas flow near the growth 
surface, reducing the surface rearrange- 
ment of the deposited atoms. In this 
case, migration of the deposited atoms 
on the growth surface is much more dis- 
turbed by the gas molecules with rela- 
tively high velocities than is the case 
with a non-rotating substrate. Thus, the 
effective mobility of the deposited 
atoms is lower. 

ii) When the crystal is grown on a non- 
rotating substrate, gas flow distribution 
on the growth surface is uniform and 
effective mobility of the deposited 
atoms is high compared with case i); 
thus ordering is enhanced. 


5.2 Ordered structures in InGaAs crystals grown 
on (110)InP by MBE 
As described in section 5.1, the ordered 
structures in crystals grown on (001) substrates 
are mostly CuPt-type. On the other hand, CuAu- 
I-type ordered structures have been found by 


Li 


a) (110) plan-view 


Kuan et al. in AlGaAs crystals grown on 
(110)GaAs substrates by MOVPE and MBE®”, 
and InGaAs crystals grown on (110)InP sub- 
strates by MBE®!). However, they only showed 
TED patterns corresponding to the CuAu-l- 
type structure; the microstructural features of 
the ordered regions and dependence of the 
degree of ordering on growth conditions have 
not yet been clarified. In this work, we describe 
a detailed TEM study on the dependence of the 
degree of ordering on the growth conditions and 
a microstructural evaluation of the CuAu-I-type 
ordered structures. 

Undoped InGaAs crystals grown on 
(110)InP substrates by MBE were examined 
in this study. The crystals were grown at 
360-485 °C on exactly (110)InP substrates, or 
on InP substrates tilted toward the [001], 
[001] or [110] direction. During growth, the 
substrates were rotated to achieve a good uni- 
formity in composition, thickness and carrier 
concentration. The V/IIl partial pressure ratio 
was kept at 400 or 60 {only for crystals grown 
on substrates tilted toward the [001] direc- 
tion}. The epitaxial layer was approximately 
1 um thick. The lattice mismatch between the 
InGaAs layer and the InP substrate was within 
tl x 10", 


b) (110) cross-section 


Fig. 40—TED patterns from the (110) plan-view and the (110) cross-section of InGaAs grown on (110)InP 


tilted 3° toward [110]. 
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5.2.1 Identification of CuAu-I-type struc- 
tures by TED and dependence of the 
degree of ordering on substrate 
orientation 

Figures 40a) and b) show the TED patterns 

from a (110) plan-view and a (110) cross-section 
of an InGaAs crystal grown on a (110)InP sub- 
strate tilted 3° toward the [110] direction. In 
both figures, superstructure spots are observed at 
positions indexed as 001, 110,110, 112, 112,... 
This particular set of superstructure spots is 
associated with a CuAu-I-type ordered struc- 
ture®®)§") In this structure, Ga atoms preferen- 
tially occupy the (000) and (1/2, 1/2, 0) sites 
and In atoms preferentially occupy the 
(1/2, 0, 1/2) and (0, 1/2, 1/2) sites in each unit 
cell. A schematic diagram of the atomic arrange- 
ment in the CuAu-I-type structure viewed along 
the [110] direction normal to the growth axis is 
shown in Fig.41. From this diagram, one can 
see that the perfectly ordered material consists 
of alternating InAs and GaAs monolayers, i.e. 
(InAs),;/(GaAs),; monolayer superlattices, when 
viewed along either the [110] growth direction 
or the [001] direction. Note that the super- 
structure spots shown in Fig. 40b) are extremely 
streaky, and that the streaks are ‘“S’’-shaped, 
tilted slightly from the [110] direction. Super- 
structure spots of a similar shape are also 
observed in the TED patterns from (110) cross- 
sections of ordered InGaP (CuPt-type) crystals 
grown on (001)GaAs tilted substrates. 

The streaky shape of the superstructure 


a) Exactly (110) 


b) (110)InP tilted 3° toward [001 ] 
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Fig. 41—A schematic diagram of the atomic arrangement 
of CuAu-I-type ordered structure in an InGaAs 
crystal viewed from the [110] direction. 


spots may be due to either 

1) the generation of anti-phase boundaries in 
the ordered region or 

2) the fact that ordered regions are plate-like 
microdomains (see Fig.46 in subsection 

5.2.3), although computer simulations of 

epitaxial growth are required to resolve this 

issue. 

Figures 42a) to c) show TED patterns from 
plan-views of InGaAs crystals grown on 
(110)InP substrates with tilting angles of 0°, 
3° and 5° toward the [001] direction. In 
crystals grown on a (110) pure InP substrate 
\see Fig. 42a)}, the intensity of the superstruc- 
ture spots is very weak, indicating very weak 
ordering. However, when a substrate tilted 3° 
toward the [001] direction is used, ordering 


c) (110)InP tilted 5° toward [001] 


Fig. 42—TED pattern from the (110) plan-view of InGaAs crystals grown at different tilting angles. 
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becomes stronger as shown in Fig. 42b). Fur- 
ther, ordering becomes much stronger in a crys- 
tal grown on a substrate tilted 5° toward the 
[001] direction {see Fig. 42c)}. Since tilting of 
the substrate orientation toward the [001] or 
the [001] direction introduces periodic arrays 
of steps on the substrate surface, one can con- 
clude that atomic steps on the growth surface 
play an important role in the formation of the 
ordered structures. It is also established that the 
ordering is stronger in crystals grown on sub- 
strates tilted toward the [001] or the [001] 
direction than those on substrates tilted toward 
the [110] direction. This may be explained as 
follows: when the (110) surface is tilted toward 
the [001], the steps are very straight, giving rise 
to strong ordering, but when the (110) surface is 
tilted toward [110], highly kinked steps are 
formed, preventing the crystal from ordering 
(in the latter case, only the kinks themselves can 
enhance the ordering). 

5.2.2 Dependence of the degree of ordering 

on growth temperature 

In order to clarify the influence of growth 
temperature on the formation of ordered 
structures, we grew InGaAs crystals at various 
temperatures in the range of 435-485 °C on 
(110)InP substrates tilted 3° toward the [110] 
direction. From TED analysis, it was found that 
in this temperature range, the degree of ordering 
increases with growth temperature. Figures 43a) 
and b) show TED patterns from (110) cross- 
sections of InGaAs crystals grown at 435 °C and 
485 °C, respectively. TED analysis has also been 
carried out for crystals grown on substrates 
tilted toward the [001] direction. It has been 
found that, in the range of 360-450 °C also, the 
degree of ordering increases with growth tem- 
perature. These results can be explained by the 
fact that the mobility of atoms deposited on the 
growth surface also increases with temperature, 
and that the formation of ordered structures 
is thought to be strongly related to the migra- 
tion and reconstruction of the deposited 
atoms®?)®?) , 
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b) Tg = 485 °C 
Fig. 43—TED patterns from the (110) cross-sections of 
InGaAs crystals grown on (110)InP tilted ai 
toward [110]. 


5.2.3. Microstructural characterization of 
CuAu-l-type structures 

HRTEM analysis was carried out to eval- 
uate the ordered structures at an atomic level. 
Figure 44 shows a typical HRTEM image of 
ordered InGaAs grown on a (110)InP substrate 
tilted 3° toward the [001] direction. Doubling 
of the (002) (see arrows X,-X3) and the (220) 
(see arrows Y,-Y,) lattice fringes, which is 
associated with CuAu-I-type structures, is ob- 
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served locally, suggesting that the crystal con- 
sists of both well-ordered and non-ordered 
regions. It should be emphasized that planar 
defects are often observed in high resolution 
images, e.g. the region denoted by APB, where 
a phase-shift of the (002) lattice fringes is 
observed. These are expected to be anti-phase 


2nm 
Fig. 44—A (110) cross-sectional HRTEM image of 
ordered InGaAs grown on (110)InP tilted 3° 


toward [001]. 
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a) A schematic diagram of the ordered structure 
(large and small dotted circles represent 
foreground and background atoms) 


b) Af (defocus) = 38 nm, 
t (thickness) = 8.3 nm 
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boundaries as suggested by Kuan et al®!. In the 


(001) cross-section HRTEM images, doubling of 
both the (220) and the (220) lattice fringes is 
observed. Figure 45 shows computer simulated 
images of a CuAu-I-type ordered structure in 
InGaAs. Figure 45a) shows a schematic diagram 
of the ordered structure viewed along the [001 ] 
direction. Figures 45b) to d) show the thickness 
dependence of a CuAu-I-type structure with an 
incident electron beam of [110] and optimum 
focusing. It is found that in each case, the In-As 
pairs cause a brighter spot than the Ga-As pairs. 

In order to clarify the cross-sectional shape 
of the ordered region, a dark-field image was 
taken with one of the superstructure spots, as 
shown in Fig. 46. In this image, ordered regions 
are observed to be brighter since the image con- 


—————— 
0.1 wm 


Fig. 46 —A (110) cross-sectional dark-field TEM image of 
an ordered InGaAs crystal from one of the 
superstructure spots. 


c) Af = —38 nm, ¢t = 16.3 nm 


d) Af = —38 nm, t = 25 nm 


Fig. 45—Computer simulated images of ordered InGaAs by multi-slice method. 
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trast is formed by only the superstructure spot. 
It is found that the ordered regions are plate-like 
microdomains 2-4 nm thick, lying on a plane 
slightly tilted from the (110) plane. It is thought 
that the streaks observed in the TED pattern |see 
Fig. 40b)} are due to the plate-like shape of 
these microdomains. The tilting of the streaks 
can also be explained by the tilting of the micro- 
domains. 

On the basis of these results, it is proposed 
that the formation of CuAu-I-type structure 
takes place by 
1) migration of the deposited atoms on the 

growth surface, and 
2) rearrangement of the migrating atoms from 

the step edges over two atomic layers. 

Further experiments and theoretical calcula- 
tions are required to completely understand the 
formation mechanism of the ordered structures. 

5.2.4 Enhancement of two-dimensional 

electron gas mobility in InGaAs/ 
N-InAlAs heterostructures 

In order to investigate the electrical proper- 
ties of ordered (110)InGaAs, we studied two- 
dimensional electron gas (2DEG) mobility in 
InGaAs/N-InAlAs heterostructures with a 
strongly ordered InGaAs layer®5)86) | We grew 
selectively doped InGaAs/N-InAIAs heterostruc- 
tures on (110)InP substrates tilted 3° and 5° to- 
ward the [001] direction. These heterostruc- 
tures consisted of a 600 nm undoped InAlAs 
buffer layer, a 600 nm undoped InGaAs channel 
layer, a 10 nm undoped InAlAs spacer layer, a 
90 nm Si-doped InGaAs cap layer with a carrier 
concentration of 3 x 10'? cm™3, and a 10 nm 
Si-doped InGaAs cap layer with a carrier con- 
centration of 3 x 10'7 cm™3. The 2DEG mobili- 
ty and sheet electron concentration (N,) are 
plotted as a function of temperature in Fig. 47. 
The 2DEG mobility is highly anisotropic: in the 
<001> direction, the mobility (denoted by 
closed squares) saturates at 100000 cm?/V:s 
and is comparable to that in the sample grown 
on a (001) substrate (see line); while in the 


<110> direction, the mobility (denoted by 
closed triangles and closed circles) is much 
higher than the alloy scattering limited mobility 


44 


InGaAs/N-InAlAs on (110) InP 


Ashaa, 
ee e0e e, 


eo 
e 
A 


Hao 3° tilted 
oe 
Y. Takeda 


5° tilted 


Ns ( 10!? cm7*) 


Lor —ty 
rp ¢ 
A 
Hoo) 


\ 


Mobility (x 10° cm?/V +s) 


Ns RRR AO 


0 
1 10 100 1 000 
Temperature (K) 


Fig. 47— Temperature dependence of the 2DEG mo! ty 
and the sheet electron concentration in the 
<110 >and < 001 > directions in selectively 
doped InGaAs/N-InAlAs hetero-structures 
grown on (110)InP substrate tilted 3° and 
toward < 001 >. 


calculated at an N, of 1 x 10!? cm™? by Takeda 
et al.®”) (denoted by an arrow) at lower tem- 
peratures. The 2DEG mobility in the <110> 
direction in a sample grown on 3° tilted sub- 
strate reached 161000 cm?/V-s with an N, of 
9.5 x 101! cm? at 6 K (101 000 cm?/V's with 
an N, of 9.6 x 10'! cm? at 77 K)®®. This is the 
highest mobility, to our knowledge, ever 
reported for lattice-matched InGaAs/N-InAlAs 
heterostructure systems. This mobility enhance- 
ment in the <110> direction is thought to be 
caused by the suppression of alloy scattering due 
to the formation of CuAu-I ordering in the 
InGaAs channel layer. On the other hand, the 
2DEG mobility in the <001> direction com- 
pares with that of the non-ordered sample. The 
difference in 2DEG mobility in different direc- 
tions could be explained as follows: Since the 
plate-like ordered and disordered regions are 
stacked on nearly the (110) plane on the tilted 
substrates, stripes of ordered and disordered 
regions are formed along the <110> direction 
at the InGaAs side of the hetero-interface where 
the 2DEG accumulates. Thus, along the <001> 
direction, i.e. normal to the stripes, electrons 
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travel across the disordered regions. Along the 
<110> direction, i.e. parallel to the stripes, the 
2DEG formed in both the ordered and disor- 
dered regions move along the stripes. The elec- 
trons in the ordered regions should be unaf- 
fected by alloy scattering in the disordered 
regions. 

The highest 2DEG mobility obtained in the 
<110> direction, however, is lower than the 
theoretical prediction in alloy scatter-free 
(InAs), (GaAs), monolayer superlattices®®). One 
of the reasons for this is that the InGaAs layers 
are not completely ordered even in the <110> 
direction. Moreover, the ordered regions include 
many anti-phase boundaries. The 2DEG mobili- 
ty in the <110> direction could be improved by 
reducing the number of such boundaries. 


6. Conclusion 

We have reviewed recent advances in the 
application of TEM to the evaluation of struc- 
tures, interfaces and defects in IJ-V compound 
semiconductor thin films, solving various mate- 
rials issues in these systems. 

First, (GaAs),,/(GaP),, strained layer super- 
lattices grown by ALE have been structually 
evaluated by HRTEM at the atomic level. It 
has been shown microscopically for the first 
time that (GaAs),/(GaP), monolayer superlat- 
tices are formed uniformly throughout the 
epitaxial layer, indicating that this method 
achieves precise layer-by-layer growth. Sim- 
ilar results have also been obtained for 
(GaAs)3/(GaP)3; and (GaAs)4/(GaP), superlat- 
tices. 

AlAs/Si heterostructures grown by ALE 
have been evaluated by TEM. The hetero-inter- 
face is found to be extremely abrupt compared 
to GaAs/Si grown by a two-step MOVPE meth- 
od. Atomic steps are clearly observed at the 
interface and no amorphous or alloyed regions 
are found. At the initial stage of growth, a three- 
dimensional island growth mode is dominant. 
However, two-dimensional growth takes place 
after depositing only 36 atomic layers of AlAs. 
Two types of defect are observed in the AlAs 
layer, stacking faults and microtwins. Although 
the density of these defects is high, they tend to 
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self-annihilate at the interface between the AlAs 
layer and the thick, upper GaAs layer. 

The nature and behavior of APBs in GaAs 
and GaP on (001)Si grown by MOVPE have 
been investigated. It has been found that APDs 
are very complicated three-dimensional polygons 
with several sub-boundaries in different orienta- 
tions. Self-annihilation of APDs during growth 
is directly observed for the first time through 
plan-view and cross-sectional observations. Based 
on these results, a mechanism of annihilation of 
these domains is proposed. 

Interface defects in MBE-grown GaAsSb 
crystals and bulk defects in MOVPE-grown InP 
doped with Fe have been characterized. In MBE- 
grown GaAsSb crystal, we have found three 
types of interface defects, i.e. stacking fault 
tetrahedra, microtwins and dissociated disloca- 
tions. FeP precipitates are found in Fe-doped 
InP. They are found to be spherical coherent 
precipitates uniformly distributed in the crystal. 
Their density increases as the doping gas flow 
rate increases from 20 ml/min to 100 ml/min. 
In crystals grown with gas flow rates lower than 
5 ml/min, in which the Fe-concentration is 
below the solubility limit, no precipitates are 
found at all. For each defect, the mechanism of 
defect generation and possible ways to eliminate 
it were discussed. 

We have also described detailed TEM studies 
on atomic ordering in MOVPE-grown InGaP 
on (001)GaAs and MBE-grown InGaAs on 
(110)InP. From TED and HRTEM analyses, 
CuPt-type and CuAu-I-type ordered structures 
have been found in InGaP on (001)GaAs and 
InGaAs on (110)InP, respectively. In each case, 
the atomic structure of the ordered regions and 
dependence of the degree of ordering on growth 
conditions were described. The generation 
mechanisms of ordered structures and their in- 
fluence on the optical and electrical properties 
of the materials were also discussed. 
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Pulsed-Jet Epitaxy of Ill-V Compounds 
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(Manscript received October 14, 1991) 


Pulsed-jet epitaxy (PJE) for IIl-V compounds (GaAs, InP, GaP, and InAs) and its possibilities 
for process technology are described. A model involving selective reactions has been pro- 
posed to explain the self-limiting process in PJE. A reduction of carbon contamination has 
been achieved for GaAs growth using PJE and carrier concentrations ranging from 10!*cm? 


to 107°cm™* for n-type GaAs, and from 10'°cm™? to 10*!cm™*? for p-type GaAs. These 
can be obtained by controlling growth conditions and doping levels. The growth of 
(GaAs), (GaP), superlattices has demonstrated the high potential of PJE in superlattice 
growth. PJE offers unique possibilities for low temperature growth, selective growth, and 


uniform thickness growth. 


1. Introduction 

Atomic layer epitaxy (ALE) has many 
advantages over conventional growth methods 
such as molecular-beam epitaxy (MBE) and 
metalorganic vapor phase epitaxy (MOVPE)!)”®. 
The current interest in quantum-effect devices 
has led to improved crystal growth now con- 
trollable at the atomic level. Such progress in 
MOVPE and MBE techniques now makes it 
almost possible to grow epitaxial layers at the 
atomic level. In this type of growth, the layer 
thickness is a function of growth time, and in 
most cases the function is proportional to 
growth time. To produce layers of the right 
thickness, the growth rate is controlled based 
on data on the preceding growth. However, 
the growth temperature, gas flow (molecular 
flux) and partial pressure of source gases change 
with growth time and from run to run. Even 
though this change is small, it makes it difficult 
to control the growth thickness within a mono- 
layer level using conventional growth methods. 
At any one time, the source gas pressure (mo- 
lecular flux) and substrate temperature are also 
different at different points on the substrate 
wafer, even at the same gas flow rate. In MOVPE 
and MBE, much work has been done to improve 
the uniformity of the layer thickness. 
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ALE provides a solution to these problems. 
Firstly, it is a self-limiting mechanism. The 
epitaxial growth automatically stops at exactly 
one monolayer. This distinguishes ALE from 
other epitaxial techniques and enables the 
thickness of the layer to be controlled precisely 
at the atomic level. The ALE technique makes 
it possible to grow epitaxial layers with an 
extremely uniform thickness with little fluctua- 
tion on large substrates. 

Metalorganic ALE (MO-ALE), which uses 
metalorganic compounds for group III sources”, 
has several advantages in a flexible reactor 
design. It leads to fast exchange of group II 
and V source gases with a small gas manifold, it 
allows a number of highly pure metalorganic 
source to be used, and the techniques and 
concepts concerning the gas system and reactor 
developed for metalorganic vapor phase epitaxy 
can be applied. 

MO-ALE has problems, however, for control 
of the self-limiting process and in crystal purity. 
Temperature ranges in which self-limited growth 
is possible are reportedly very narrow for many 
Il-V compounds, and carbon contamination 
from metalorganic sources has been serious in 
GaAs ALE. 

A new MO-ALE technique called pulsed-jet 
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epitaxy (PJE) was developed to solve many of 
the problems of conventional ALEs for III-V 
compounds. PJE enables high purity II-V thin 
films with a clear self-limiting mechanism to be 
grown under a wide range of conditions. This 
paper describes the principle of PJE and possible 
applications of the technique to device pro- 
cesses. 


2. Self-limiting mechanism 

A source molecule (precursor) introduced 
into a reactor arrives at the growth surface 
after thermal decomposition and chemical 
reactions with other molecules in the gas phase. 
The arrival rate of the molecule, r,, per-second 
at a lattice site, i.e. adsorption site, is assumed 
to be proportional to the vapor pressure of the 
source molecule py, and can be written as 


Yar = KarPv, 


where K,,; is a constant proportional to the 
inverse square root of temperature T'/* , and 
derived from detailed balance and statics. 

The stability of a molecule at the adsorption 
site is determined from the total potential 
energy of the system including the adsorbed 
molecule and the adsorption site. When an 
adsorbed molecule in thermal equilibrium with 
the surface lattice desorbs from the surface in 
the gas phase, its thermal desorption rate, 
ras is written”? 


ras = Kas(T)vexp(—Eqs/kT), ou (2) 


where v is the thermal vibrational frequency 
of the adsorbed molecules, Fg, the dissociation 
energy (or band enthalpy), k the Boltzmann 
constant and Kg, the prefactor which involves 
the rotational and vibrational energies of the 
adsorbed molecule and is a function of tem- 
perature. 

In the ALE of a binary compound composed 
of cation A and anion B, let A®° be the source 
molecule for A (see Fig. 1). In some cases, 
A° may partly decompose or react with other 
chemicals in the gas phase. If the system can be 
described as a simple process where the source 
molecule A° is adsorbed (chemisorbed) at 
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Eya(B*) << E,a(A*) 
E,;(B*) >> E4.(A*) 


ee a 
VV AAA 


Fig. 1—Self-limiting process in ALE of crystal AB. 


surface site B* with a relatively large adsorption 
energy E,qg, then a new adsorption site A* will 
be formed instead of site B*. Molecule A° 
reacts chemically with adsorption site B*, 
and usually partly or completely decomposes 
on the surface. 

When the adsorption energy E,q of molecule 
A® is sufficiently large at the new adsorption 
site A*, A° does not adsorb onto the site. 
If the surface chemical species involving cation 
atom A has a small desorption energy Eg, at 
adsorption site A*, it will easily desorb from the 
surface. In both cases, the ‘effective’? adsorp- 
tion rate of molecule A®° appears to decrease 
as the surface density at site A* increases. 
This saturation effect of adsorption and desorp- 
tion may be used to make self-limitation possi- 
ble in ALE of compound AB. 

A typical example of the saturation effect 
in effective adsorption is where A* has a stable 
methyl group or a chlorine which prevents gas 
phase molecules from approaching the A®* site 
due to steric hindrance. However, the methyl 
group or chlorine is not necessarily needed for 
the self-limiting mechanism if site A*, which 
consists of only cation atom A, does not react 
with molecule A° (or has a large adsorption 
energy or a small desorption energy) at the 
growth temperature. 

When a molecule approaches the crystal 
surface, the chemical bond of the molecule is 
sometimes broken and new chemical bonds are 
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formed between the adsorbed molecule and 
surface atoms. This type of adsorption, or 
dissociative adsorption, is an important chemi- 
sorption process and it plays an important role 
in ALE growth. 

In dissociative adsorption, the valence 
electrons which form a chemical bond in the 
molecule transfer to the surface atom, while 
antibonding electrons transfer from the surface 
to the adsorbed molecule as the molecule 
approaches the surface. This weakens the 
strength of the chemical bond in the molecule. 
As the molecule approaches the surface, a 
strong chemical bond forms between surface 
atoms and atoms in the molecule, and finally 
the valence bond of the molecule is broken. 

Generally, the electron distribution, and 
therefore the bond length in the adsorbed 
molecule, charge, even though the change is 
not big as in the case of dissociative adsorption. 
Some molecules tend to decompose more easily 
on the surface than in the gas phase, and de- 
composition of the adsorbed molecule in ALE 
often cannot be distinguished from dissociative 
adsorption. 


3. Experiment 

Epitaxial growth was carried out by deposit- 
ing a monolayer of group III] atoms followed 
by a monolayer of group V atoms. The gas 
sources for the group III atoms were trimethyl- 
gallium (TMGa), triethylgallium (TEGa), tri- 
methylaluminum (TMAI) and _ trimethylindi- 
um (TMIn), and the gas sources for growth were 
contained in a low-pressure reactor designed 
so that thermal decomposition of the source 
gases was suppressed in the gas phase and 
proceeded only on the substrate surface. The 
reactor has a chimney-type structure where 
the carrier gas and reactants are introduced 
from the bottom and exhausted through the top 
by high capacity vacuum pumps. The reactor 
structure was determined by simulation to 
give smooth gas flow streamlines without 
vortices. The substrate was positioned in the fast 
gas stream emitted from a jet nozzle to prevent 
thermally decomposing source molecules build- 
ing up in stagnant layers. The exchange of 
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reactant gases was achieved using a pressure- 
balanced vent-and-run manifold. Alternating 
pulses of group II] metalorganic and group V 
hydride gases were separated by a hydrogen 
purge pulse which prevented the reactants from 
mixing. 

Growth thickness was estimated by measur- 
ing the step height between the epitaxial surface 
and a substrate whose surface was masked with 
SiO, which was removed before measurement. 
The growth rate was calculated by dividing the 
growth thickness by the gas cycle. A growth of 
500 cycles was usually used for the growth rate 
experiments. 


4. GaAs growth 
Figure 2 shows the variation of growth rate 
of GaAs(001) as a function of the TMGa pulse 


bh 


GaAs (001) O: 550°C O: 500°C 
@: 530°C m: 480°C 
A: 460°C 


o— 
lcycle 
oe 


TMGa | | | | | | 


AsH, 


Growth rate (layer/cycle) 


0 5 10 15 20 25 
TMGa pulse duration (s) 
Fig. 2—Dependence of GaAs growth rate as a function 
of the TMGa pulse duration. 


GaAs 
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Growth rate (layer/cycle) 
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TMGa flow rate, jrmoa (sccm) 

Fig. 3—Dependence of GaAs growth rate as a function 
of TMGa flow rate j/7MGa for different pulse 
durations. /~MGa corresponds to TMGa partial 
pressure. 


FUJITSU Sci. Tech. J., 28, 1, (March 1992) 


duration from 460 °C to 550°C. The growth 
rate is defined as the epitaxial thickness by the 
gas cycle and is expressed in monolayers/cycle. 
The thickness of one monolayer corresponds to 
0.283 nm (2.83 A) for GaAs(100). The rise time 
to one monolayer thickness decreases as growth 
temperature increases. This occurs because the 
decomposition rate of TMGa molecules adsorbed 
on the surface increases as the temperature 
increases. 

The GaAs growth rate was studied for 
different TMGa pulse duration as a function of 
TMGa flow rate jrmco, defined as the flow 
rate of hydrogen carrier gas through the TMGa 
cylinder (see Fig.3). TMGa_ partial pressure 
PTMGa 1S proportional to flow rate jpyoa, When 
the total flow rate of hydrogen carrier gas is 
constant in the reactor. The growth rate in- 
creases with jrmo, for different saturation 
thicknesses, then levels off. The saturation 
thickness continuously increases with TMGa 
pulse duration up to one monolayer (2.83 A). 
Figure 3. shows that the TMGa partial pres- 
sure corresponding to a TMGa flow rate above 
10 standard cm? min™! (sccm) is so high that 
the surface is completely covered with TMGa. 

The steric hindrance effect due to the 
methyl group is significant when the dissocia- 
tion energy Eg, of TMGa is very small at a sur- 
face site and TMGa therefore desorbs easily 
from the site at the growth temperature. After 
TMGa is supplied to the arsenic-coated surface, 
a TMGa molecule adsorbed on one site soon 
cannot occupy the nearest sites, because the 
molecule has a diameter of about 0.7 nm (7 A) 
which involves the van der Waals repulsive force. 
The surface density of the TMGa molecule 
cannot reach the density of — surface 
arsenic sites. 

Figure 3 shows that the saturation thickness 
is one monolayer for a TMGa pulse duration 
over 6s and no TMGa molecules exist on the 
surface. This suggests that the TMGa molecule 
is thermally decomposed or dissociatively 
adsorbed on the surface into di- or mono- 
methylgallium, or atomic gallium. This decom- 
posed state of TMGa is hereafter referred to as 
Ga*. 
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A rate equation can be used to describe the 
growth rate dependence in GaAs PJE. One of 
two things may happen to a TMGa molecule 
which arrives at the surface: one is desorption 
to the vapor phase, and the other is chemisorp- 
tion (dissociative adsorption) onto surface 
arsenic atoms where TMGa decomposes into 
dimethyl, monomethyl or atomic gallium. We 
assume that TMGa molecules cannot chemisorb 
onto the Ga*. 

Note that adsorbed TMGa quickly desorbs 
even on the arsenic atom around 500 °C. The 
desorption time constant Tg, was reported to 
be 0.6s for the TEGa molecule at 350 °C and 
much faster for TMGa®). This is consistent with 
the fact that, when the purge time following 
the TMGa pulse was changed from 0.1s to 
120s, no variations in the growth rate were 
observed. These results show that TMGa mole- 
cules cannot be effectively observed on the 
surface at 500 °C after a TMGa gas pulse. The 
surface density of adsorbed TMGa nypygc, and 
that of thermally decomposed methyl gallium 
or atomic gallium ng,* are calculated by the 
rate equations, 


dntmca/dt = Tar (1 — ntmGa —"%Ga*) 
= (Tes +ras)MTMGa> 
dng, */dt = res NTMGa: 


Neglecting the adsorbed TMGa at the end of 
the TMGa pulse, the growth thickness per gas 
cycle becomes ng,*. 

The temperature dependence of growth 
rate was quantitatively compared with the rate 
equations using ry, as the fitting parameter. 
razr is proportional to the mole fraction of 
TMGa, and a flow rate of 1 sccm through the 
TMGa bubbler corresponds to 5.6 x 10? s' for 
rar. The value of rgg was of the order of 10° s"! 
and its effect on the growth rate was small, 
because the growth carried out so that 7,, was 
much larger than rg,. This condition may be 
satisfied if the TMGa flow rate jpmcg, is above 
10 scem at 500 °C (see Fig. 2). The solid curves 
in Figs. 2 and 3 were calculated from the rate 
equations. The fitting parameter Tc, = in 8 
plotted as a function of growth temperature in 
Fig. 3, and an activation energy of 42 kcal/mol 
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Fig. 4—Comparison of XPS spectra between AsH3- 
exposed surface A and TMGa-exposed surface 
B. The dominant peak is due to the Auger 
LMM peak. 
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Fig. 5—Adsorption and desorption processes of TMGa 
molecule for GaAs (100) surface. 


is obtained for the dissociative chemisorption 
of TMGa on a GaAs(100) surface from 450 °C 
to 560°C. The saturation thickness below one 
monolayer for short pulse durations (see Fig. 2) 
is produced from the fast desorption of un- 
decomposed TMGa after a TMGa gas pulse. 
Surface adsorbates were identified using 
X-ray photoelectron spectroscopy (XPS) using 
an analysis chamber and growth chamber so that 
free movement from one chamber to another 
was possible out without exposure to air®?. 
Figure 4 compares XPS spectra between an 
AsH,-exposed surface A and a TMGa-exposed 
surface B. The GaAs(100) surface prepared by 
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GaAs (100) 
500°C 


Growth rate (layer/cycle) 


TEGa pulse duration (s) 
Fig. 6—Growth rate dependence of GaAs layer as a 
function of TEGa pulse duration (pulsed-jet 
epitaxy). 


ALE growth was exposed to AsH3; or TMGa 
vapor at 500 °C and cooled to room tempera- 
ture for the XPS analysis. The observed peak 
was not due to XPS, but was due to a Ga LMM 
Auger emission. If there are carbon or carbon- 
related bonds on the surface, XPS emission 
associated with carbon should be observed at 
285 eV on the higher energy side of the GaLMM 
Auger peak. This shows that there were few 
methyl radicals on the TMGa exposed surface. 
Figure 5 shows the adsorption and desorp- 
tion of TMGa on the arsenic-coated (100)GaAs 
surface, which explains the atomic Ga model®?, 
TMGa dissociatively adsorbs onto the GaAs 
surface and decomposes into monomethylgalli- 
um or dimethylgallium. These methylgalliums 
may desorb from the surface (the surface life- 
time is very short, in the order of milliseconds), 
but some methylgalliums further decompose 
into Ga atoms (activation energy 42 kcal/mol). 
TEGa is a typical gallium source, which 
decomposes at a lower temperature compared 
with TMGa and easily decreases the carbon 
contamination of the GaAs layer in conventional 
MOVPE or MOMBE growth. Figure 6 shows the 
growth rate variation in the GaAs layer as a 
function of the TEGa pulse duration in pulsed- 
jet epitaxy. Single monolayer saturation tended 
to occur for shorter pulse durations, but the 
growth rate increases beyond one monolayer 
when the pulse duration is above 7s. A large 
fluctuation in the thickness profile was also 
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observed for an epitaxial layer. Error bars repre- 
sent thickness variations within the epitaxial 
layer. When the growth rate increased beyond 
one monolayer of GaAs, the surface morphology 
of the layer changed drastically from a mirror- 
like to a rough surface. The surface was covered 
with gallium droplets when the growth rate was 
above one monolayer. Growth for the self- 
limiting mechanism was studied using TEGa 
under various growth conditions, which involved 
low temperature growth below 400 °C, but no 
self-limiting was observed in the growth rate. 
However, there is no definite model to explain 
this. Given that the vapor pressure of source 
molecules is closely related to the self-limiting 
mechanism, these phenomena may be produced 
more easily for source materials with a higher 
vapor pressure. 

Tri-isobutylgallium (TIBGa), which is liquid 
at room temperature, was studied as a possible 
gallium source for low temperature ALE. TIBGa 
molecules decompose on the As-coated surface 
even at 200 °C. The vapor pressure of TIBGa is 
200 Pa at 54°C, thus it can be used in the 
conventional organometallic gas system. A GaAs 
layer was grown using TIBGa and AsH; in the 
pulsed-jet reactor. When growth was at 400 °C, 
the thickness per cycle increased with TIBGa 
pulse duration. Below 350 °C, growth was partly 
self-limited and saturated at one monolayer. 
GaAs layers with excellent surface morphology 
were grown by ALE using TIBGa and AsH,, 


GaAs (100) 
Ts: 300°C 
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Fig. 7—Growth rate of GaAs as a function of AsH3 pulse 
duration. 
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even at 250 °C. 

In ALE, there has been little focus on the 
supply condition of AsH3 because of the high 
As vapor pressure and the weak bonding be- 
tween As atoms at the growth temperature. 
Below 350 °C, the growth rate strongly depends 
on the AsH3 supply. Figure 7 shows the growth 
thickness per cycle (growth rate) as a function 
of AsH,; pulse duration. After a quick rise in the 
growth rate at pulse durations below 1 s, the 
growth rate decreased with AsH, pulse duration. 
The surface morphology was very smooth and 
mirror-like for longer AsH; pulse durations. The 
same result was obtained when the AsH; partial 
pressure was changed, i.e. increasing the AsH; 
pressure decreased the growth rate. Another 
interesting feature is that the growth rate 
decreases for incremented steps of one mono- 
layer of GaAs. This indicates that TIBGa is very 
sensitive to the atomic structure of the arsenic- 
coated surface. 


5. InP, InAs and GaP 

Figure 8 shows the dependence of the 
growth rate of InP on the TMIn pulse duration. 
At 600°C, the thickness per cycle increased 
beyond one monolayer/cycle as pulse duration 
increased. At 350 °C, the growth rate remained 
constant for a pulse duration above 1s. The 
saturation thickness was, however, 0.5 mono- 
layers on the (100) face. A similar dependence 
was observed as the TMIn mole fraction was 


Growth rate (layer/cycle) 


TMIn pulse duration (s) 


Fig. 8—Growth rate of dependence of InP as a function 
of the TMIn pulse duration. 
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varied. The thickness per cycle increased with 
TMIn mole fraction beyond 1 monolayer above 
400°C, while the thickness saturated at 
0.5 monolayer below 400 °C!”. 

XPS studies were made for the InP surface. 
No methylindium was observed on the surface 
around 350°C as the GaAs surface. This sug- 
gests that the self-limiting process is due to 
selective adsorption and desorption on indium 
and phosphorus atoms, as in the pulsed-jet 
epitaxy of GaAs. The saturation thickness of 
about half a monolayer is probably due to a 
stable surface reconstruction of InP around 
350 °C. 

InAs and GaP ALE was studied using PJE. 
Figure 9 shows the growth rate dependence of 
GaP(100) as a function of the TMGa flow rate at 
500°C. A clear self-limiting condition was 
observed the growth rate saturated at just one 


GaP 500°C 
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Fig. 9—Growth rate dependence of GaP on the TMGa 
flow rate. 
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Fig. 10—Temperature dependence of saturation thick- 
nesses for InAs, GaP and GaAs by PJE. 
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monolayer when the pulse duration was higher 
than 5s with increased TMGa flow rate. A self- 
limiting condition was also observed for the 
TMGa pulse duration dependence. Similar re- 
sults were obtained for InAs ALE, where the 
growth rate saturated at just one monolayer 
with increased TMIn pulse duration and partial 
pressure. Figure 10 shows the temperature 
dependence of the saturation thickness for InAs, 
GaP, and GaAs. One monolayer saturation 
occurs over wide temperature ranges for these 
compounds’ a 


6. Purity and impurity doping 
Purification and impurity doping techniques 
have been developed for PJE'?!®). GaAs 
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Fig. 11—Carrier (hole) density variation with TMGa pulse 
duration. 
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Fig. 12—Carrier (hole) density variation with AsH3 pulse 
duration. 
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epitaxial growth under normal conditions 
exhibited p-type conductivity due to carbon 
acceptors. Figures 11 and 12 show the carrier 
(hole) densities of GaAs layers as a function of 
TMGa and AsH3 pulse duration in PJE. The 
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Fig. 13—Photoluminescence spectra for GaAs layers 
A: 550°C, TMGa pulse duration of 0.8 s, 
B: 500°C, TMGa pulse duration of 20 s. 
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carrier density decreases drastically as TMGa 
pulse duration decreases or AsH, pulse duration 
increases. High-temperature growth above 
550°C can achieve one-monolayer saturation 
for a short TMGa pulse for a duration below 1 s 
(see Fig. 2). This allows highly pure GaAs layers 
to be grown by PJE. Figure 13 compares 4.2 K 
photoluminescence spectra for growth at 500 °C 
and TMGa pulse duration of 20s, and 550 °C, 
and TMGa pulse duration of 0.8s. The first 
sample was a_ p-type (hole density: 
3 x 108 cmv?) and it showed broad lumines- 
cence at 830 nm due to donor-acceptor (carbon 
acceptor) pair recombination. The second 
sample was an n-type (carrier density: 
3x 10!4 cm). It exhibited sharp emission 
peaks due to the free and bound excitation 
peaks near the band gap energy but no emissions 
associated with the carbon acceptor, which 
reflects the high purity of the epitaxial layer. 

It was possible to adjust the carrier density 
from 3 x 10!4 cnr? to 1 x 10?! cn‘? for p-type 
GaAs by varying the growth conditions, and 
from 1x10 ex? to 210" cx? for 
n-type by doping control using H, Se, resulting 
in good surface morphology (see Fig. 14). These 
show that doping control in PJE is sufficient for 
growing GaAs layers for conventional devices. 
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Fig. 14—Carrier concentrations of ALE-GaAs. 
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7. Selective growth and thickness uniformity 

In ALE, ideal selective growth is possible, 
since growth proceeds  two-dimensionally. 
Figure 15 compares the thickness variation of 
GaAs epitaxial layers from the SiO, mask edge 
for MOVPE and ALE growth. The epitaxial 
thickness of ALE shows the same surface level, 
while that of MOVPE gradually deceases for the 
SiO, edge because of faceted or lateral growth 
at the edge. It was also found that growth on 
stripes narrower than | wm was with a complete- 
ly flat surface with no abnormal growth. 
Epitaxial patterns had perfectly mirror-like 
surfaces, although the substrate surface under- 
went several processes. This is important for the 
precise masking process using planar techniques. 

It is easy for ALE to grow an epitaxial layer 
with an extremely uniform thickness. Figure 16 
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Fig. 15— Thickness variations of MOVPE and ALE (PJE) 
grown GaAs epitaxials from SiO, edge. 
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Fig. 16—Thickness variations of MOVPE and ALE (PJE) 
grown GaAs epitaxials on 2-inch substrates. 
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shows the thickness variations of ALE and 
MOVPE epitaxial layers on 2-inch wafers. Both 
epitaxial layers were grown in the same 
chamber, without any rotation of the GaAs 
substrate. The thickness variation for ALE 
grown layer was within 1.5 percent which is 
comparable to error in measuring thickness. 


8. Heteroepitaxy and superlattices 

Heteroepitaxy requires a much more com- 
plex procedure for ALE than homoepitaxy. The 
growth mechanism of the first layer for 
compound AB is not the same as that of the 
homoepitaxy for AB, because desorption, migra- 
tion, and decomposition of precursors on the 
surface differ for crystal AB and substrate A'B’. 
Interface stress due to lattice mismatch greatly 
affects the growth mode of epitaxy and some- 
times compromises the self-limiting mechanism. 

The growth of strained-layer superlattices by 
ALE is an important test to examine the 
potential of ALE, because many complex proce- 
dures are involved. We grew (GaAs), (GaP), 
strained-layered superlattices using PJE!”, 
because both crystals are self-limited at just one 
monolayer in the same temperature region, 
around 500 °C (see Fig. 10). 

In ALE, the partial pressure and pulse dura- 
tion of the sources were changed depending on 
the nature of the top surface. When the surface 
was coated with arsenic (phosphorus), TMGa 
was supplied to the surface under conditions 


(GaAs), (GaP),/GaAs 
500°C 


Growth rate (layer/cycle) 


TMGa pulse duration (s) 
Fig. 17—Growth rate dependence of monolayer 
superlattice (GaAs), (GaP), as a function of 
TMGa pulse duration after a PH3 pulse. 
One monolayer is equal to 5.55 A (0.555 nm). 


FUJITSU Sci. Tech. J., 28, 1, (March 1992) 


(GaAs)m(GaP), 


(m,n) 


Experimental thickness (nm) 


0 1 


Calculated thickness (nm) 


i) 


GaAs 


v4 


+nxX 


cap 
mX GaP 
2 


Fig. 18—Comparison of one-period experiment and 
calculated thicknesses for various superlattices 
(GaAs) (GaP), . 


similar to GaAs (GaP) homoepitaxy. AsH3; and 
PH; were supplied with a pulse duration and 
partial pressure high enough to attain one 
monolayer for GaAs and GaP homoepitaxies. 

Figure 17 shows the growth rate of one 
period of (GaAs),(GaP), for different pulse 
durations of TMGa onto the arsenic-coated 
surface. One monolayer is defined as the sum of 
a GaAs and GaP single-monolayer thickness and 
is equal to 0.555 nm (5.55 A). A self-limiting 
mechanism was maintained for growth of the 
strained-layered superlattice. A self-limiting 
condition in the growth rate was also observed 
for variations in other growth parameters (TMGa 
partial pressure, ASH, pulse duration and partial 
pressure) involving the gallium-coated surface. 
Figure 18 compares the calculated and experi- 
mental one-period thickness of (GaAs),,, (GaP), 
superlattices. The calculated thickness is defined 
as m x (GaAs lattice constant)/2 +n x (GaP lat- 
tice constant)/2 according to Vegard’s law. The 
two sets of values agree well, indicating that the 
growth of each layer is self-limited to one mono- 
layer. 

We analyzed layer structures for superlat- 
tices using X-ray rocking curves and Raman 
scattering measurements. (GaAs), (GaP), , 
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Fig. 19—X-ray rocking curves for (GaAs), (GaP), 
(m=n=1, 2,3). 
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Fig. 20—4.2 K photoluminescence spectrum for 
(GaAs), (GaP). 


(GaAs) (GaP), , and (GaAs)3(GaP)3; show many 
diffraction peaks in addition to the (400) and 
(200) peaks which are allowed for the 
GaAso.5Po.5 alloy (see Fig. 19). Their diffraction 
angles agree precisely with the calculated ones, 
as do their intensities. Raman scattering showed 
new vibrational modes resulting from Brillouin 
zone folding at the zone center for 
(GaAs),,(GaP),. The periodicity of the 
superlattice causes Brillouin zone folding of the 
acoustic phonon branch and the confinement of 
optical phonons. The observed peak energies 
agree with those calculated when based on an 
elastic continuum model!®?. 
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Electronic properties were studied by photo- 
luminescence and photo-reflectance measure- 
ments!°?, Figure 20 shows a photoluminescence 
spectrum for the monolayer superlattice 
(GaAs); (GaP),. A strong luminescence peak 
occurs at about 700 nm with a relatively narrow 
half-width. Other superlattices (m, n = 2, 3,4, 5) 
and a GaAsg5Po9.5 alloy exhibit no such sharp 
luminescence. Photo-reflectance measurements 
also showed that the band-gap energy abruptly 
decreased with layer number from m = n= 2 to 
1, although the band-gap energy slowly 
increased as the layer number decreased 
above m, n > 2. These findings suggest that 
(GaAs); (GaP); has a direct band-gap ,energy 
differing from other superlattices and 
GaAso 5Po.5. 


9. Conclusion 

We showed PJE to be effective in ALE for 
Ill-V compounds. A self-limiting mechanism in 
PJE is explained by dissociative adsorption of 
TMGa molecules onto surface arsenic atoms. 
Atomic gallium may be responsible for the self- 
limiting condition in ALE. The activation energy 
of TMGa decomposition on the GaAs surface 
was found to be 42 kcal/mol around 500 °C, 
which is smaller than that in the gas phase 
because of the surface catalytic effect. Very 
pure GaAs films were produced, together with 
a wide range of carrier concentrations for 
n-type and p-typed doping of GaAs layers. 
Selective growth on fine patterns and uniform 
growth on large sugstrates show great promise 
for process techniques. Optical studies by photo- 
luminescence and reflectance measurements of 
strained-layered superlattices (GaAs), (GaP),, 
grown on GaAs substrates by PJE suggest 
that the monolayer superlattice (GaAs), (GaP), 
had a direct-band gap structure. 
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Safer MOVPE growth of GaAs, AlGaAs, and InGaAs on 3-inch GaAs substrates was investigat- 


ed using organic group-V sources. 


Tertiarybutylarsine (TBAs) and monoethylarsine (EtAs) were used with a horizontal atmos- 
pheric pressure reactor. GaAs and AlGaAs grown using TBAs were of high purity and extremely 
uniform, equivalent to those using arsine. With either silane or disilane doping, Si incorpora- 
tion efficiency using TBAs or EtAs was always higher than that using arsine. The group-V 
pyrolysis simulation implies that H,AsSiH3 dominates in the gas phase when using TBAs, 
corresponding to high doping efficiency. High quality, high electron-mobility transistors were 
achieved with TBAs. These results verify that epitaxial layers grown using TBAs are of suffi- 


ciently high quality for device applications. 


1. Introduction 

Metalorganic vapor phase epitaxy (MOVPE) 
is a key technique in growing large-area III-V 
semiconductor materials for high-speed and 
optoelectronic devices, because it provides 
highly uniform, top quality materials having 
abrupt compositional changes on an atomic 
monolayer scale. The drawback is the danger 
posed by the use of highly toxic group-V hy- 
drides, i.e. arsine (AsH;3) and phosphine (PH; ). 
Since they are stored in high-pressure cylinders, 
careful handling, especially when changing 
cylinders is required to avoid dangerous leakage. 
Safety apparatus, i.e. special exhaust equipment, 
is also required in case of accidental release. 

Given the danger involved, many studies 
have been made in the search for less toxic 
As and P sources alternative to arsine and 
phosphine, as Stringfellow has reviewed!?. Re- 
cently, less hazardous organic group-V pre- 
cursors have attracted much attention. They 
are less toxic than arsine and phosphine, and 
are liquid at room temperature. Should a large 
leakage occur, the danger of exposure should 
be limited due to the lower toxicity and lower 
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cylinder pressure. Early studies with trimethyl- 
arsine (TMAs) were hampered by significant 
impurity in the source materials”, Using 
'3C-labeled TMAs, Lum et al. showed that 
carbon from TMAs was incorporated at a 
concentration of 5 x 10!® cm in GaAs layers”). 
The carbon background concentration could be 
reduced by using a graphite thermal precracker 
for TMAs*). Some researchers used TMAs with 
arsine as a carbon dopant source for the growth 
of heavily p-type GaAs°?. Using triethylarsine 
(TEAs), GaAs layers were contaminated with 
carbon in the 10!7 cm™3 range’)"?), Using layers 
grown with TEAs, Fujita et al. fabricated a 
double-heterojunction laser diode that operated 
in pulse mode at room temperature!” . 

In contrast with fully substituted organo- 
arsines, i.e. TMAs and TEAs, partially hydrogen- 
substituted organoarsines have given more 
promising results. High quality GaAs layers were 
obtained using diethylarsine (DEAs)'!), phenyl- 
arsine (PhAs)!??, monoethylarsine (EtAs)3)14) 
and__ tertiarybutylarsine (TBAs)!*) 29), Using 
TBAs, Haacke et al. obtained the purest GaAs 
layers among non-arsine group-V sources?*)»?®) 
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The GaAs layers had an electron mobility of 
160000 cm? V's"! at 77 K. These results verify 
that partially substituted organoarsines, as As 
sources safer than arsine, are promising candi- 
dates for obtaining high-quality GaAs films with 
negligible carbon incorporation due to the effect 
of the As-H bond?”?. However, there has been 
little research into the growth of AlGaAs and 
InGaAs on GaAs substrate under various condi- 
tions, or into uniformity across a 3-inch sub- 
strate, due to limited quantities of organoarsine. 

Several researchers have reported differences 
in Si-doping characteristics in GaAs between 
TBAs and _ arsine?®)-3?). However, the mecha- 
nism responsible for these differences in doping 
efficiency was unknown. With arsine, many 
studies have been reported using silane (SiH, ) 
and disilane (Si, H, )°*?”*”?. 

Experimental device fabrication using a 
hydrogen-substituted organoarsine is essential 
for developing MOVPE manufacturing. Many 
successful results were achieved with TBAs, such 
as GaAs metal semiconductor field effect 
transistors (MESFETs)°°">>»5*), p-i-n photo- 
diodes?) 55),56) lasers®’)“5®), and solar cells”, 
Recently, Sundaram et al.°!) reported GaAs 
MESFETs exhibiting high performance. How- 
ever, no high-speed devices consisting of both 
AlGaAs and GaAs have yet been fabricated by 
MOVPE using TBAs. 

In this paper, we describe the use of TBAs as 
an arsenic source for the growth of high-quality 
uniform GaAs, AlGaAs, InGaAs, and selectively 
doped heterostructures on 3-inch GaAs sub- 
strates in an atmospheric pressure reactor under 
various conditions, and compare the properties 
of epitaxial layers grown using TBAs with those 
grown using arsine?>»?*) EtAs was also used, 
and appears to be another promising source for 
GaAs growth??. We investigated silane and 
disilane doping of GaAs with organoarsine and 
arsine*!) to clarify the role of As-precursors, 
because doping control is indispensable to the 
reproducible production of high-performance 
devices. Furthermore, we demonstrated effec- 
tive planar doping of Si into GaAs with TBAs 
to eliminate the effect of group-II sources*”), 
The comparison between our results and pre- 
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vious reports is addressed. We applied a kinetic 
simulation on the pyrolysis of the group-V 
sources, using the results to propose a precise 
Si-doping model. Electron traps in GaAs grown 
using TBAs were measured by deep level tran- 
sient spectroscopy (DLTS)°”), because the 
stoichiometric deviation in a crystal also pro- 
vides useful information about the type of 
reactant on a substrate, and deep levels, such as 
EL2, are strongly affected by stoichiometric 
deviation. 

We fabricated n-AlGaAs/GaAs high electron 
mobility transistors (HEMTs) and n-AlGaAs/ 
InGaAs/GaAs pseudomorphic HEMTs_ with 
wafers grown using TBAs°*. The growth of a 
selectively doped heterostructure®*) and the 
fabrication of HEMTs are landmarks in feasibility 
checks for device applications. High quality 
GaAs and AlGaAs films and abrupt composition- 
al transitions are required to obtain high quality 
two-dimensional electron gas (2DEG) at the 
heterointerface°*?. HEMTs are potentially su- 
perior to GaAs MESFETs due to the superior 
electron transport in HEMT channels and the 
unique electrical properties of the HEMT 
structure, which results in high-speed _per- 
formance and low power dissipation®®?. Further- 
more, the HEMT structure suppresses the short 
channel effect®®’. The n-AlGaAs/InGaAs/GaAs 
pseudomorphic HEMTs are more promising than 
n-AlGaAs/GaAs conventional HEMTs because 
they perform better in dc and microwave 
applications, due to a higher sheet carrier 
concentration and higher saturated velocity in 
the InGaAs channel layer® 768) | 

The natural convection induced by buoyancy 
which occurs in all conventional horizontal 
atomspheric pressure reactors promoted deposi- 
tion onto the ceiling, degrading film uniformity. 
We designed an inverted horizontal atmospheric 
pressure reactor to demonstrate the growth of 
highly-uniform epitaxial layers. The inverted 
reactor eliminated buoyancy-induced flow, 
reducing source consumption due to _ wall 
deposition. 


2. Experimental 
We used two types of horizontal atmos- 
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Fig. 1—Schematic drawing indicating the expected flow 
line?) 


pheric pressure reactor using rectangular cross- 
section fused quartz of 10 cm to 14 cm wide by 
2 cm tall as shown in Fig. 1. The first, Fig. 1a), 
was a common conventional reactor’). The 
other, Fig. 1b), was an inverted reactor we 
developed in this study. The stream lines in 
Fig. 1 schematically represent the gas flow we 
speculated. 

Epitaxial layers were grown on stationary or 
rotated 2- or 3-inch GaAs substrates oriented 
2.5° off (100) toward <110>, which were 
brought to growth temperature from 425 °C 
to 700 °C by radio frequency (RF) induction. 
Group-III sources were trimethylgallium (TMGa), 
trimethylaluminum (TMAI), and trimethylindi- 
um (TMIn). Organometallic precursors were 
contained in stainless steel bubblers in cooling 
baths. The TBAs bubbler (American Cyanamid 
Company) was kept at 4°C or 10°C with a 
vapor pressure of 63 Torr or 96 Torr’®’. Those 
for EtAs (Morton International) were at 10 °C 
and 311 Torr’). Three different lots of TBAs 
were used in this section, to investigate growth 
reproducibility from batch to batch of TBAs?”. 
Lot A, lot B and lot C were produced in 1987, 
1989, and 1990. The group-V flow rate ranged 
from 4x 1074 mol/min to 6 x 10°? mol/min. 
GaAs, Alp »gGda9.72 As, and Ing ;s;Gap.g5As were 
grown with V/II ratios of 5 to 40, 20 to 120, 
and 2.5 to 90, respectively. Typical growth 
rates at 650°C for GaAs, Alo gGag.72As, and 
Ing.15Gag.gsAsS were 1lyum/h to 3.6 um/h, 
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1.5um/h, and 1.2 wm/h, respectively. Silane or 
disilane was used as the n-type dopant source. 
Hydrogen gas diffused through a palladium cell 
was used as the carrier gas. All species were 
mixed just before putting them into the reactor. 
The flow velocity was varied from 8.5 cm/s to 
30 cm/s. 

The surface morphology was observed by 
Normarski interference contrast optical micro- 
scopy, and by a surface contamination analyzer 
(Surfscan 4500, Tencor Instruments) that could 
count the density of irregular particles on the 
surface as a function of the laser scattering cross 
section. The donor concentration (Ng) was 
determined using the capacitance-voltage (C-V) 
and van der Pauw-Hall methods. GaAs growth 
rates were determined using a stylus profiler 
to measure GaAs thickness after selectively 
removing an unmasked GaAs layer on an 
Alp.65Gao.35As stopper layer using citric acid’?), 
AlGaAs growth rates were determined by the 
C-V method while observing the depth of the 
peak which appeared at the heterointerface of 
the n-AlGaAs/n-GaAs structures. The optical 
properties of the epitaxial layers were evaluated 
by photoluminescence (PL) spectroscopy per- 
formed at 77K and 4.2 K. The samples were 
excited using multiline emissions of an argon ion 
laser and PL signals were detected with a Si p-i-n 
photodetector or a photomultiplier. Electron 
traps were measured by the DLTS method’*? 
from —180°C to +150°C. A double-crystal 
X-ray diffractometer was also used to evaluate 
the quality of the AlGaAs. Both the 77 K PL 
and X-ray diffractometer measurements were 
performed to determine the AlAs mole fraction 
of the Al,Ga,., As. We used the results of Baliga 
et at.’ for the relationship between band 
gap energy and the InAs mole fraction for 
In,Ga,.,As growth (y <0.2). The residual 
impurities in the epitaxial layers were charac- 
terized by secondary ion mass spectroscopy 
(SIMS) using Cs* ions as the primary ion source. 
The self-compensation of Si into GaAs was 
estimated using PL measurements at 77 K, 
observing the self-activated (SA) center peak 
appearing around | pm75), 
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3. Results and discussion 
3.1 Undoped GaAs 

First, we investigate carbon incorporation 
and impurity contamination in undoped GaAs. 
Differences in growth rate due to different 
growth mechanisms between TBAs and arsine 
will be addressed. 

The surface morphology of GaAs grown 
using TBAs was excellent, as good as that grown 
using arsine. There was no difference among the 
TBAs lots. Figure 2 shows the growth rates of 
the GaAs epitaxial layers for TBAs and arsine as 
a function of growth temperature (7,). The 
growth rate using TBAs increases with increasing 
temperature below 550°C, but is fairly inde- 
pendent of temperature between 550°C and 
650°C. The behavior at temperatures below 
650 °C is similar to that observed using arsine. 
In the former region (below 550 °C), called the 
kinetically controlled region, the growth rate is 
controlled by the decomposition rate of group- 
Ill species. In the latter region (between 550 °C 
and 650°C), called the transport-controlled 
region, it is controlled by the diffusion of the 
decomposed metalorganic components’®). The 
activation energies of the growth rates observed 
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Fig. 3—Growth rates as a function of the V/III ratio for 
GaAs layers grown at 475°C using TBAs and 
arsine. 


in the kinetically controlled region were 145 kJ/ 
mol using arsine. The activation energies for 
both sources are extremely close, although the 
value of the growth rate using TBAs in the 
kinetically controlled region is about twice 
that of the rate using arsine. Above 650 °C, the 
growth rate for TBAs decreases slightly with 
increasing temperature. In contrast, the growth 
rate for arsine is independent of temperature. 
There was no difference in growth rate among 
the TBAs lots. 

The growth rate dependence on the V/II 
ratio at 475 °C is shown in Fig. 3. At 475 °C, 
the growth rate using TBAs increases with 
increasing V/III ratio. If this is due only to the 
difference in decomposition temperatures of 
TBAs and arsine!’)"!9?”), the growth rate at 
a higher V/III ratio should saturate at the same 
value for both TBAs and arsine. Since these 
values seem to be different, as shown in Fig. 3, 
we suppose that there must be a different 
growth mechanism for each source under these 
growth conditions. At low temperatures, the 
reaction between TMGa and arsine takes place 
heterogeneously and may proceed via_ the 
Langmuir-Hinshelwood mechanism’®. If this 
is the case, we must consider the reversible 
excluded-volume chemisorption mechanism 
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Fig. 4—Growth rates for GaAs films. Growth tempera- 
ture is 600 °C and V/II ratio is 109), 


presented by Aspnes et al’). In this mechanism, 
an undecomposed TMGa molecule chemisorbed 
at a specific lattice site on an initially arsine- 
stabilized surface hides several adjacent sites, 
causing decreased growth efficiency. TBAs is 
expected to be a stronger adduct former than 
arsine due to the electron-donating nature of the 
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tertiarybutyl substitute. Therefore, the TMGa- 
TBAs adducts chemisorb on the surface. Either 
the reversible excluded-volume mechanism does 
not play an important role in the growth, or the 
excluded volume of the TMGa-TBAs adducts is 
smaller than that of TMGa due to the deforma- 
tion of TMGa from a planar to a pyramidal 
structure with adduct formation. We believe 
that the difference in growth rates of both 
sources should be observable at 475 °C. 

In addition, we obtained two significant 
results pertinent to GaAs growth rates at 600 °C. 
First, the growth rate using arsine is a linear 
function of the TMGa flow rate, as shown in 
Fig. 4a). On the contrary, the growth rate using 
TBAs saturates at a high TMGa flow rate over 
5 x 10° mol/min. Second, growth rates depend 
strongly on the total gas flow rate at 600 °C, and 
different for each source, as shown in Fig. 4b). 
The growth rate for TBAs varies in inverse pro- 
portion to flow rate, while that for arsine varies 
inversely with the square root of the flow rate. 
It is remarkable that the growth rates for higher 
total flow rates are the same for TBAs and arsine. 
This suggests that the species residence time for 
gas flow over the wafer has a significant effect 
on GaAs growth. According to the boundary 
layer theory, the growth rate should be approxi- 
mately proportional to v~'/?, where vu is the 
velocity of total gas flow®®’. The growth rate for 
arsine can only be predicted theoretically. We 
speculate that gas phase reactions, i.e. adduct 
formation and homogeneous nucleations among 
metalorganic components, play an important 
role in the growth of GaAs using TBAs. Lower 
TMGa flow rates and higher total gas flow rates 
reduce the gas phase reaction between TMGa 
and TBAs. In addition, the results of growth 
uniformity agree with this interpretation. We 
discuss AlGaAs growth uniformity in sections 
3.3 and 3.4. 

Figure 5 shows the free carrier concentration 
and electron mobility at 77 K, using lot A, as a 
function of growth temperature for a V/III ratio 
of 10. The epitaxial layers showed n-type 
conductivity under all growth conditions. The 
best crystal was obtained at 700 °C. The highest 
electron mobilities were 6700cm*V!s! at 
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Fig. 5—Free carrier concentration at 300 K and electron 
mobility at 77 K of GaAs films using lot A with 


a V/IIl ratio of 10 as a function of growth tem- 
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Fig. 6—Photoluminescence (PL) spectra in the exiton 
region measured at 4.2 K for GaAs layers grown 
at 700 °C with a V/III ratio of 10. 


c) Arsine 


room temperature and 30 000 cm? V"!s! at 77 K 
with a carrier concentration of 3 x 10!5 cm™? 
at room temperature for the 2 um-thick GaAs 
layer. The carrier concentration decreases and 
the electron mobility increases with increasing 
growth temperature. The residual dopant was 
identified as sulfur only by SIMS analysis. The 
compensation ratio, defined as the ratio of 
ionized acceptor concentration to ionized donor 
concentration, was assessed as 0,3 at growth tem- 
peratures of 600°C to 700°C by Hall measure- 
ments®!). This ratio for films grown using arsine 
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Fig. 7—PL spectra at 4.2 K for GaAs layers. 


with an electron mobility of 120000 cm? V7!s"! 
at 77 K was about 0.28. 

The 8 wm-thick epitaxial layer grown at 
650°C with a V/III ratio of 10, using lot B 
TBAs, showed p-type conductivity with hole 
mobilities of 400 cm?V~!s7! at room temper- 
ature and 7600cm?V's!?! at 77K, with a 
carrier concentration of 3 x 10!4 cm™3 at room 
temperature. The other epitaxial layers grown at 
temperatures between 600 °C and 700 °C were 
2pum thick and were too thin for electrical 
property measurements, indicating high purity 
films. 

Figure 6 shows PL spectra in the exciton 
region measured at 4.2 K for films using lot A, 
lot B, and arsine, and grown under the same 
conditions. The various sharp lines which appear 
just below the band gap correspond to the 
recombination of free and bound excitons. The 
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spectrum for lot A is insufficient compared with 
that of arsine due to the contamination of lot A 
with sulfur. The excition peaks of the lot B 
film with high resistivity are as well resolved as 
those of the arsine film with an electron mobil- 
ity of 120000 cm?V"!s7? at 77K. The emis- 
sions from free excitons (FEs) for the lot B 
layer are very strong. These results show that the 
lot B GaAs epitaxial layers are of high purity. 

The best PL spectrum using lot C which 
contains acceptor-related peaks is shown in 
Fig. 7a). This film was grown at 600 °C with a 
V/IIl ratio of 20. The intensity of carbon- 
acceptor-related peaks are comparable to that of 
exciton peaks. The full-width at half-maximum 
(FWHM) of the emission corresponding to band 
to acceptor (e, C°) transition is 1.7 meV for the 
films using TBAs and 1.5 meV for those using 
arsine. The Zn detected by PL comes from the 
TMGa. This is known because Zn also remains 
in GaAs films grown using arsine. Therefore, we 
conclude that lot B is much purer than lot A and 
carbon incorporation into GaAs grown using 
TBAs is negligible. As well as lot B, GaAs 
layers grown using lot C (the most recent 
batch) exhibit high quality with excellent 
reproducibility from run to run. 

GaAs layers were grown using EtAs at 
650°C. The surface morphology was good. 
The 5 um epitaxial layer grown with a V/II 
ratio of 10 had p-type conductivity showing a 
hole concentration of 2.5 x 10!4 cm? with a 
hole mobility of 400 cm?V-!s at 300 K and 
5700 cm?V~!s* at 77 K, and that grown with 
a V/IIl ratio of 20 showed high resistivity. 
Figure 7b) shows the PL spectrum at 4.2 K for 
the film with a V/III ratio of 20. Exciton and 
acceptor peaks were well resolved, suggesting 
that highly pure GaAs equivalent to arsine- 
grown and TBAs-grown GaAs were obtained 
using EtAs. Through the carrier concentration 
depends on the V/III ratio, we suppose that a 
higher V/III ratio will reduce the amount of 
carbon incorporation. 


3.2 Si-doped GaAs 
The purity of undoped GaAs is not the 
only important factor for device fabrication. 
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Controllability of n-type doping of GaAs is 
also required. In this section, we describe 
detailed doping research using silane or disilane 
as a Si dopant source. Differences in Si doping 
efficiencly between TBAs and arsine will be 
discussed, proposing a doping mechanism using 
simulation of group-V source decomposition. 
Only pure lot C TBAs {see Fig. 7a)} was used in 
the doping study. 

3.2.1 Bulk doping 

Our experiments compared several group-V 
sources in respect of the dependence of the 
donor concentration on Si mole fraction, growth 
temperature, V/III ratio, growth rate, total gas 
flow velocity and substrate misorientation®!). 
We used an inverted horizontal reactor to 
eliminate apparatus issues related to non- 
uniform distribution?’*®. Detailed discussion 
about the inverted horizontal reactor will be 
addressed in the section on AlGaAs. The growth 
rate and donor concentration varied within 
+1-2 % across a 3-inch substrate. The reproduci- 
bility of the measured donor concentration 
value from run to run was within +1 %. 

Figure 8 shows the dependence of the donor 
concentration (Ng) on the disilane flow rate 
using TBAs or arsine for GaAs at 650 °C. Using 


Disilane 
T,: 650°C 


Donor concentration Ng (em) 


10" 
10-5 10~* 107* 


[Disilane] / [TMGa] 

Fig. 8—Dependence of the donor concentration Ng 
on the disilane flow rate divided by the TMGa 
flow rate using lot C TBAs and arsine. The 
V/Il ratio is 20. The group-V flow rate is 
8 x 10 mol/min?)), 
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TBAs, Ng is always about three times higher 
than when using arsine. As the disilane flow in- 
creases, Ng increases in proportion for both 
group-V sources. A similar difference was ob- 
served when silane was used. We confirmed by 
SIMS that Ng was the same as the Si concentra- 
tion. Hall measurements for samples having the 
same Ng demonstrate that the electron mobility 
and compensation ratio values are the same for 
each group-V source®!), These results indicate 
that the difference in Ng was attributable to the 
Si incorporation itself and that most Si atoms 
reside at Ga sublattice sites. SIMS also con- 
firmed that no other residual impurity such as a 
carbon acceptor is present in the epitaxial 
layers. This proves that TBAs enhances Si 
incorporation. 
The effective distribution coefficient, 


ka = {Nq/2.21 x 1022} / { [Si input] / 
[Ga input]} , 


is estimated at 0.1 for growth with TBAs and 
0.03 with arsine. In the treatment of [Si input] 
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Fig. 9—Donor concentration as a function of the growth 
temperature. The solid lines show Ng for layers 
doped with silane, and the broken lines with 
disilane. The Si mole fraction is 2x 10° mol/min 
for both silane and disilane. The V/III ratio is 20. 
The growth rates are constant at 1 um/h in this 
temperature range for both group-V sources?!) 
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in Equation (1) when using disilane, it is postu- 
lated that two atoms per molecule of disilane 
contribute to Si doping reactions. This is be- 
cause in the case of arsine-grown GaAs at a high 
growth temperature*”?, when the disilane flow 
was set to half the value of the silane flow, the 
doping efficiency with disilane was the same as 
with silane. 

Figure 9 shows Ng as a function of the 
growth temperature. 
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Fig. 10—Dependence of Ng on the V/III ratio at 650 °C 
with a growth rate of 1 um/h. The Si mole 
fraction is 2 x 10-8 mol/min. The results for 
EtAs using disilane are included for com- 
parison?!), 
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Fig. 11—Effect of growth rate: The dependence of Ng 
on the V/III ratio at 650 °C with a growth rate 
of 3.6 um/h. The group-V flow rate and the 
total flow rate are the same as in Fig. 3, while 
only the TMGa flow rate is varied to increase 
the growth rate?!) 


Using TBAs, Ng is always higher than when 
using arsine with either silane or disilane. When 
using arsine, the incorporation efficiency with 
silane increases exponentially with the growth 
temperature, approaching that for TBAs. 

Figure 10a) shows the dependence of Ng on 
the V/IIIl ratio at 650 °C when silane is used as a 
dopant source. With TBAs, Ng is always higher 
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Fig. 12—Difference between Si dopant sources derived 
from Figs.10a) and b). The Ng for a layer 
doped with silane subtracted from one doped 
with disilane is shown as a function of the V/III 
ratio?) 


than with arsine for all V/IIl ratios. Ng decreases 
with either group-V precursor as the V/II ratio 
increases. In contrast, when we use disilane | see 
Fig. 10b)} , Ng decreases with TBAs or EtAs, 
and increases with arsine, with an increasing 
V/Ill ratio. Interestingly enough, the trends for 
EtAs match those using TBAs. 

The growth rate should drastically affect 
doping efficiency?*”3539)42)59 When the 
growth rate was increased in this study, the H, 
carrier gas and dopant source flow rates held con- 
stant while the TMGa flow rate was increased. 
With silane, using either TBAs or arsine, Ng 
decreases with an increasing V/III ratio when the 
growth rate is 3.6 wm/h, as shown in Fig. 1 1a). 
The same phenomena are found when the 
growth rate is 1 wm/h {see Fig. 10a)}. With 
disilane, the Ng for arsine-grown layers decreases 
as the V/II ratio increases to 20, as shown in 
Fig. 11b). When the V/III ratio exceeds 20, the 
Ng increases with the V/III ratio. This behavior 
for the arsine-disilane combination differs from 
that observed when the growth rate is slower 
{see Fig. 10b)}. Using TBAs or EtAs, the 
behavior of Ng as a function of the V/III ratio is 
independent of the growth rate. When disilane 
is used, the increase of growth rate from | wm/h 
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to 3.6 um/h elevates the effective distrubution 
coefficient ky from 0.1 to 0.17 with TBAs and 
from 0.03 to 0.08 with arsine. 

Figure 12 shows the difference between the 

silane results {see Fig. 10a)} and the disilane 
results {see Fig. 10b)}. We found that, at a 
higher V/TII ratio, the Ng in the case of disilane 
is higher than that of silane using either TBAs or 
arsine. This trend is independent of the growth 
rate. 
Figure 13 shows the Ng dependence on the 
V/I ratio for epitaxial layers grown using 
disilane at 700 °C. Using TBAs, Ng decreases as 
the V/III ratio increases. The Ng with arsine is 
independent of the V/III ratio in Fig. 13. We 
surmise that this is the same behavior as that 
observed when the V/III ratio exceeds 100 at 
650°C | see Fig. 10b)} , because the pyrolysis 
rate of the As-species increases at higher temper- 
atures!®). 

Figure 14 shows the dependence of Ng on 
flow velocity. The flow velocity is varied by 
modifying only the hydrogen carrier gas flow 
under atmospheric pressure, while all precursor 
mole flow rates are constant. As the flow 
velocity increases, Ng with TBAs_ increases 
and that with arsine decreases using either silane 
or disilane, as shown in Fig. 14. 
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Fig. 13—Ngq using disilane at 700 °C as a function of the 
V/Il ratio. Except for the growth tempera- 
ture, the growth parameters are the same as in 
Fig. 10b)>)), 
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Figure 15 shows Ng as a function of the 
substrate misorientation towards <110> for 
(100) substrates using disilane. The Si incorpora- 
tion for the epitaxial layer grown on the on-axis 
substrate is more efficient than on the misori- 
ented substrate using both TBAs and arsine. 
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Fig. 14—The solid lines show Ng as a function of the 
flow velocity for layers doped with silane at 
600°C. The broken lines show that with 
disilane. The V/III ratio is 20. The Si mole frac- 
tion is 2 x 10°§ mol/min. The growth rate de- 
creases as the flow velocity increases?!) 
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Fig. 15—Dependence of Ng on substrate misorientation 
from the (100) surface at 650°C. The dopant 
source is disilane only. The V/III ratio is 20. The 
disilane flow rate is 2x 10° mol/min. The 
growth rate is 1 um/h?)), 
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3.2.2 Planar doping 

It is important to know whether this differ- 
ent efficiency is due to the group-V precursor 
alone or to the reactions between group-V and 
group-IIl precursors. We also examined planar 
doping of Si into GaAs using TBAs or arsine to 
eliminate the effect of the group-III precursor, 
and we discuss the role of the group-V precursor 
in the Si doping process. Si planar doping 
is a promising technique for obtaining high 
carrier concentrations for high performance 
FETs®?)®*). Therefore, before using it practi- 
cally, we must clarify the effect of TBAs on 
planar doping?” , 

GaAs layers consisted of a 100 nm undoped 
GaAs layer, a planar-doped atomic layer, and a 
100 nm undoped GaAs layer. The growth rate 
of GaAs was lyum/h using either group-V 
precursor. Planar doping was carried out during 
an interruption in the GaAs growth?”, and was 
conducted by supplying disilane under a group- 
V flow without a Ga-precursor for 2 min to 
10 min. 

Figure 16 shows the dependence of the sheet 
carrier concentration (N,) on the disilane flow 
rate when either TBAs or arsine is used for a 
two-minute planar-doped layer. The group-V 
flow rate is 8x 10% mol/min. Although a 
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Fig. 16—Dependence of sheet carrier concentration N, 
on the disilane flow rate for two-minute planar 
doped layers using lot C TBAs or arsine. The 
group-V flow rate is 8 x 10+ mol/min>2), 
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monotonic increase in N, is observed for disilane 
flows of up to 1.5 x 10°* mol/min, N, seems 
to saturate at higher disilane flow rates for both 
group-V precursors. The N, for TBAs is always 
higher than that for arsine. The higher N, of 
2.2 x 10!% cm? is obtained for TBAs, even 
though, because of the increase in the number of 
L-band electrons®*), Hall measurements could 
not account for the extremely high electron 
carrier concentration. For arsine, a high disilane 
flow rate exceeding 10~° mol/min is required to 
obtain a high N, exceeding 1013 cm~?. From 
the C-V method, the FWHM in the depth profile 
of carrier concentration for TBAs-grown layers 
with an N, of 6 x 10!% cm™~ is only 5 nm. The 
peak wave-length of the PL spectrum for this 
layer is the same as that for undoped GaAs 
layers grown with TBAs, indicating an abrupt 
doping profile and no Si diffusion into the 
GaAs layers. 

As the growth temperature increased from 
650°C to 670 °C, N, increased for TBAs from 
155% 10" env’? to. 2% 10 cnr’; As the 
planar doping time increased from 2 min to 
10 min under the same disilane flow rate, N, for 
TBAs increased from 1.6 x 10!* cm™ to 2.2 x 
1012 cm™?, and that for arsine from 7.5 x 
1012 cm? to 1.6 x 1013 cm. Either the high 
growth temperature or long planar doping time, 
however, degrades the surface morphology for 
both group-V precursors, indicating that it is 
difficult to obtain high sheet carrier concentra- 
tion with good morphology using arsine. The 
optimum growth conditions for heavy planar 
doping are restricted for arsine. Thus, we con- 
clude that TBAs is very useful in heavy planar 
doping of Si into GaAs. 

Figure 17 shows N, as a function of group-V 
flow rate for a disilane flow rate of 2 x 10°? mol/ 
min for a two-minute planar-doped layer. 
As the As-precursor flow rate increased from 
8 x 10 mol/min to 4x 10 mol/min, the N, 
for TBAs decreased from 8 x 10!2 cm? to 
6.9 x 101? cm™, while that for arsine increased 
from 2 x 10** cm™ to 5.9 x 10** emi. This is 
the same trend as observed in the dependence of 
Ng on the V/III ratio in bulk doping as shown in 
Fig. 10b), confirming that the reactions with 
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Fig. 17—N, as a function of the group-V precursor 
flow rate with a disilane flow rate of 
2x 10-7 mol/min for two-minute planar doped 
layers>2), 


group-III species have little effect on the carrier 
concentration dependence. We investigated the 
Si concentration in the epitaxial layer using 
SIMS profile as a function of depth from the 
surface. The Si concentrations measured by 
SIMS correspond exactly to the differences in 
N, shown in Fig. 17. This indicates that the 
difference in N, is only due to Si incorporation 
itself and that most of the Si atoms reside at 
Ga sublattices. Since only the group-V precursor 
and disilane were supplied during the planar 
doping sequence, the reactions between the 
group-V species and the Si-species should 
affect the Si doping mechanism, regardless of 
the group-II] precursor. This could be influ- 
enced by the pyrolysis rate of the group-V 
precursor. 

The dependence of N, on the gas flow veloc- 
ity for a two-minute planar doped layer is 
shown in Fig. 18. The gas flow velocity greatly 
affects N,. This behavior is very similar to that 
for bulk doping. The Si concentrations for high 
flow velocities are relatively high when using 
TBAs for both planar and bulk doping. Since the 
planar doping was carried out under atmospheric 
pressure, there is no effect of pressure. 

The flow velocity determines the duration that 
the precursors remain adjacent to the surface. 
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Fig. 18—Solid lines show N, as a function of gas flow 
velocity for two-minute planar doped layers. 
Broken lines show carrier concentration as a 
function of flow velocity for bulk doped layers. 
The disilane and group-V flow rates are 
2x10? mol/min and 8x 10% mol/min for 
both planar doping and bulk doping. The 
growth rate for bulk-doped layers is 1 pm/h>2), 


The duration affects gas phase reactions between 
precursors. As flow velocity increases, the con- 
centrations of As-species and Si-species in the gas 
phase decrease, reducing gas phase reactions. 
Therefore, flow velocity should significantly 
control gas phase reactions. The opposite trend 
in behavior between the group-V precursors in 
Fig. 18 verifies that the progress of the gas 
phase reactions between As-species and Si- 
species will also differ greatly between group- 
V precursors. 

3.2.3, Doping model 
1) Conventional models of Si-doping 

To predict the Si-doping of GaAs, a simple 
equilibrium model represents all Si dopant 
reaching the surface to be incorporated into the 
epitaxial layer due to the low vapor pressure of 
the Si dopant*!?; this is called the mass-transport 
limited case®®. Si incorporation efficiency is 
predicted to be independent of the growth tem- 
perature and group-V flow rate in this model®®. 
From experimental results, however, incomplete 
decomposition of silane causes a strong temper- 
ature dependence?” ”38):40),43),46).47) String- 


fellow has already pointed out that the simple 
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equilibrium model is too simple and inade- 
quate®®). 

Many researchers have reported experiments 
on silane and disilane decomposition®”)*?%). 
Roenigk et al.°5) calculated the silane decom- 
position using the Rice, Ramsperger, Kassel and 
Marcus (RRKM) theory. Field et al.*”’, De 
Croon et al.*®) and Moffat et al.*® reported 
numerical modeling for the growth of Si-doped 
GaAs using two reversible reactions A and B, 
as listed below: 


SiH, @ SiH, + SiH,, san th) 
SiH, © SiH, +H). 


They obtained good agreement with experi- 
mental results. Reactions A and B describe the 
doping process by SiH, (silylene) insertion into 
epitaxial layers*3)»45),46),48),50) These reactions 
can be applied to model Si epitaxy using silane 
and disilane?®)-?®) 

2) Arsine as the group-V source 

First, we discuss the silane doping into GaAs 
using arsine as a group-V source. Using silane, 
we observed that Si incorporation shows a 
strong temperature dependence (see Fig. 9). 
The  temperature-dependent behavior for 
arsine-grown layers is the same as previously 
reported? 7»38),49).43),46),47) | The pyrolysis of 
silane is incomplete for growth temperatures 
between 600°C and 650°C%*®. Thus, -the 
increased incorporation efficiency for arsine- 
grown layers reflects the more effective de- 
composition of silane at higher growth tem- 
peratures®®. The apparent activation energy for 
the arsine-silane process evaluated from the 
slope in Fig.9 is 1.5 eV, similar to previously 
reported values??)?®4% In our reactor, Si 
incorporation is inversely proportional to the 
V/Ml ratio at constant group-IIl partial pres- 
sure | see Figs. 10a) and Ila)}. Bass?” and 
Sakaguchi et al.*”), using atmospheric pressure 
reactors, reported the same behavior as in our 
results. On the contrary, Kuech et al.?®) and 
Moffat et al.*®, using reduced pressure reactors, 
observed that V/II ratio had only a weak effect 
or no effect on Si incorporation. This difference 
could be attributed to the growth pressure*?), 

Using disilane, we observed that, below 


74 


650 °C, the Si incorporation is independent of 
the growth temperature (see Fig. 9). But above 
650 °C, Siincorporation efficiency increases with 
increasing growth temperature as when using 
silane. Furuhata et al.**) reported the same 
trends as ours. However, Shimazu et al.*3) and 
Moffat et al.*° reported different trends. Field 
et al.4”? and Shimazu et al.** observed that the 
dependence of Si incorporation on the growth 
temperature strongly depends on the growth 
pressure. The growth pressure affects the Si 
incoporation activation energy*?). Under atmos- 
pheric pressure, the reactants have a long resi- 
dence time and the gas phase reactions probably 
play an important role. In our reactor, we specu- 
late that below 650°C, reaction A proceeds 
rapidly in the gas phase and the diffusion of SiH, 
radicals controls Si incorporation efficiency. 
Reaction B proceeds at a temperature above 
650 °C, increasing Si incorporation efficiency*”?. 

For the V/II ratio dependence using disil- 
ane, we found that, contrary to when using 
silane, Si incorporation efficiency increases with 
an increasing V/II ratio | see Figs. 10b) and 
11b)} or group-V partial pressure (see Fig. 17). 
This trend is opposite to the silane results 
{see Figs 10a) and 11a)}. Furuhata et al.*®? 
also reported the same trends with disilane as in 
our results. Figure 12 shows the difference 
between silane and disilane. With an increasing 
V/Ill ratio, the difference increases and tends to 
saturate. Thus, we speculate that there is an 
enhancement effect with increased V/II ratios 
when using arsine with disilane. Moffat et al.*® 
introduced reaction C when using arsine as the 
group-V source. 


AsH, + SiH, > H, AsSiH3. lO 


They consider H, AsSiH, (silylarsine) to be 
an additional doping component which can 
explain the V/II ratio dependence. 

3) TBAs as the group-V source with silane 

When using TBAs as the group-V source, Si 
incorporation efficiency of GaAs is always 
higher than that using arsine with either silane 
or disilane under all growth conditions. Serreze 
et al.?®) and Stolz et al.?®) reported high doping 
efficiency using TBAs with silane. Lum et al.°° 
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also reported both high Si doping efficiency 
with (CH; ),Si, and high sulfur doping efficien- 
cy with H,S for TBAs-grown GaAs layers. 
Below 650 °C, the Si incorporation using TBAs 
is independent of the growth temperature, and 
remarkably enough at 600 °C, Ng using TBAs is 


almost ten times that using arsine. In this 
temperature range, no previous research has 


observed such mass-transport-limited behavior 
using silane with arsine. This behavior with 
TBAs indicates that components contributing 
effectively to the doping chemical reactions are 
being formed in the gas phase, in spite of the 
low pyrolysis rate of silane, and that their diffu- 
sion to the surface controls Si incorporation 
efficiency. At a higher growth temperature, 
silane decomposes due to thermal energy, result- 
ing in an increased Ng at 700°C using TBAs 
(see Fig. 9). 

Si incorporation decreases with an increasing 
V/Ill ratio as observed for the TBAs-silane 
combination | see Figs. 10a) and 11a)} . This 
trend is independent of the growth rate. Serreze 
et al.?®) also found the same trend as in our 
results with TBAs, but Stolz et al.?% reported 
the opposite trend. We conclude that there is a 
suppression effect with increased V/II ratio 
using silane as a dopant source with either 
group-V source under our conditions. 

4) TBAs and EtAs as the group-V source with 
disilane 

Using disilane, the dependence on the 
growth temperature is similar for each group-V 
source, but with different incorporation effi- 
ciencies (see Fig. 9). Si incorporation decreases 
with an increasing V/III ratio (group-V partial 
pressure) using TBAs or EtAs {see Figs. 10b) 
and 17){. This trend is opposite to the arsine- 
disilane combination. Figure 12 shows, however, 
that the difference between silane and disilane 
using TBAs has the same trend as in the case of 
arsine. With an increasing V/III ratio, Si incorpo- 
ration efficiency with disilane becomes higher 
than that with silane, using either TBAs or arsine. 
Thus, we conclude that, with disilane, enhance- 
ment occurs in addition to the suppression 
observed with silane. 

Figures 14 and 18 show an opposite depend- 
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ence of Si incorporation on the total gas flow 
velocity between TBAs and arsine using either 
silane or disilane. In this study, flow velocity 
was varied by controlling the hydrogen carrier 
flow rate. The flow velocity affects the growth 
rate and boundary layer thickness. However, 
these dependences were the same for both group- 
V sources?®). The flow velocity determines the 
residence time of precursors adjacent to the 
surface, affecting the reaction time in the gas 
phase. Thus, according to Fig. 14, we conclude 
that rates of gas phase reactions between As- 
species and Si-species should differ greatly be- 
tween the group-V sources. 

Figure 15 shows that Si incorporation effi- 
ciency is high for on-axis substrates. Tang et 
al.5% reported the same dependence on misori- 
entation as in Fig. 15 in the case of the arsine- 
silane combination. Smooth epitaxial growth re- 
quires atomic steps and kink sites at the growth 
fronts and few steps contributing to the growth 
process could exist on the on-axis substrate 
during growth, increasing the Ga-vacancy*°”>??, 
Consequently, the probability of Si adsorption at 
the Ga sublattice increases, leading to increased 
Si incorporation efficiency®®?. The relative 
increment of the on-axis substrate sample 
compared to the misoriented sample is the 
same for both group-V sources, supporting our 


interpretation. 
To verify the influence of self-compensation, 


we investigated SA center peaks in PL spectra, 
which are identified as a mixture of Ga-vacancies 
and Si-donors at the Ga-site”). An SA peak 
appears at a wavelength of around of | um in 
PL spectra for Si-doped GaAs layers. We ex- 
amined disilane-doped GaAs layers with an 
Ng of 5x 10!7cm™3 grown using TBAs or 
arsine. We observed a broad peak at | wm in 
77K PL for both layers. The peak wavelength 
and integrated intensity are the same for both 
group-V sources. Since no difference arises 
between group-V sources and the evaluated 
layers have the same Ng, the probability of 
Ga-vacancy formation during growth is the same 
for both group-V sources. We conclude that the 
difference in Ng is not attributable to a differ- 
ence in self-compensation, but to a difference in 
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the amount of active Si-species on the growth 
surface which could contribute to the doping 
process. 
5) Decomposition of group-V sources 

The strong dependence on the V/III ratio 
and flow velocity supports the concept that the 
simple equilibrium model cannot be adapted to 
Si doping with either silane or disilane, and that 
non-equilibrium conditions control doping fea- 
tures. We speculate that differences in doping 
efficiency among the group-V sources are due 
to the difference in gas phase reactions between 
As-species and Si-species. Even though group-Ill 
sources are supplied during growth, and the 
effect of reactions associated with group-Ill 
sources should be taken into account in dis- 
cussing the doping mechanism, the most signifi- 
cant difference in the growth process with the 
different group-V sources is the pyrolysis of the 
group-V sources, not of group-IIl sources! 7-21). 
We simulated the group-V source pyrolysis 
kinetically, to illuminate the relationship be- 
tween the doping efficiency and the product of 


As-species, and to represent realistic conditions. 

Recently, MOVPE models have been ex- 
panded to fully three-dimensional descriptions 
of horizontal reactors due to the rapid develop- 
ments in computation speed and memory. Model 
predictions have shown good agreement with 
experimental data, such as growth rate, impurity 
levels and composition variation across a wafer 
aligned with the gas flow in horizontal reactors. 
Tirtowidjojo et al.?°”!°° performed the equilib- 
rium calculations for the growth of GaAs from 
TMGa and arsine. They considered the effects of 
temperature, flow rate and inlet gas composition 
at the different reactor pressures. The theoretical 
calculations showed that reaction pathways and 
rate-controlling steps are sensitive to reactor 
pressure. Buchan et al.'°!) calculated the uni- 
molecular rate constants using the RRKM 
theory. Giling et al.!°?) focused on the maxi- 
mum coverage values for the growth species 
adsorbed on the surface under near-equilibrium 
conditions. Coltrin et al.'°?) treated rotating-disk 
reactors based on ideal flow geometry. This 
simulation revealed the relative roles of the 
different reactions, but could not describe the 
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spatial variations in films. Jensen et al.!°) and 
Fotiadis et al.!°°!°® demonstrated fully three- 
dimensional GaAs growth predictions in both 
horizontal reactors and vertical reactors. They 
concluded that the gas phase and surface chem- 
istry, as well as transport descriptions, must be 
included in models to obtain accurate predic- 
tions. 

Our simulation addressed the As-species 
concentration adjacent to the surface as a func- 
tion of reaction time and input group-V flow 
rate, including gas phase and surface reactions, 
and diffusion of the species in the gas phase 
perpendicular to the surface, as reported pre- 
viously !07)»!°8) | The inverted horizontal reactor 
was used in the simulation. We assumed that the 
gas flow in the reactor is well developed and 
unchanging. Since the mole fractions of pre- 
cursors in the reactor are small, the influence of 
energy release on gas flow dynamics due to 
growth is negligible. The change in the concentra- 
tion of species in 20 mesh elements in the 
boundary layer is calculated with gas phase reac- 
tions and the diffusions of species between 
meshes. The height of each mesh element is 
fixed at | mm and the thickness of the bounda- 
ry layer is 20mm. It is postulated that the 
temperature of the inverted substrate and the 
floor of the reactor are 650 °C and 300 °C. The 
gas flow velocity and gas temperature in each 
mesh element increase and decrease linearly 
from the substrate to the floor. The simple 
finite differential method was used to solve 
the differential equations. Calculations of 
the 50 us gas phase and surface reactions are 
iterated ten times followed by 500 us diffusion. 
The initial concentration of precursors is the 
same in all mesh elements. The rate constant k, 
for kinetic simulation, is written as: 


k=AT® exp (—E,/RT). on) 


The reactions considered are given in Table 1. 
Reactions T-1 to T-6 are taken into account 
only for TBAs pyrolysis. We assume that the 
first step of TBAs pyrolysis is a first order 
unimolecular reaction, dissociating into tertiary- 
butyl radicals and AsH, radicals!°)»®?) | and that 
the rate constant of this fission is the same as for 
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the dissociation of C,H,NO into tertiarybutyl 
radicals and NO radicals!®®). Li et al.!*) and 
Speckman et al.''® reported that the first step 
of EtAs pyrolysis is the dissociation into an 
ethyl radical and an AsH, radical. On the 
contrary, Mountziaris et al.!!! reported the 
kinitic simulation of TBAs pyrolysis and postu- 
lated that productionof both AsH and AsH, was 
the first step of TBAs pyrolysis. They treated 


the gas flow, temperature distribution, and sur- 
face reactions very carefully, but considered 
only two gas phase reactions corresponding to 
TBAs pyrolysis, and ignored the effect of C,H,. 
Our simulation treats the gas phase reactions 
more precisely, but other factors are simplified. 
Li et al.” reported that the reaction between 
the TBAs and C,H, produces ditertiarybutyl- 
arsine. However, we could not introduce this 


Table 1. Gas phase reactions for As-species considered in this study 


_ Gas phase reactions for As-species eae ae s) B pecan 
T-1: TBAs > AsH, + C4Hg 16 0 52 
T-2: AsH, °F Gi Hy == AsH3 — C4 Hg 12 0 0 
T-3: C4Ho9 + Cy Hy = Cg Hi, 12 0 0 
T-4: C4 Hy + C4Ho9 = C4Hg + Cq4Hyio 12 0 0 
T-5: C4Ho9 + H, = C4H;o + H rit 0 10 
T-6: C4H9 => C3Hg + CH; 16 0 46.1 
A-1: AsH3 + M > AsH, + H + M LTS 0 75 
A-2: AsH, + M =e AsH + H + M LTA 0 69 
A-3: AH, #8 > Ae + BH, 13 0 2 
A-4: AsH, + AsH, > Assit H 13 0 36.5 
AS: AsH, + AsH, + Ash, + Asi 12 0 10 
A6: As) Hy + Ast, + Asi, 15 0 34.6 
A-7: AsH, + H, + i +H 11 0 27.5 
A-8: AsH3 + CH; = AsH +: (CH. 11.8 0 1.65 
A-9: AsH, + CH; = AsH_ + CHg A 0 0 
A-10: AsH, + AsH ‘AeA 13 0 2 
A-11: AsH + Hy, sg AsH3 11 0 0 
A-12: H + +M = H, + M iy —0.6 0 


The units of A depend on the reaction order, but are given in terms of mols, cubic centimeters, Kelvin, and seconds. 
E, is in kcal/mol. M means third-body species that are buffers of thermal energy?!), 


Table 2. Gas phase reactions for Ga-species considered in this study 


al Gas phase reactions for Ga-species eat ok s) B .. 
Giz  (CH,),Ga > (CH3),Ga + CH; 17.5 0 59.5 
G-2: (CH3 )2 Ga <> CH3 Ga + CH; 7.94 0 35.4 
G-3: CH3 Ga = Ga + CH; 16 0 717.6 
G-4: (CH )3 Ga i > (CH3 )2GaCH,+ Hy, 13.3 0 5 
G-5:  (CH3),GaCH, + Hy +  (CH3);Ga +H 11 0 6.6 
G6: (CH3;),Ga + H + CH;GaCH, + H, 13.1 0 5 
G-7: CH;GaCH, + H, +  (CH3),Ga +H 11 0 6.6 
G-8: CH; Ga + H + CH, Ga + H, 12.8 0 5 
G-9: CH, Ga + H, > CH; Ga + H 11 0 6.6 
G-10: (CH; ),GaCH, = CH3;GaCH, + CH; 17.5 0 59.5 
G-12: CH3GaCH JEL = CH, Ga + CH; 175 0 40 
Gll: GH; + CH; 3 C2H¢ 13.9 —0.4 0) 
G-13: CH; + H + M > CH + M 36.9 3 0 
G-14: CH; + H, > CH, + H 11.9 0 10.9 
G-15: CH, + H > CH; + Hy, 4.34 3 8.75 


The simulation for As-species pyrolysis calculates Ga-species pyrolysis at the same time3). 
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reaction into the simulation due to the absence 
of knetic data for it. Since the rate of AsH, 
formation is predicted to be fast for TBAs, 
many gas phase reactions must be considered, as 
listed in Table 1. The unimolecular reaction for 
arsine is calculated using the RRKM theory!°)). 
Reactions A-8 and A-9 with the CH; radical are 
introduced to consider the effect of group-Ill 
sources?”?. To estimate the concentration of 
CH, radicals, TMGa pyrolysis is also simulated 
considering the reactions in Table 2. We assume 
that no adduct formation between Ga- and As- 
species contribute to the growth reactions in 
the inverted horizontal reactor, as mentioned in 
the following AlGaAs section!®*®”). The surface 


Table 3. Surface reactions for As-species considered 
in this study 


Surface reactions for As-species 


AsH3 a AsH) +H 
AsH, > AsH +H 
AsH3 +H =F AsH), +H, 
AsH +H —F AsH3 

AsH) +H = AsH + Hy 
AsH > As (growth) 


Interactions between As-species are not involved?!) ; 
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reaction is considered only for the As-species 
with hydrogen (see Table 3). The adsorption 
and desorption steps are also considered using 
rate equations. 

Figure 19 shows the results of our simula- 
tion; in it we plot the time dependencies of the 
chemical species concentration in the gas phase 
adjacent to the surface. We calculated that As- 
species radicals are abundant when using TBAs, 
suggesting to us that As-species radicals are 
responsible for the enhanced Si incorporation. 
6) Doping mechanism with TBAs and EtAs 

We propose the following scheme: 


AsH, + SiH, > H, AsSiH; +H, ew GD) 
AsH_ + SiH, > H, AsSiH3. 
In our proposed mechanism, reactions D and 
E control Si doping efficiency rather than reac- 
tions A and B. The following reactions are intro- 


duced to explain the difference between group- 
V sources from our simulation: 


TBAs > AsH, +C,H,, (T-1) 
AsH, + C4Hy > AsH,;+CyH,, (T-2) 
AsH, +M-—> AsH, +H+M. (A-1) 


0 100 200 300 
Time (ms) 


b) Arsine 


Fig. 19—Results of simulation: As-species concentration in the gas phase adjacent to the surface as a function of 


the reaction time for GaAs growth. 


The initial input condition represents a V/I ratio of 20 and a growth temperature of 650 “eh. 
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Our simulation indicates that large amounts 
of As-species radicals exist adjacent to the sur- 
face, resulting in rapid rates of reaction for reac- 
tions D and E. These reactions produce a large 
amount of silylarsine. Moffat et al.*® has modi- 
fied a doping model associated with an addition- 
al effect of reaction C. In contrast, we speculate 
that silylarsine formed via reactions D and E 
mainly controls doping efficiency, resulting in 
a higher donor concentration using TBAs (see 
Figs. 8 and 16), because silylarsine is stable in 
the gas phase due to its low standard heat of 
formation of 37 kcal/mol and its high disso- 
ciation energy of 73 kcal/mol!!”), whereas 
the standard heat of formation for AsH, is 
41 kcal/mol!» and for SiH, is 63.5 kcal/mol?®), 

As-species radicals can adsorb at the Ga- 
sublattice during As-rich growth, resulting in 
an As-stabilized surface!!”. Si-precursors must 
decompose and then break the As-dimerized 
bond to incorporate into the epitaxial layers. 
This difficulty causes a low distribution coeffi- 
cient, typically 0.03 when arsine is used. When 
silylarsine is formed in the gas phase and reaches 
the growth surface, dissociation should occur 
immediately, producing activated As- and Si- 
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Lo? 
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10-” 


150 ms 


0 100 200 


V/Ill ratio 
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species radicals at the same place. Moffat et 
al.*®) assumed that the surface reaction proba- 
bility of silylarsine is 1. After the dissociation 
reaction on the surface, we speculate that most 
As-species radicals incorporate into an As-site 
and that Si-species should adsorb on the Ga- 
sublattice sites without difficulty, resulting in 
a higher distribution coefficient. However, 
silylarsine can also adsorb effectively near Ga- 
sites due to the nature of the Si-atom, producing 
ASc, in the epitaxial layer. We observed a large 
amount of Asc, for TBAs-grown GaAs layers 
using the DLTS (deep level transient spectros- 
copy) method, suggesting the existence of 
silylarsine on the surface!°®). Even in this case, 
we speculate that few Si atoms can incorporate 
into the As-sublattice sites due the intrinsic 
property of GaAs. Thus, we concluded that the 
silylarsine, predicted to be abundant with TBAs, 
should enhance Si incorporation into the Ga- 
sublattice sites. The rates of reactions D and E 
are expected to be fast at 600 °C because reac- 
tion T-1 occurs at quite a lower temperature. We 
surmise that the doping process between 600 °C 
and 650 °C is in the silylarsine transport limited 
region, consistent with the observed tempera- 


0 100 200 
V/Ill ratio 
b) Arsine 


Fig. 20—Results of simulation: Dependence of As- species concentration in the gas phase adjacent 


to the surface for the V/III ratio. 


The reaction time is fixed at 150 ms. The growth tempera- ture is 650 at A 
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ture-independent trends for silane doped TBAs- 
grown layers (see Fig. 9). 

We also propose the surface reaction mech- 
anisms related to the AsH; molecule. Most silane 
decomposes on the surface via adsorption®” 
similar to arsine. Tang et al.5” also reported 
that the decomposition of silane happens very 
close to the hot surface. Silane should compete 
with AsH3 in adsorbing®*” °*), Excess AsH3 
adsorption inhibits Si-species adsorption to 
break the As-dimerized bond on the As-stabi- 
lized surface and to incorporate into the epitax- 
ial layer. We speculate that this mechanism is 
responsible for the depletion effect with the V/II 
ratio when using silane {see Figs. 10a) and 1 la)}. 
Figure 20 shows the As-species concentration 
adjacent to the surface after 150 ms as a func- 
tion of the V/III ratio from our simulation. The 
AsH3 concentration adjacent to the surface 
increases with increasing As-precursor flow 
rate with either arsine or TBAs. We speculate 
that reaction T-2 is predominantly responsible 
for the production of a large amount of AsH3 
in TBAs pyrolysis. 

Using disilane produces more complicated 
situations. Disilane rapidly dissociates into SiH) 
and SiH, via reaction A. Thus, due to a large 
amount of SiH», Si incorporation when using 
arsine is controlled by the diffusion of SiH» (see 
Fig.9). We speculate that, even using disilane, 
reactions D and E proceed effectively with TBAs, 
resulting in enhancement {see Figs. 10b) and 
17}. Moreover, we have to consider reaction C 
to be significant for the case of disilane. In 
reaction C, an SiH, radical reacts with AsH3, 
producing silylarsine*®?. AsH3 increases with an 
increasing V/III ratio (see Fig. 18). Therefore, 
using only disilane, as the V/III ratio increases, 
the rate of reaction C increases, producing a 
large amount of silylarsine and enhancing Si 
incorporation not only using TBAs but also 
using arsine (see Fig. 12). 

The only difference in growth conditions 
between Figs. 10b) and 11b) is the TMGa flow 
rate. The group-V flow rate and disilane flow 
rate are the same. According to our proposed 
mechanism, the concentrations of both silylarsine 
and silylene are the same. However, the effective 
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distribution coefficient kg in Fig. 11b), higher 
growth rate case, is higher than that in Fig. 10b). 
Sakaguchi et al.*”) reported similar results. On 
the contrary, Tang et al.°° observed the 
opposite trend. In our experiment, the increase 
in growth rate corresponds to a decrease in the 
V/IIl ratio, increasing the probability of Ga- 
vacancy formation during growth. We speculate 
that the increase of Ga-vacancies is responsible 
for the high doping efficiency under a constant 
doping component concentration. 

Our simulation indicates that the rate of 
formation of the As-species in the gas phase 
differs between the group-V sources by about 
150 ms (see Fig. 19). Precursors pass beyond 
the hot susceptor of 46mm before reaching 
measurement point: in this study 153 ms under 
atmospheric pressure with a flow velocity of 
30 cm/s. As the reaction proceeds, AsH, pro- 
duction with TBAs increases, while that with 
arsine decreases (see Figs. 14 and 18). A higher 
flow velocity corresponds to a shorter residence 
time of reactants in the gas phase; a slower flow 
to a longer residence time. Thus, comparing 
Fig. 19 with Figs. 14 and 18, we conclude that 
different dependencies on the flow velocity 
(or on the reaction time) are due to differences 
in the rates of As-species formation between 
the group-V sources. 

We speculate that the same discussion as for 
TBAs should be used for the EtAs-disilane com- 
bination because of the similarity in the decom- 
position reactions. Li et al.'*) and Speckman et 
al.'!°) have noted that the dominant reaction in 
the decomposition of EtAs is: 


EtAs > C)Hs5 + AsH>. 


This reaction occurs at a lower temperature 
than for the pyrolysis of TBAs!!® 

3.2.4 Deep level concentration 

Two lots of TBAs, i.e. lot A (contaminated 
by sulfur) and lot B (as pure as the recent batch) 
as mentioned in section 3.1, were used in the 
DLTS”*? study. Disilane was used as a dopant 
during growth using lot B. Doping was inten- 
tionally not done during growth using lot A. The 
V/Ill ratio was 10 for all samples. The growth 
temperature was between 575 °C and 700 °C. 
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The growth rate was almost constant in this 
temperature range as shown in Fig.2. The 
EL2 concentration was thus not affected by the 
growth rate in this study. 

In a typical DLTS spectrum of GaAs with 
TBAs, one negative peak was detected. The 
activation energy is 0.83 eV and the capture 
cross-section is 6.5 x 107-3 cm?. The Arrhenius 
plot for this trap falls near those for EL2!!*), 
ETI''®), and EB2!!”. These are identified as 
the same trap''®), the so-called EL2. EL2 is a 
dominant electron trap in GaAs, and we believe 
that the trap observed here is EL2. EL2 was 
also detected in GaAs with arsine. 

Figure 21 shows the growth temperature 
dependence of the EL2 concentration in GaAs 
with TBAs or arsine. With TBAs, the EL2 
concentration decreases with an_ increasing 
growth temperature below 675°C, and _ is 
constant above 675 °C. This does not depend on 
the type of donor impurity. The EL2 con- 
centration is constant for GaAs with arsine in 
this growth temperature range. 

Much work and time have gone into clarify- 
ing the origin of EL2, a native defect in GaAs. 
It has been attributed to substitutional As, 
(Asga)'™, or a complex of Asg, and inter- 
stitial As, (As;)'?%. EL2 concentration is 
affected by excess incorporation of As into the 
crystal'?!) and, therefore, depends on the type 
of reactant on the substrate. We speculate that 
no EL2 is introduced when the As-As bond 
dissociates on the surface. When using TBAs, 
AsH) radicals dominate the gas phase as men- 
tioned in subsection 3.2.3. Thus, at tempera- 
tures below 650 °C, excess As is supplied on the 
surface with TBAs, resulting in high EL2 con- 
centration (see Fig. 21). EL2 is completely absent 
in GaAs grown by molecular-beam epitaxy 
(MBE), where the dominant arsenic source is 
As, **??, Taking the origin of EL2 into con- 
sideration, this fact shows the possibility that 
excess incorporation of As decreases when As) 
molecules are supplied on the surface of the 
epitaxial layer. Thus, at temperatures above 
650°C, As) molecules are formed by AsH) 


consumption with TBAs°??, resulting in low 
EL2 concentration (see Fig. 21). 
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Fig. 21—Growth temperature dependence of the EL2 


concentration®”), 


The 50 percent decomposition temperature 
of arsine is 675 °C'*. Since this is much higher 
than that of TBAs, the gas phase reaction does 
not seem to contribute to the As formation 
before reaching the substrate, resulting in the 
constant EL2 concentration between 575 °C 
and 700 °C (see Fig. 21). 


3.3 Undoped AlGaAs 

Growth of ternary or quaternary alloys is 
required to fabricate many kinds of device. Not 
only purity of epitaxial layers but also uniformi- 
ty of solid composition, growth rate, and carrier 
concentration are key features in this growth. 
In this section, we discuss AlGaAs growth on 
GaAs. AlGaAs has proven to be a more difficult 
material to grow by MOVPE than GaAs due to 
the high reactivity of aluminum, forming strong 
bonds with carbon and oxygen and degrading 
uniformity. 

3.3.1 Conventional horizontal reactor 

The surface morphology of the AlGaAs 
epitaxial layers was excellent under all growth 
conditions?>?. Figures 22a) and b) show the uni- 
formity of the growth rate and the AlAs mole 
fraction (x value of Al, Ga,_, As epitaxial layer) 
along a wafer grown using TBAs and arsine at 
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Fig. 22—Uniformity of AlGaAs along a wafer parallel to 
gas flow for Alo.23Gao.72As epitaxial layers 
grown at 680 °C with a V/II ratio of 60. AlAs 
mole fraction in vapor is fixed at 0.24. The uni- 
formity of films using arsine is independent of 
the total gas flow rate?>), 


680 °C. The growth was less uniform for TBAs 
than for arsine. We found that y, defined as the 
xX value in the solid form divided by the x value 
in the vapor form, varied from 1.1] to 1.5 across 
the wafer for TBAs and was constant at 1.1 for 
arsine under the same growth conditions. 
Increasing the flow rate, decreasing the growth 
temperature, and cooling the ceiling of the 
reactor improved y uniformity across the wafer. 
The growth uniformity for films using arsine 
does not depend significantly on the total gas 
flow rate. Figures 23a) and b) show the 
dependence of the growth rate and x value on 
growth temperature at the center of the wafer. 
The growth rate decreases and the x value 
increases with increasing temperature using 
TBAs with a fixed x value of 0.24 in vapor. On 


the other hand, for arsine, the growth rate 
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Fig. 23—Dependence of AlGaAs growth on growth tem- 
perature for Alg.2gGao.72As epitaxial layers 
grown with a V/III ratio of 60. AlAs mole frac- 
tion in vapor is fixed at 0.24. Only data for the 
center of the wafer is given in this figure?>) 
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g. 24—Calculated AlAs growth rate and calculated 
GaAs growth rate profile along a wafer for films 
grown at 680°C. AlAs and GaAs growth rates 
are determined by multiplying the Al,Ga,_,As 
growth rate and the x value or 1-.x value 
shown in Fig. 337°). 
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decreases slowly with increasing temperature 
and the x value is quite independent of tempera- 
ture. 

The calculated growth rate profile of AlAs 
along the wafer parallel to the gas flow for TBAs 
is the same as that for arsine, as shown in Fig. 24. 
Here, the AlAs growth rate was determined from 
the product of the AlGaAs growth rate and the 
AlAs mole fraction shown in Fig.22. The 
growth rate of GaAs using TBAs decreases more 
rapidly in the direction of flow than that using 
arsine. A higher y and a less uniform growth rate 
for TBAs compared with arsine are caused by 
depletion of Ga-sources in vapor. Black deposits 
on the ceiling of the horizontal reactor were 
notable. A gas phase reaction resembling adduct 
formation should occur as mentioned above. In 
addition, TMAI is unlikely to react with TBAs to 
form an adduct in vapor, since TMAI is a 
dimmer in vapor!??). The gas phase reactions 
between TMGa and TBAs promoted Ga-group 
deposits, i.e. (GaCH3 \n’, to form on the 
ceiling of the reactor, causing the GaAs growth 
rate to decrease along the wafer. Thus, we have 
to reduce gas phase reactions to control growth 
strictly. Possible options using a conventional 
horizontal reactor are to increase the flow 
velocity or to reduce the reactor pressure. 

The FWHM of the X-ray diffractometer 
rocking curves for AlGaAs films using TBAs at 
700°C with a V/III ratio of 60 was 27 arcs. 
This value is equivalent to that of 30.4 arcs for 
arsine. This suggests that films grown with TBAs 
are good crystals. Undoped AlGaAs epitaxial 
layers grown at 670°C to 700°C with an x 
value of 0.24 in the vapor phase showed n-type 
conductivity with a carrier concentration of 
10'> cm to 10'® cm, and high resistivity 
at 650 °C. Si was detected as an n-type dopant 
with a concentration of 10'®° cm by SIMS 
analysis. The Si concentration is close to the 
carrier concentration. Thus, we conclude that Si 
is the main n-type dopant. The behavior of the 
carrier concentration is similar to that using 
arsine. We speculate that Si comes from TMAI 
since the GaAs epitaxial layers grown using 
TBAs and TMGa are of high quality. 

Figure 25 shows the carrier concentration 
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Fig. 25—Carrier concentration and full-width at half- 
maximum (FWHM) of PL band edge emission 
at 77K for Alo.2.3Gao.72 As grown using lot B 
at 680 °C as a function of the V/II ratio23). 


and the FWHM of the PL spectra measured at 
77 K for films grown at 680°C, as a function 
of the V/III ratio. The carrier concentration 
increased slightly with increasing V/III ratio for 
films using not only TBAs but also arsine. The 
FWHM decreased with increasing V/II ratio. 
From the FWHM results, we conclude that a 
larger V/I ratio produces better AlGaAs 
epitaxial layers. 

Figure 26 shows the PL spectra at 4.2 K of 
undoped Alp 2g Gag.72As layers grown at 680 °C 
with V/IlI ratios of 20, 60, and 120. The spectra 
for films using TBAs and arsine are quite similar 
for the same V/III ratios. The carbon-related 
peak is prominent for epitaxial layers grown 
with small V/IIl ratios. We suppose that the 
higher binding energy between aluminum and 
carbon forces the carbon in organometallics to 
incorporate into AlGaAs epitaxial layers, while 
carbon incorporation is negligible in the GaAs 
epitaxial layers. We speculate that hydrogen 
radicals hardly remove hydrocarbons that bind 
with aluminum from the surface. However, 
the ratio of the exciton peak intensity to the 
carbon acceptor peak intensity increases with 
increasing V/III ratio. The exciton emission for 
films grown with a V/III ratio of 120 is quite 


strong. The SIMS analysis showed that the 
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Fig. 26—PL spectra measured at 4.2 K for 
Alo.23Gao.72 As layers grown at 680°C with 
three V/III ratios?>). 


concentration of residual carbon had a magni- 
tude of 10'® cm“? for films with a V/III ratio of 
20. Since undoped AlGaAs films showed n-type 
conductivity with a free carrier concentration of 
10'° cm=3, we suppose that carbon incorpora- 
tion as a p-type dopant is as low as 10'* cm™ 
to 10'® cm73. An increased gas flow rate also 
improved the quality of the AlGaAs. At the 
beginning of this work, we expected that high 
quality AlGaAs layers would be grown with even 
a small V/III ratio when using TBAs, since the 
decomposition temperature of TBAs was quite a 
lot lower than that of arsine. From this work, 
however, it seems that a higher V/II] ratio is im- 
portant for the growth of high quality AlGaAs 
when using TBAs as well as when using arsine. 

3.3.2 Inverted horizontal reactor 

In AlGaAs growth using TBAs with conven- 
tional horizontal atmospheric pressure reactors, 
epitaxial layers showed poor growth uniformity 
as discussed in subsection 3.3.1. Increasing the 
total gas flow rate improves uniformity slightly 
but decreasing the growth temperature has no 
apparent effect. This was mainly due to the 
atmospheric pressure apparatus. Because black 
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deposits on the ceiling of the reactor just above 
the substrate consumed group-III species and 
decreased the source supply over the entire 
substrate surface, the uniformity in growth be- 
came worse. Deposits were more prominent for 
TBAs than for arsine. The natural convection 
induced by buoyancy in all conventional atmos- 
pheric pressure reactors promoted deposition 
onto the ceiling, especially when used TBAs. 

We estimated the influence of natural con- 
vection using the Rayleigh number (Ra), defined 
as follows!?*?: 


Ra = (gBh3 AT)/vDr, 


where g is the acceleration due to gravity, 6 the 
coefficient of thermal expansion, h the typical 
length of flow, AT the typical difference in tem- 
perature, v the kinetic viscosity, and Dy the 
coefficient of thermal diffusion. Natural convec- 
tion strongly affects the flow when Ra exceeds 
1708!75). Ra was evaluated at 311 for our 
horizontal reactor at atmospheric pressure. Con- 
vection does not disorder the flow completely, 
but the influence cannot be ignored for an Ra 
of 311. Thus, we developed an inverted horizon- 
tal atmospheric pressure reactor, as shown in 


Fig. 1b), to eliminate the buoyancy-induced 
flow and reduce parasitic deposition so that 


epitaxial growth affected only by the nature of 
the materials themselves could be carried out 
easily under atmospheric pressure. No parasitic 
deposition was observed in the inverted reactor 
for either source. Puetz et al. obtained uniform 
InGaAs with an inverted horizontal reduced- 
pressure reactor 126) We found an additional 
advantage to using the inverted reactor. Excel- 
lent morphology could be obtained easily for all 
epitaxial layers, because few particles were ad- 
sorbed on the substrate due to face-down 
growth. We obtained epitaxial layers with a total 
particle density of 3 cm~ larger than 0.24 um? 
using TBAs. 

The uniformity of Al,Ga,_,As films was 
greatly improved using the new inverted reactor. 
In the PL spectrum at 77 K, the variation in the 
wavelength of the band edge emission line was 
within +0.5 nm across the substrate, showing 
excellent x uniformity. Figure 27 shows the 
distribution of thickness across a_ stationary 
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Fig. 27—Variation of Alg.23Gao.72 As thickness using lot 
B TBAs and arsine across a 3-inch stationary 
substrate parallel to the gas flow at 650 °C with 
a V/III ratio of 60 in a conventional and in- 
verted reactor?”). 


3-inch substrate at 650 °C with a V/III ratio of 
60. A high V/III ratio is important for AlGaAs 
growth as shown in Fig. 26. With the inverted 
reactor, the uniformity of AlGaAs using TBAs 
and arsine was virtually the same. 


3.4 Si-doped AlGaAs 

An n-type AlGaAs with a carrier concentra- 
tion greater than 10!7 cm” is essential to fabri- 
cate HEMTSs or lasers consisting of AlGaAs. 

An initial study on _ Si-doping for 
Alo.2gGao.72 AS has been completed using 
disilane. Undoped Alop..gGag.72As grown with 
a high V/III ratio using TBAs or arsine possesses 


a low carrier concentration as shown by the Hall 
measurement and low carbon incorporation ac- 


cording to the PL spectra (see Fig. 26). AlGaAs 
grown using EtAs at 650 °C, however, resulted 
in poor morphology due to intermediate reac- 
tions and to residual impurities such as water. 
Figure 28 shows the dependence of the Si in- 
corporation efficiency on the disilane flow rate 
for epitaxial layers grown using TBAs or arsine 
at 650°C. The Ng using TBAs is always 1.5 
times higher than using arsine. This occurs in the 
case of GaAs growth, indicating that the same 
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Fig. 29—Distribution of Ng and layer thickness of Si- 
doped Alg.2g Gag.72 As layers grown using TBAs 
at 650°C with a V/III ratio of 60 across a 
rotated 3-inch substrate. The average Ng is 
1.81 x 108 cm. The average layer thickness 
is 90.4 nm?4), 


mechanism controls Si incorporation for both 
GaAs and AlGaAs. 

AlGaAs layers were grown on substrates 
rotated at 12 r/min. Figure 29 shows the distri- 
bution of Ng and the layer thickness using 
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TBAs. Despite atmospheric pressure growth and 
a vapor phase reaction between TMGa and TBAs 
caused by the low decomposition temperature 
of TBAs!°”?%), variations in Ng and the layer 
thickness were less than +1.2 % and 40.9%. An 
edge effect due to protuberances on the suscept- 
or for holding substrates slightly decreased 
layer thickness at the edge region. This result is 
not optimized but already similar to that using 
arsine, suggesting that a uniformity sufficient for 
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Fig. 30—Distribution of band edge emission in PL 
spectra at 77K for In,Gaj,-)As layers across 
a stationary GaAs substrate using TBAs grown 
at 630 °C. Only TMIn flow rate is varied?%), 
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device fabrication can be obtained using TBAs 
with this inverted reactor even at atmospheric 
pressure. Thus, we believe that improving the 
reactor to eliminate deposits should produce 
excellent AlGaAs uniformity using TBAs. 


3.5 InGaAs on GaAs 

In this section, In,Ga,_,As layers were 
grown on GaAs _ substrates aiming to 
AlGaAs/InGaAs/GaAs_ pseudomorphic hetero- 
structures. All layers were 1.2 wm thick in order 
to minimize the distortion caused by lattice 
mismatch. The V/IIl ratios of InGaAs grown 
were 2.6 to 10 for TBAs and 5 to 95 for arsine. 
The surfaces of InGaAs layers grown using arsine 
with a V/III ratio of 5 were milky due to insuf- 
ficient decomposition of arsine. 

Figures 30 and 31 show the distributions of 
band edge emissions in PL spectra at 77 K (cor- 
responding to InAs mole fraction; v) across 
stationary 3-inch substrates grown at 630 °C and 
650 °C, respectively. In Fig. 30, only TMIn flow 
rate was varied and the V/IIl ratio was around 
five. Excellent y uniformity was obtained for 
the layers with y up to 0.15. In Fig. 31, InAs 
mole fraction in gas phase (),,;) was kept 
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Fig. 31—Distribution of band edge emission in PL spectra at 77 K for In, Ga,-,As layers across a 
stationary GaAs substrate with various V/III ratios grown at 650 ea 
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to 0.2. An increase of only 20 °C deteriorated 
the y uniformity significantly using TBAs. There 
are substantial differences between sources. The 
y uniformity for TBAs became worse with in- 
creasing V/III ratio. However, y was almost con- 
stant in the region upstream from the center and 
the InGaAs thickness was also confirmed to be 
uniform in that region from the growth of 
GaAs/InGaAs/GaAs quantum well structures. 
For arsine, y increases slightly downstream from 
the top, regardless of V/III ratio, up to 95. 

No deposit was left on the reactor ceiling 
just under a substrate nor on the injector for 
either source. Depositions on the cool wall 
downstream for TBAs were more abundant than 
those ‘for arsine. It should be noted that, in 
AlGaAs growth at 650 °C, only reducing source 
consumption due to wall deposition resulted in 
highly uniform AlAs mole fraction and highly 
uniform layer thickness. Thus, InGaAs growth 
distributions across a substrate are affected only 
by the material itself. 

Since there was no deposition on the injec- 
tor, there is no apparent reaction between TMIn 
and arsenic sources at room temperature. There- 
fore, we speculate that different reaction mecha- 
nisms operate in the vapor phase only around 
substrates including TMIn at atmospheric 
pressure. They may be reactions involving 
adduct formation or a parasitic reaction, which 
does not contribute to the growth reaction 
directly like triethylindium (TEIn)-arsine adduct 
formation!?”): either may eliminate In-contain- 
ing species from the vapor just below the sub- 
strate. This is because the temperature of the 
vapor around the substrate is only 200 °C to 
500°C, even if the growth temperature is 
700 °c!°5). TBAs begins to decompose at 
300 °C!9))?1)_ and is expected to react with 
TMIn in vapor at that temperature. This unclari- 
fied vapor phase reaction must be eliminated 
to fabricate pseudomorphic structures. Since 
InGaAs layers are already uniform in the up- 
stream region tor TBAs (about 40 mm length), 
we believe that improved source injection or 
reduced pressure should produce uniform 
InGaAs across a 3-inch substrate. 


FUJITSU Sci. Tech. J., 28, 1, (March 1992) 


3.6 Selectively doped heterostructures 

The quality of selectively doped heterostruc- 
tures is very sensitive to the abruptness at the 
heterointerface, the purity of materials, and 
doping concentration. Thus, we investigated the 
growth of selectively doped heterostructures 
with TBAs, comparing that with arsine. These 
results were utilized for HEMT fabrication as 
described in section 3.7. 

3.6.1 n-AlGaAs/GaAs 

Selectively doped n-AlGaAs/GaAs_hetero- 
structures consist of a 700 nm undoped GaAs 
layer, a Onm to 7 nm undoped Alo..gGao.72 AS 
spacer layer, a 37 nm Si-doped Alp. Gag,72 AS 
layer, and a GaAs cap layer, The electron con- 
centration of the n-AlGaAs layer was maintained 
at 1.4 x 1018 cm™3. The V/III ratios of GaAs and 
AlGaAs for TBAs were 40 and 60 and those for 
arsine were 10 and 60. 

Strong oscillations of the Shubnikoy-de Haas 
(SdH) effect in a transverse magnetoresistance 
were observed at 4.2 K for the heterostructures 
of the 5 nm spacer layer grown using TBAs as 
shown in Fig. 32. The magnetic field was applied 
perpendicular to the heterointerface. These 
oscillations confirm the existence of high-quality 
2DEG at the heterointerface. 


Magnetoresistance (arb. units) 


~ 


0 1 2 3 4 5 6 
Magnetic field (7) 
Fig. 32—Shubnikov-de Haas oscillations observed at 4.2 K 
for a selectively doped n-Alo.2g Gao.72 As/GaAs 
heterostructure grown using TBAs®). 
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Fig.33—Electron mobility and sheet carrier con- 
centration at 77K _ for selectively doped 
n-Alo 28 Gao.72 As/GaAs heterostructures grown 


at 650 °C as a function of AlGaAs spacer layer 
thickness2*). 


Figure 33 shows electron mobility (wu) and 
sheet carrier concentration (N,) at 77K as a 
function of spacer layer thickness for films 
grown using lot C (produced in 1990). The 
inverted horizontal reactor was used. There is no 
difference between the group-V sources. We ob- 
tained high-quality n-AlGaAs/GaAs heterostruc- 
tures using TBAs with a uw of 98000 cm? V!s71 
and an N, of 6.9 x 10!! cm at 77K for the 
7 nm spacer layer. This proves that an abrupt 
interface between the AlGaAs spacer layer and 
GaAs buffer layer was formed using the inverted 
reactor. The excellent performance of the 
heterostructures also verifies the purity of these 
epitaxial layers. 

3.6.2 n-AlGaAs/InGaAs/GaAs 

In the growth of n-AlGaAs/InGaAs/GaAs 
pseudomorphic heterostructures, a 15nm 
InGaAs channel layer was grown between the 
undoped GaAs layer and the undoped AlGaAs 
spacer layer. In the atmospheric pressure re- 
actor, the InAs mole fraction in the InGaAs layer 
using TBAs was less uniform than that in the 
layer using arsine as described in section 3.5. 
A V/IIl ratio of less than 5 was required to 
obtain better uniformity with TBAs. However, 
with arsine, InGaAs grown with such a low V/III 
ratio resulted in a milky surface due to insuffi- 
cient arsenic source supply on the substrate. 
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Fig. 34—Electron mobility and sheet carrier concen- 
tration of n-Alo.23Gao.72 As/In, Ga,-, As/GaAs 
pseudomorphic heterostructures at 300 K as a 
function of the InAs mole fraction in the 
InGaAs channel layer®?), 


Thus, when pseudomorphic HEMTs were grown, 
the InGaAs channel layers had a V/III ratio of 
2.5 for TBAs and 10 for arsine. The growth 
parameters (e.g. growth rate) were optimized at 
the center of the substrate. 

Figure 34 shows yw and WN, of selectively 
doped n-AlGaAs/InGaAs/GaAs pseudomorphic 
heterostructures at room temperature as a func- 
tion of the InAs mole fraction (vy). The AlGaAs 
spacer layer was 2.5 nm thick. As the InAs mole 
fraction increased, N, increased monotonically, 
while the uw remained almost constant when 
the InAs mole fraction was below 0.15. At 
300 K, N, and pw were 1.4.x 10!2 cm? and 
7200 cm?V~!s!_ respectively for the film 
grown using TBAs with an InAs mole fraction of 
0.15. The N, and pw at 77 K for this film were 
1.25 x 10!* cm? and 45000cm?V™"s" re- 
spectively, demonstrating an abrupt interface be- 
tween the AlGaAs spacer layer and the InGaAs 
channel layer. The uw resulting from using TBAs 
was higher than that when arsine was used. On 
the other hand, when TBAs was used, N, was 
lower than that for arsine. In the growth of 
pseudomorphic heterostructures, the V/II ratio 
for the growth of the InGaAs channel layer is 
different between the group-V sources, because 
the V/IIl ratio for optimum uniformity is differ- 
ent, as shown in Fig. 31. The Si incorporation 
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efficiency also differs between sources?®)™3?). 


Thus, we speculate that the difference in electri- 
cal properties is attributed to the difference in 
the quality of both Si-doped AlGaAs and 
InGaAs layers. 


E-HEMT D-HEMT 


| | | | 


Gate metal 


Insulator 


Pa 


n-AlGaAs 


Ohmic metal 


GaAs substrate 
2DEG 


Fig. 35—Circuit structure of m-AlGaAs/GaAs HEMT 
DCFL. 


Ver. : 1mA/div 
Hor.: 200 mV/div 


Fig. 36—Source-drain /-V characteristics at 300 K for 
n-Alo.28 Gao,72 As/GaAs E-HEMTs with a 
0.5-um gate using TBAs. The step of source- 
gate voltage V,, is 0.1 vo). 


Ver.: 1 mA/div 
Hor.: 200 mV/div 


Fig. 37—Source-drain /-V characteristics at 300 K for 
n-Alo.28 Gao.72 As/Ing 15 Gao.g5 As/GaAs 
E-HEMTs with a 0.5-um gate using TBAs. 

The step of source-gate voltage Vg, is 0.1 vy), 
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3.7 HEMT fabrication 

In this chapter, we describe the fabrication 
of HEMTs using the heterostructures grown as 
described in section 3.6. 

3.7.1 n-AlGaAs/GaAs HEMT 

Figure 35 shows an n-AlGaAs/GaAs HEMT 
direct-coupled-FET-logic (DCFL) circuit struc- 
ture. The gate length is 0.5 wm and the width is 
60 wm. A very thin AlGaAs layer is embedded 
in the GaAs cap layer to stop selective dry 
etching during the fabrication of enhancement- 
mode HEMTs (E-HEMTs) and depletion-mode 
HEMTs (D-HEMTs). The source-drain current- 
voltage (/-V) characteristics of n-AlGaAs/GaAs 
HEMTs grown using TBAs are shown in Fig. 36. 
Excellent saturation and pinch-off were ob- 
served using TBAs. We achieved high perform- 
ance HEMTs, with an extrinsic transconduc- 
tance (gm, ) of 324mSmm7! and K-valve of 
418 mSV!mm'?., 

3.7.2 n-AlGaAs/InGaAs/GaAs HEMT 

We also fabricated n-AlGaAs/In, Ga,_, As/ 
GaAs pseudomorphic HEMTs using TBAs and 
arsine. A 15 nm Ing.,5Gag.g5 AS layer was grown 
for this structure, because the dependence of the 
InAs mole fraction on various properties of the 
pseudomorphic heterostructures, as shown in 
Fig. 34, implies that Ing ;;Gag.g;As is most 
suitable for the 15 nm channel layer. Figure 37 
shows the J/-V characteristics, measured at 
300 K, of n-AlGaAs/InGaAs/GaAs pseudomor- 


Arsine 


Source-drain current 


=f 0 ] 
Source-gate voltage Vg. (V) 
Fig. 38—Square-root of source-drain current versus 
source-gate voltage for pseudomorphic HEMTs. 
The source-drain voltage V gg is 1 es 
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phic HEMTs grown using TBAs. Excellent satu- 
ration characteristics were observed. 

In addition to the good electrical properties 
of selectively doped heterostructures, complete 
isolation of the active region is indispensable for 
the development of HEMT LSIs. A high-resis- 
tance buffer layer is required to obtain complete 
isolation. It is a significant issue whether the 
buffer layer grown using TBAs is resistive 
enough to isolate the adjacent transistor. 
Figure 38 shows the square root of the source- 
drain current (/Tas) plotted as a function of the 
source-gate voltage (V,,) for pseudomorphic 
HEMTs grown using a) TBAs and b) arsine. The 
source-drain voltage (Vqg,) was 1 V. Figure 38 
verifies that the pinch-off characteristics are 
excellent for both sources, which indicates 
that the TBAs-grown buffer layer has sufficient 
isolation capability for sub-micron gate HEMTs. 
For pseudomorphic HEMTs, the values of g,, 
and the K-value were 350mSmm7™! and 
480 mSV!mm'! for TBAs; and 308 mSmm! 
and 423 mSV!mm! for arsine. The results for 
pseudomorphic HEMTs confirm better perform- 
ance than conventional AlGaAs/GaAs HEMTs 
with both TBAs and arsine. The performance for 
TBAs was superior to that for arsine. However, 
the Si incorporation efficiency into AlGaAs and 
GaAs?®)-32) and the quality of the InGaAs 
layer?” were quite different between precursors 
having the same V/III ratio. We speculated that, 
since the performance of HEMTs depends on Ng 
of the AlGaAs layers and on the InGaAs purity, 
slight differences in HEMT performance should 
occur between group-V sources. Detailed re- 
search into InGaAs growth on GaAs with TBAs 
must improve the performance of pseudomor- 
phic HEMTs. 


4. Conclusion 

We have investigated undoped and Si-doped 
GaAs and AlGaAs, InGaAs, and selectively 
doped AlGaAs/GaAs_ heterostructures using 
TBAs in conventional or inverted horizontal 
atmospheric pressure reactors. EtAs was also 
used for comparison. 

The surface morphology of the GaAs and 
AlGaAs epitaxial layers grown using TBAs was 
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comparable to that of layers grown using arsine. 
PL measurements at 4.2 K revealed well-resolved 
emissions in the exciton region of the GaAs 
films which suggested that the GaAs films grown 
using TBAs and EtAs were of high quality, and 
that carbon incorporation into the GaAs epitax- 
ial layer was as low as achieved with arsine. 
X-ray analysis showed that AlGaAs epitaxial 
layers grown using TBAs were as good as those 
grown using arsine. Carbon incorporates into the 
AlGaAs epitaxial layers with a concentration of 
1015 cm™ to 10!® cm’. We interpret this as 
resulting from the stronger binding energy be- 
tween aluminum and carbon. However, increas- 
ing the V/III ratio reduces carbon incorporation, 
producing high-quality AlGaAs. In this study, 
the purity of EtAs was insufficient for AlGaAs 
growth. 

It was difficult to assure a uniform AlGaAs 
layer across the wafer using a conventional 
horizontal atmospheric pressure reactor. Gas 
phase reactions between TMGa and TBAs form 
adducts and other intermediates, causing the 
growth efficiency to vary markedly over the 
wafer surface. However, with an _ inverted 
horizontal atmospheric pressure reactor, we 
eliminated buoyancy-induced flow, thus reduc- 
ing the source deposition onto the reactor wall, 
and obtained high quality, highly uniform 
AlGaAs layers using TBAs. The uniformity of 
AlGaAs films using TBAs across a 3-inch sub- 
strate was within +1.2% for Ng and +0.9% for 
the layer thickness with an inverted horizontal 
reactor. 

Our study (bulk doping and planar doping of 
Si) revealed Si-doping to have seven significant 
key features: 

1) Using TBAs and EtAs enhances Si incorpora- 
tion. 

2) Diffusion of the doping component controls 
Si-doping using TBAs, even using silane. 

3) Using silane, Si-doping efficiency decreases 
as the V/III ratio increases. 

4) Using disilane, Si-doping efficiency decreases 
as the V/III ratio increases in the case of 
TBAs and EtAs, but increases in the case of 
arsine. At a higher V/III ratio, the doping 
efficiency using disilane is higher than that 
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using silane not only for arsine but also for 

TBAs. This increment is similar with TBAs 

and with arsine. 

5) The growth rate affects Si-doping efficiency. 

6) The dependence of Si-doping efficiency on 
the gas flow velocity is opposite for TBAs 
and for arsine, indicating that the gas phase 
reactions between As-species and Si-species 
dominate the Si-doping process and that the 
rate of reaction in the gas phase is different 
for different group-V sources. 

7) The dependence of Si-doping efficiency on 
substrate misorientation is the same for 
TBAs and arsine. 

To verify the relationship between the 
pyrolysis of As-species and the doping efficien- 
cy, we applied a kinetic simulation on group-V 
decomposition including the gas and surface 
reactions. The TBAs pyrolysis produces a large 
amount of As-species radicals. We propose that 
a large amount of silylarsine is formed by reac- 
tion between As-species radicals and SiH,. 
Silylarsine is predicted to be prominent when 
using organoarsine, and silylarsine rather than 
silylene is the most effective contributor to Si 
doping reactions. We also introduced the 
hindrance of Si-species adsorption on _ the 
surface due to AsH3. The mixture of enhance- 
ment and suppression could explain the compli- 
cated doping features. Using TBAs, EL2 concen- 
tration of GaAs decreases as growth temperature 
increases. This also suggests that AsH, forma- 
tion dominates the gas phase when using TBAs. 

Uniform Inp.,;Gag.g; AS on GaAs was ob- 
tained using TBAs at 630 °C. Increasing the InAs 
mole fraction, growth temperature, and V/III 
ratio resulted in poor uniformity, signifying an 
active vapor phase reaction between TMIn and 
TBAs. 

The electron mobility of the selectively 
doped n-AlGaAs/GaAs heterostructures 
with a 7nm spacer layer was as high as 
98000 cm?V~!s"!. The existence of 2DEG 
at the heterointerface was confirmed by the 
observation of Shubnikov-de Haas oscillations. 
Selectively doped n-AlGaAs/InGaAs/GaAs pseu- 
domorphic heterostructures were grown using 
TBAs or arsine, resulting in comparably good 
electrical properties. We demonstrated the first 
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fabrication of 0.5 um-gate n-AlGaAs/GaAs 
HEMTs and n-AlGaAs/InGaAs/GaAs_ pseudo- 
morphic HEMTs using TBAs. Sharp pinch-off 
and excellent saturation characteristics were 
observed for HEMTs made from either source. 
Vlas vs Vgs characteristics verified that the 
buffer layer had high resistance. We obtained 
high-performance HEMTs using either TBAs or 
arsine. These results demonstrate that GaAs, 
AlGaAs, and InGaAs layers grown using TBAs 
are very pure and are suitable for HEMT LSI 
fabrications. 
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Parallel processing techniques have been applied to the molecular dynamics calculation 
module of AMBER developed by the University of California. The target machine is the 
AP1000 distributed memory parallel computer developed at Fujitsu Laboratories Limited. 
AP1000 consists of up to 1024 processor elements connected with three different networks. 


To obtain a higher degree of parallelism and better load balance between processors, a 


particle division method was developed to randomly allocate particles to processors. Experi- 
ments showed that a problem with 41095 atoms can be processed 226 times faster with a 


512-processor AP1000 than by one processor. 


1. Introduction 

Molecular dynamic simulation is becoming 
an important technique in physics and chemistry 
and now plays a key role in protein engineering”). 
This technique will become more important as 
computer technology improves, eventually 
enabling the handling of larger molecules with 
added precision. Computing power is crucial to 
the analysis of proteins, and the speed and 
memory size of computers currently limit the 
size of proteins that can be studied. If we were 
also to consider the water molecules and counter 
ions surrounding a given protein and simulate 
biological phenomena, we would need a much 
faster machine. 

Molecular dynamics uses Newton’s equation 
of motion and performs an integration for each 
particle. This is very time consuming especially 
as the number of particles increases. Moreover, 
tens of thousands of integration steps are often 
necessary. For example, a 10000-step calcula- 
tion for a 40000-atom problem would take 
seven days on the large-scale FUJITSU M-780 
mainframe. The most time-consuming part is the 
force calculation, which accounts for over 90 
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percent of the total runtime. Vector processors 
have been used to speed up this calculation by 
factors of five to fifteen times, and have made 
simulation practical for small to medium-size 
problems. However, this is still insufficient for 
larger problems consisting of over ten-thousand 
particles. 

In this paper, we study a parallel processing 
approach to speed up the force calculation 
where certain problems such as communica- 
tion and load balancing must be solved. We 
parallelized a molecular dynamics module of 
AMBER??™?, a widely used program for analyz- 
ing molecular dynamics and mechanics of 
proteins and DNAs. AMBER was developed at 
the University of California in San Francisco. 
The machine used is a highly parallel computer, 
a 512-processor AP1000°)®, developed by 
Fujitsu Laboratories. AP1000 accommodates up 
to 1024 processor elements and has a peak 
performance of 8.53 GFLOPS. 


2. Molecular dynamics 


2.1 Molecular dynamics method 
In molecular dynamics, the force on each 
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particle is the sum of the forces from all mutual- 
ly interacting particles. Each particle is usually 
treated as a mass point, often having an electrical 
charge. 

The position and velocity of each particle 


are calculated by numerically integrating 
Newton’s equations of motion 
ad*x; 
Lie ie Pe (1) 


where F; is the force vector of particle i, m; is 
the mass of particle i, X; is the position vector of 
particle i, and ¢ is the time. Force vector F; is 
calculated as the summation of the gradient of 
potential LE. 


= EAE IXY}, 


where, X;;(t) is the distance between particle 7 
and particle j at time f¢. 
The potential energy function of AMBER is 


described as follows: 


Protas = % Ba t® — Ro ‘i 7 = ao (8 — 9)? 


+ > “4 11 0+ costae) 


dihedrals 2 
_ Ai Bij 4 UG 
i<j Ri'* Rij° eRij 
‘eer Des 
+> ( Bey ee (apie 


H-bonds Ry Ry 


The first three terms are for the energy of bond- 
ed atoms, and the remaining terms for non- 
bonded interaction by van der Waals or Coulomb 
energy (see Fig. 1). 

Newton’s equation of motion is time- 
integrated using the Leapfrog method”). The 
Leapfrog method is an alternate version of 
Verlet’s method’) and has been adopted for 
AMBER. It is described as below. 

First, the velocity of each atom is obtained 
by 

F; PG) 


vite y= 0(¢ S54 at (4) 
m; 


Then, coordinates are obtained by 


A 
X,(t + At) = X;(t) + AtV;(t +=). sae CD) 
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Fig. 1—Force field parameters. 


The time step Afr is usually 0.5 x 107!5s 
(0.5 femtoseconds), an increment obtained by 
considering the periodicity of the vibrating 
motion of bonds involving hydrogen atoms. 
However, this time step can be maximized to 
2 femtoseconds by stopping the stretching 
motion of the hydrogen atoms using the SHAKE 
method of Ryckaert et al?). 

The scaling method!® is used to maintain 
the pressure and temperature at desired values. 
Temperature is adjusted by scaling the velocity 
of atoms. Pressure is adjusted by scaling the 
length of each side of the basic cell that is the 
cubic region in which particles are held. 

Initial velocity is given according to the 
Boltzman distribution for a given temperature. 
Initial coordinates are given at those of a 
minimum energy structure obtained through 
molecular mechanics. 


2.2 Molecular dynamics module of AMBER 

The module used to execute molecular 
dynamics is called MOLDYN. Figure 2 shows 
the program flow of the original MOLDYN 
module. 

For calculating nonbonded forces, 
MOLDYN first makes a list of atom pairs, where 
the atom of each pair lie within the cutoff 
radius. This pair list generation is usually done 
every 20-30 time-steps. Bond connection 
information is used for bonded force calculation, 
and is constant throughout the calculation. Most 
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Nonbonded pair list generation 


Force calculation for nonbonded 
and bonded pairs 
Pressure control by scaling 
atom coordinates 
Time integration and SHAKE 
processing 


Time-step loop 


Fig. 2—Flow of original MOLDYN. 


Other 3.07 % Th P 
Pair list generation 


5.19 % 


Force 
(bonded) 
0.43 % 


20 665 atoms 


200 steps 

112 min 

Force Pe Cutoff 1.0 nm 
(nonbonded) Pair list generation/ 
91.31 % 25 steps 


Fig. 3—Timing result of original MOLDYN on M-780. 


of the computing time is spent calculating non- 
bonded forces. Ninety-one percent of the calcula- 
tion time is spent executing 200 time-steps of 
dynamics for a 20 665-atom problem (see Fig. 3). 


3. Parallelization of molecular dynamics 
3.1 AP1000 architecture 

Figure 4 shows the architecture of the 
AP1000 and Table 1 lists the specifications. 
Processor elements, called cells, are connected 
with three networks, called the broadcast 
network (B-net), torus network (T-net), and 
synchronization network (S-net). A SUN 
4-workstation host controls the AP1000. 

The host, or any cell, can broadcast data to 
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mw nchronization network (S-net) 


buffers and 
RTC hard disks 
Two-dimensional torus network (T-net) 


Cell : Cell processor 
RTC: Routing controller 


BIF : Broadcast network interface 


Fig. 4—AP1000 architecture. 


Table 1. AP1000 specifications 


Item Specifications 
Architecture Distributed memory, 
System MIMD parallel computer 
Number of 64-1 024 cells 
processors (0.53-8.53 GFLOPS) 
Structure | Tree 
S-net Function Barrier synchronization 
and satus (40 events) 
Structure Ring + hierarchical bus 
Woes (50 Mbyte/s) 
Function Broadcast, scattering, 
and gathering 
Structure Two-dimensional torus 
(25 Mbyte/s) 
Routing Wormhole routing 
T-net method with structured buffer 
pool 
Special function] Broadcast to specified 
cells 


MIMD: Multiple instruction stream multiple data stream 


all cells through the B-net, which is implemented 
with a ring and hierarchical connections. It also 
scatters and gathers data to and from cells 
directly by hardware. The broadcast network 
interface (BIF) chip handles these functions. 
The peak transfer rate is 50 Mbyte/s. 

The T-net is used for communication 
between any two cells. The topology is a two- 
dimensional mesh with boudary wires connected 
to the other sides. The routing controller (RTC) 
on each mesh point performs automatic message 
routing using wormhole routing and a structured 
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SPARC 
IU + FPU 


BIF : B-net interface 

FPU : Floating point unit 
MSC: Message controller 
RTC: Routing controller 


Fig. 5—Cell processor configuration. 


buffer pool algorithm. The transfer speed is 
25 Mbyte/s per channel. 

Figure 5 shows the configuration of a cell 
processor, which consists of a 25 MHz SPARC 
integer unit (IU) with a Weitek 3170 floating- 
point processor, 128 Kbyte of cache memory, 
and communication devices. The message 
controller (MSC) supports DMA transfer func- 
tions. 

The AP1000 software runs under Cell OS® 
which resides in all cells and schedules multiple 
tasks in cells. The application software is written 
in C and FORTRAN, and communication and 
synchronization functions are provided as 
callable functions. An application program con- 
sists of the host program and various cell tasks. 


4. Implementation of particle division method 

The MOLDYN module was parallelized using 
the particle division method which partitions the 
ensemble of atoms into groups and assigns them 
to processors. This method makes it possible to 
use more than several hundred processors for 
large problems which contain more than several 
thousand atoms. A degree of parallelism (the 
number of atoms) much greater than the 
number of processors is used to balance the load 
by assigning randomly selected multiple atoms 
to processors. The calculation load is evenly 
distributed in a statistical sense by this alloca- 
tion method. 
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Fig. 6—Program structure of MOLDYN. 


4.1 Structure of the parallelized MOLDYN 

Figure 6 shows the program structure of 
AMBER converted for use on parallel machines. 
The original MOLDYN program was split into 
host and cell task programs. The host program 
sets up initial data and broadcasts them to all 
cells. The time-step loop is performed in cells 
without host program intervention. 


4.2 Data division and allocation. 

Data is broadcast from the host program to 
all cells. Each cell selects the data necessary for 
itself and discards any unnecessary data. Model 
data consists of protein molecules and many 
water molecules as solvent. We chose residues of 
proteins and water molecule as units of data 
allocation because they are treated as units 
in MOLDYN. These are called submolecules. 
Residues and water molecules are distributed 
to cells in such a way that each cell holds an 
equal number of residues and water molecules. 

A bond list of atoms is read at the beginning 
of calculation from the input file and distributed 
to cells. This information does not change 
throughout the molecular dynamics calculation. 


4.3 Nonbonded pair list generation 


A nonbonded pair list is generated for 
assigned residues and water molecules. A cutoff 
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Fig. 7—Selection of submolecule pairs (original). 


distance check is not done between individual 
atoms but instead is done between submolecules 
(residues or water molecules). If the minimum 
of- the distances of all atom pair between two 
submolecules is less than the cutoff radius, all 
the pairs are registered in the pair list. This is 
called residue-based pair list generation and is 
often used in the MOLDYN program. 

In the original code, the pair-list is generated 
by nested DO loops indexed by submolecule 
numbers as shown below (see Fig. 7). 


DOI=1,NR 
DO J=1,NR 


Check minimum distance between 
submolecules I and I 
IF distance < cutoff 
THEN generate pair list between 
atoms in two submolecules 


END DO 
END DO 
where NR is the total number of residues and 
water molecules. 
Before parallelizing MOLDYN, this code is 
changed to 
DOT=1,NR 
DO J=1,NR 
Check minimum distance between 
submolecules I and J 
IF distance < cutoff 


THEN generate pair list between 
atoms in two submolecules 
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Fig. 8—Selection of submolecule pairs (parallel version). 


1 


Residues 


Water 
molecules 


END DO 

END DO 

When this nested loop is executed, each sub- 
molecule pair is checked twice. We avoided this 
using a selection rule determined by considering 
even load distribution. It is shown in Fig. 8. 

In Fig. 8, circles show pairs of submolecules 
I and J. Pairs shown as black circles are selected 
and white ones are discarded for distance check- 
ing. For residue-residue and water-water pairs 
(circles in the top-left and right-bottom rectangle 
areas in Fig. 8), even numbered pairs (I + J is an 
even integer) are selected when I is less than J, 
odd numbered pairs (I+J is an odd integer) 
are selected when I is greater than J. When I is 
equal to J, the pair is always selected. 

The changed code is then parallelized by 
breaking up outer DO loop index I into 
segments and assigning them to processors. In 
other words, each processor executes a different 
set of rows in Fig. 8. 

The distance check of submolecules I and J 
is done by the processor assigned to the sub- 
molecule. If the distance is less than the cutoff, 
atom pairs between the submolecule pair are 
registered in the pair list in each cell for later 
force calculation. Cyclic allocation of submole- 
cules to processors (submolecule I is assigned to 
processor K =I] modulo P, where P is the total 
number of processors) allows the number of 
nonbonded atom pairs to be distributed fairly 
evenly. 
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The coordinate data of all atoms are stored 
in each cell for pair list generation (performed 
every 20-30 time steps) and force calculation. 
Coordinate data of all atoms can be collected 
efficiently by broadcasting data assigned to each 
cell to all other cells through the B-net. As a 
result, each cell holds a copy of the coordinate 
data for all atoms. The amount of coordinate 
data sent is 48 bytes (6 double precision words) 
per atom. 


4.4 Force calculation and summation 

The force between two atoms is accumulat- 
ed to form the force data for each individual 
atom. One of these two atoms is possibly assign- 
ed to another cell. After force calculation is 
finished, this force data is sent to cells to which 
interacting atoms are assigned (see Fig. 9). This 
point-to-point (cell-to-cell) communication is 
done using the T-net. Force data to be sent to 
the same atom is summed up before being sent. 


Force calculation 
and summation 


Cell A 


Force summation 


Cell B 


Atom j 


O 


Receive Fi; 


Atom 7 
O 


Calculate Fi); 
Ppa Ppceelty 


Bye Py iy 


T-net 


Fig. 9—Force calculation and summation. 


This reduces the amount of data sent. 

There are 24 bytes (3 double-precision 
words) of force data per atom. The amount of 
force data sent to other cells depends on the 
number of interacting atoms allocated to other 
cells. 


4.5 Visualization 

The AP1000 can be used to display simula- 
tion results. A parallel display program was also 
developed to show atoms and bonds as full-color 
images. 

The visualization software consists of two 
types of tasks. OBJECT tasks receive geometric 
data from MOLDYN tasks and generate images 
of shaded spheres and lines with z values 
(distance from the eye point). Generated images 
are redistributed to IMAGE tasks which remove 
hidden surfaces using the z-buffer algrorithm. 
Image space is divided and processed in the 
IMAGE tasks!» 1?) (see Fig. 10). Running this 
program concurrently with MOLDYN makes 
possible the visual simulation of molecular 


Host program 


Initial data Image data 


Cell | | MOLDYN OBJECT 
Cell 0 | | | 
Molecular Image Hidden 


Fig. 10—Visual simulation program. 


Table 2. Test model data 


: Tess monte Lysozyme-A Lysozyme-B Hemoglobin 
Data item 
Protein Lysozyme Lysozyme Hemoglobin 
Atoms (total) 20665 41095 88 384 
Atoms (protein) 2041 2041 9 082 
Atoms (water) | 18624 39024 79302 
Residues 130 130 578 
Water molecules 6 208 13 008 26 434 
Box size (mm) 7.03 x 5.89 x 5.44 T5351 x 7.51% 751 41.26: x 11.26% 11226 
Cut-off radius (nm) 1.00 1.00 1.00 
Cut-off checking residue base residue base residue base 
Pair list generation interval every 25 steps every 25 steps every 25 steps 
Number of nonbonded pairs (initial) 4711181 8 965917 18 888 552 
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dynamics. Computer animation can easily be 
produced by recording generated images. 


5. Experiments 
5.1 Test data 

We used three problems to evaluate the 
performance of MOLDYN. Two of them contain 
a human lysozyme and the other a hemoglobin 
at the center of a pressure control cubic box 
filled with water molecules. Human lysozymes 
consist of 130 residues and hemoglobin 578. 
The sizes of cubes and the number of water 


molecules are different between these three data 
items (see Table 2). These problems were also 
solved using an M-780 for comparison. Average 
execution times per time step are listed in 
Table 5. 


5.2 Results and performance analysis 

We experimented with these model problems 
using the AP1000 with 512 cells. Tables 3 and 4 
list the timing results for a single time step. 
Experiments using less than 512 cells were also 
done using subsets of the 512 cells. Communica- 


Table 3. Timing results of Lysozyme-A 


aston nares aeigatial (A 16* 32 64 128 256 512 
Pair list generation (s) 971.80 60.74 30.37 16.56 8.65 4.96 231 
Nonbonded force calculation (s) 265.30 16.58 | 8.29 4.15 2.08 | 1.04 0.52 
Other calculation (s) 10.60 0.66 0.33 0.20 0.13 0.07 0.04 
Communication (total) (s) 0 0.79 0.77 0.72 0.68 0.78 | 0.96 
Coordinate data (s) 0 0.30 0.30 0.31 0.32 0.35 0.41 
Force summation (s) 0 0.40 0.38 0.31 0.24 0.21 0.13 
Other (s) 0 0.09 0.09 0.10 0.11 0.22 0.42 
Idle time (s) | 0 0.50 0.49 0.55 0.30 0.36 0.32 
Step time (with pair list gen.) (s) 1 247.70 79.27 40.21 22.18 11.80 7.21 4.16 
Step time (s) 275.90 18.53 9.89 5.62 3.19 2.25 1.85 
Average step time (s) 314.77 20.96 11.10 6.28 3.53 2.45 1.94 
Speed-up ratio 1.00 15.02 28.35 50.10 89.06 iL 128.56 162.05 
Efficiency (%) 100.00 0.94 0.89 0.78 0.70 0.50 0.32 
Table 4. Timing results of Lysozyme-B 
re Serer) 4s 16* 32" 64 128 | 286 512 
Pair list generation(s) | 3 764.00 235.00 118.00 58.81 30.20 il 11.56 8.24 
Nonbonded force calculation (s) 506.00 31.60 15.80 791 3.96 Loe 0.99 
Other calculation (s) 19.20 1.20 0.60 0.31 0.19 0.10 0.07 
Communication (total) (s) 0 1.48 1.41 1.36 LT Lag 1.39 
Coordinate data (s) 0 0.60 0.60 0.61 | 0.62 0.64 0.70 
Force summation (s) o| o80 | 0.73 0.59 0.42 0.31 0.32 
Other (s) 0 0.08 [ 0.08 0.17 0.12 | 0.24 0.38 
Idle time (s) 0 0.80 0.78 | 0.50 0.38 0.32 0.22 
Step time (with pair list gen.) (s) | 4289.00 270.00 I 137.00 68.89 35.90 15.14 10.90 
Step time (s) 525.00 35.08 18.59 10.08 5.70 3.58 2.66 
Average step time (s) 675.56 44.48 23.33 12.43 6.91 4.04 | 299 
Speed-up ratio 1.00 15.19 28.96 | 54.34 97.79 167.12 225.97 
Efficiency (%) 100.00 0.95 0.91 0.85 0.76 0.65 0.44 


*: The data for the number of cells indicated was not executable due to memory limitations; values are estimated by 


extrapolating experimental data. 
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Fig. 11—Timing results of Lysozyme-B and Hemoglobin. 
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Fig. 12—Performance of MOLDYN on the AP1000. 


tion time was measured for coordinate data, 
force data, and other communication for global 
summation of scalar values such as energy, 
weight, and minor communications. 

Figure 11 shows the timing analysis of 
Lysozyme-B and Hemoglobin. Pair list genera- 
tion and nonbonded force calculation time were 
efficiently decreased by increasing the number 
of cell processors. When a small number of cells 
are used, each atom interacts with atoms in 
almost every cell. This is because we allocated 
atoms randomly to processors irrespective of 
their coordinates. The time for sending force 
data decreases slightly as the number of cells 
increases. This shows that the number of atoms 
for which force data is sent decreases as the 
number of cells increases. 

Communication — time 
increases slightly with the 


for coordinates 
number of cells 
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Table 5. Performance comparison 


Test model | Lysozyme-A | Lysozyme-B | Hemoglobin 
Test Time(s) Time(s) Time(s) 

system (Prfm. ratio) |(Prfm., ratio) | (Prfm. ratio) 
M-780 33.57(1.00) 66.95 (1.00) |100.00(1.00) 
AP1000-64 | 6.28 (5.35) [12.43 (5.39) NA | NA 
AP1000-128] 3.53(9.51) | 6.91(9.69) | NA | NA 
AP1000-256| 2.45(13.70)| 4.04(16.57)| NA | NA 


AP1000-512 


1.94 (17.30) 


2.99 (22.39) 


7.02 (14.25) 


Fig. 13—Example of a generated image (Lysozyme-A). 


ball-and-stick models. 


buttons for controlling molecular dynamics 


The right hand side shows an image of the 
configuration of a human lysozyme molecule 
(center) with water molecules displayed with 


The left 


hand side are 


simulation and visualization. Users can start and 
stop molecular dynamics simulation and can 
change the point of view, select display items 
(protein or/and water, main chain only or all, 
etc.), and display styles (ball-and-stick or stick 


only 


Ds 


because of an increase in the broadcast com- 


munication 


overhead. The amount of 


data 


received by each cell is constant because coordi- 


nate data 


communication 


for all 


the 


atoms are collected. This 
reduces 


performance 


improvement when over 256 cell are used (see 


Fig, 12). 


The 512-cell AP1000 executes Lysozyme-A 
17.3 times faster the M-780 (see Table 5), and 


Lysozyme-B 


runs 


22.4 


times 


faster. One 


problem is the large memory requirement for 
coordinate data. Model data containing more 
than 100000 atoms cannot be processed. This 
problem is left for future work. 
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Figure 13 shows a generated image for a 
Lysozyme-A model containing 20 665 atoms. 


6. Conclusion 

This paper examined parallel processing of 
molecular dynamics. The AMBER -molecular 
dynamics calculation program was parallelized 
using a particle division technique for execution 
by a distributed memory parallel computer 
consisting of several hundred processor elements. 
The parallelized program was implemented on 
an AP1000 and tested with reasonably large 
problems. Using 512 processors, we improved 
performance by factor of 162 fora 20665 atom 
problem and 226 for 41095 atom problem. 
Performance was respectively 17.3 and 22.4 
times faster than on the FUJITSU M-780 main- 
frame. 

Two problems remain to be solved. One is 
improving performance for more than 256 
processors. The sequential part of coordinate 
data communication must be minimized. The 
other is to reduce the memory requirements for 
coordinate data. To solve these problems, we are 
studying a method that avoids copying coordi- 
nate data for each cell. 
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This paper discusses an automated visual inspection system for bonded IC wires. The system 
uses high-contrast image capture, an algorithm for accurately measuring the bonding balls, 
and a wire inspection algorithm. The high-contrast image capture system consists of a 
wide-area camera and two types of lighting optics, one for ball sensing and the other for 
wire sensing. The algorithm for. measuring the bonding balls is based on morphological 
techniques, and the wire inspection algorithm is based on a border following algorithm. 
The automated inspection system measures ball diameters to within +5 zm accuracy, which 
corresponds to half a pixel of the captured picture. The system takes 0.2 seconds to inspect 
a wire and ball. Combining the inspection system with an automatic wire bonder enables 


a fully automatic bonding system. 


1. Introduction 

Automated visual inspection is vital to wire 
bonding to ensure a high product yield. In the 
wire bonding process, gold wire 20 to 38 wm in 
diameter is bonded to IC pads and package 
leads. Gold balls 70 to 110 um in diameter are 
bonded to IC pads ranging in size from 90 x 90 
to 120 x 120 um?. One IC chip may have as 
many as 300 bonding points. A wire bonder 
bonds one wire within 0.2 seconds. 

Visual inspection can detect defects in a 
wire and ball, and can evaluate bonding quality. 
The detected wire defects include broken 
wires, wires that are too close together and 
incorrect wiring paths. The ball defects that 
can be detected include deformation and incor- 
rect placement. The inspection can also measure 
bonding location and ball deformation, which 
reflect the performance of a bonder. Thus, the 
inspection can find symptoms which cause the 
performance of the bonder to deteriorate. The 
detected symptoms are analyzed and bonding 
parameters, such as bonding location and force, 
are modified. This prevents a large number of 
defective ICs from being produced. Because the 
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location tolerance is within +10 wm, a measure- 
ment accuracy of +5 jm is needed. 

The bonding quality can be monitored 
by combining the inspection system and bond- 
ing system in-line. This guarantees the bonding 
quality, resulting in improved yield. 

To achieve the objectives, a visual inspection 
system must meet two basic requirements: 

1) Accurate evaluation of bonding quality 

The ball measurement must be accurate 
to within +5 um. 

2) Real-time inspection 

The inspection of a wire and ball must 
be completed within 0.2 seconds to keep up 
with the wire bonder. 


2. Background 

Three technical problems must be solved: 

1) Obtaining clear images of both the wires and 

balls (see Figs. 1 and 2) 

The gold balls are bonded onto aluminum 
pads. Because both the balls and the pads are 
shiny, it is difficult to distinguish the ball 
images from the pad images. 

2) Measuring ball size and shape 
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10 mm 


Wiring defects 


Ball defects 


Fig.2— Typical wire and ball defects. 


Wire images overlapping the ball patterns 
make it difficult to measure ball shapes accu- 
rately. 

3) Inspecting wire shapes 

Wire paths are not completely straight, 
and their deviation makes it difficult to measure 
the distance between the wires. 

To solve these problems, the authors de- 
veloped high-contrast image capture, an accurate 
ball measurement algorithm, and a versatile 
wire shape inspection algorithm. 


3. Image capture system 

The authors developed a_ high-contrast 
wire and ball image capture system consisting 
of a wide-area camera and two types of lighting 
optics. 
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Table 1. Megacamera specifications 


Item Specifications 

Sensor 2 048-element CCD linear sensor 
Resolution 10 um/pixel 

Density 256 gray levels 

Area 2048 x 2048 pixels (20 x 20 mm’) 


Input time 2s 


Megacamera 


Vertical _ 
illumination 


Side 
illumination 


20 mm 


Se 


20 mm 


Fig.3— High-contrast image capture system. 


3.1 Wide-area camera 

A wide-area camera called a megacamera’ 
can capture a 2048 x 2048 pixel image using 
a mechanically operated CCD linear image 
sensor. Table 1 shows the specifications of 
this camera. The resolution is 10 um, and the 
capture area is 20 x 20 mm7?. This camera can 
capture an image in two seconds. 


) 


3.2 High-contrast lighting system 

High-contrast imaging is essential for correct 
inspection. The system uses two types of light- 
ing optics. They are vertical illumination and 
side illumination (see Fig. 3). 

Vertical illumination is used for ball pattern 
sensing. The light goes down to the pad and 
ball vertically. While the pad is flat, the ball is 
round, thus only the pad reflects the light back 
to the camera. This makes the ball image dark, 
and the pad image light. 
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Side illumination is used for wire sensing. 
The side illumination apparatus consists of 
optical fibers and a halogen lamp. Figure 4 
shows the side illumination optics. Each optical 
fiber is arranged in a circle facing inward. 
The camera is set vertically above the side 
illumination optics and views wires through 


Camera 


Lighting system 


Lighting angle VY 


Fig.4— Side illumination optics. 
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Lighting angle (degrees) 


Fig.5— Lighting characteristics. 


a) Vertical illumination 


the circle. The wires are illuminated from all 
sides. The pad does not reflect light, but the 
wires do, regardless of their direction. 

Figure 5 shows the lighting characteristics of 
the side illumination optics. The wire image 
clarity depends on the lighting angle. This 
angle is defined as the angle between lighting 
and the horizontal plane. The solid line is a 
clear image boundary, and is based on experi- 
ment. As the lighting angle decreases, the range 
of clear images increases. Usually, the wire 
inclinations range from O to 40 degrees. The 
lighting angle was therefore set at 10 degrees 
to cover the whole angle. 

Figure 6 shows the signal-to-noise ratio 
(SNR) of the wire images. The signal-to-noise 
ratio of wire images also depends on the lighting 


20 
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Fig.6— SNR of wire-image capture. 
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b) Side illumination 


Fig.7— Ball and wire image. 
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angle. In this case, noise is caused by a chip’s 
surface pattern. As the lighting angle decreases, 
the signal-to-noise ratio increases. This result 
indicates that a lighting angle of 10 degrees is 
sufficient for high contrast at an acceptable 
signal-to-noise ratio as well as wire image clari- 
ty. 

Figure 7 is IC image taken using this illumi- 
nation system. In Fig. 7a), only vertical illumi- 
nation was used. The ball patterns are dark, 
while the pad patterns are light. In Fig. 7b), 
side illumination was used. The wires are clearly 
visible against the dark background. 


4. Inspection algorithm 
4.1 Ball measurement algorithm”? 

4.1.1 Requirements for ball measurement 

There are two requirements for ball measure- 
ment. The first is accurate edge detection from 
a low-resolution image. The resolution of the 
original image is 10 wm. The required accuracy 
for ball measurement is +5 wm, which corre- 
sponds to half a pixel on the original image. 

The second requirement is ball outline 
reconstruction. With vertical illumination, as 
shown in Fig.8, the wire pattern overlaps 
the ball. This wire pattern must be eliminated 
before the ball pattern diameter can be meas- 
ured. 

4.1.2 Edge detection 

Subpixel accuracy is needed for edge detec- 


Wire 
Ball 


Diameter 


Fig.8— Side and top view of ball and wire. 
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tion®”*), Figure 9 shows a ball edge detection 


technique the authors developed. 

The original image, which has a resolution 
of 10 um, is sensed at 256 gray levels. A finer 
edge pattern is produced by _ interpolating 
pixels between pixels. Given digital error, 
the intervals of interpolated pixels must be 
2.5 um to measure ball diameters within a 
+5 um accuracy. The gray-level value of the 
interpolated pixels is calculated by the gray 
levels of the four nearest neighbors. Applying 
a Laplacian operator to the interpolated image, 
it is possible to produce a second-derivative 
image. Zero-crossing points”? in the second- 
derivative image indicate edge locations, and 
are detected from the ball center. 

These zero-crossing points correspond to ball 
edges with a +2.5 wm accuracy. As a result, the 
diameter can be measured within an accuracy 
of +5 um, 

4.1.3 Reconstruction of ball outlines 

Preprocessing for ball outline detection is 
shown in Fig. 10. The orthogonal coordinates 
of the ball image are transformed to polar 
coordinates. The origin is at the center of the 
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Fig.9— Edge detection algorithm. 
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Fig.10—Ball outline detection preprocessing. 
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ball. In Fig. 10, plane envelopes are made to 
correspond to ball outlines. The jutting envelope 
is a wire outline. The solution to the problem of 
deleting the wire pattern is to delete the large 
convex of the envelope. 

To distinguish the ball outline from the wire 
image, a morphological transformation®’?) is 
used, called opening. Opening compares an 
object pattern with a disk structuring element 
and removes the sharp convex part where the 
outline does not fit the structuring element. 

Figure 11 shows ball outline reconstruction. 
Applying this to polar coordinates, the disk 
structuring element rolls beneath the envelope 
of the ball outline. The transformed edge 
in polar coordinates is transformed back to 
orthogonal coordinates, reconstructing only 
the ball outline. The ball shape is then inspected 
by measuring the ball diameters in all directions. 


4.2 Flexible wire shape inspection algorithm 
The authors used the border following 
technique to inspect the wire shape (see Fig. 12). 
Masks consisting of a 3 x 3 pixel matrix detect 
the right dnd left border locations. The current 
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Fig.1 1—Morphological ball outline reconstruction. 
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Fig. 12—Wire inspection algorithm. 
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left border position is pixel L. The next border 
position is detected by checking the contents 
of the pixels in the matrix. The pixel is checked 
counterclockwise starting from the last border 
pixel L’. Position 5, where the space changes to 
the pattern, is detected as the next border 
position in this case. The position of the right 
border is detected the same way in the clockwise 
direction. 

The wire width, the distance between 
neighboring wires, and deviation from the 
standard line are calculated during border 
following. The distance is checked against the 
standard at every line, and wires that are either 
shorted or too close together are detected. 
Deviation is checked to detect broken wires, 
lifted balls, and lifted crescents. Missing wires 
are detected at the beginning of border follow- 
ing. 


4.3 Results 

The results of the ball measurement algo- 
rithm and wire inspection algorithm are given 
below. 

Figure 13 shows an original image and its 
interpolated image. The original image, Fig. 13a), 
is a 256-level image of the bonding ball on the 
pad, magnified four times. The resolution of 
the original image is 10 wm. The interpolated 
image, Fig. 13b), is obtained by using bilinear 
interpolation. The gray-level values on every 
2.5 wm grid are computed from the values of 
the neighboring pixels. 

Figure 14 shows an example of the ball 
outline reconstruction process. Figure 14a) is 
the gray level representation of the interpolated 


a) Original image b) Interpolated image 


Fig. 13—Pixel interpolation. 
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image. The resolution is 100 x 100, coded in 
8 bits. Figure 14b) shows the second-derivative 
image, obtained by applying a 3 x 3 Laplacian 
mask to Fig. 14a). Figure 14c) shows the zero- 
crossings detected at quarter-pixel accuracy. 
These zero-crossings indicate the pad, wire, 
and ball outlines. Figure 14d) shows the coordi- 
nate transformation of the ball zero-crossings. 
The radius vectors represent the lines radiating 
from bell center to the zero-crossings on the ball 
outline. Figure 14e) shows the results of ball 
outline reconstruction, obtained by applying 
morphological opening operators to the polar 
coordinates. The cross in Fig. 14e) is the ball 
center detected from the ball outline. Finally, 
the ball diameters are measured by using the 
reconstructed outlines and the detected center. 

The experimental results for ball measure- 
ment indicate that the measurement algorithm 
measured the ball diameters to an accuracy of 
+5 wm (see Fig. 15). 

Figure 16 shows some examples of defective 
wires detected by the wire inspection algorithm. 


Fig. 14—Ball outline reconstruction. 


a) Broken wire 


In these applications, three wires are inspected 
at a time. While following the wire borders, 
the wire width and the distance between wires 
are measured. Figure 16a) shows a broken wire. 
The mark indicates where the wire width be- 
came zero in the middle of the wire. Figure 1 6b) 
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Fig.15—Ball measurement accuracy. 
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b) Shorted wires 


Fig.16—Detected wiring defects. 
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Fig. 17—System configuration. 


shows shorted wires. The mark indicates where 
the distance between neighboring wires became 
zero. 


5. Automated inspection system 
5.1 System configuration 

In the system configuration, the mega- 
camera captures the wire image by side illumi- 
nation and the ball pattern by vertical illumina- 
tion (see Fig. 17). The video signal is converted 
by an 8-bit analog-digital converter. Of the 
total 8-Mbyte image memory, 4 Mbytes are 
used for wire images and 4 Mbytes for ball 
images. The system is controlled by 32-bit 
microprocessors and co-processors. The wire 
bonder and inspection system exchange bond- 
ing addresses and bonding quality. The wire 
bonder parameters can be adjusted based on 
the bonding quality. 

Figure 18 shows the wire inspection system. 
A wire bonder can be combined with a frame 
feeder. A CRT displays inspection results. 


5.2 Specifications 

The system is used to inspect bonded 
wires on lead-frame ICs and can detect wires 
which are broken, too close together or shorted. 
The system can also detect incorrect wiring, 
deformed balls, and incorrect ball placement. 
Bonding reliability is calculated by measuring 
the size and shape of the deformed balls. The 
measurement accuracy is +5 wm which corre- 
sponds to half a pixel. It takes 0.2 seconds to 
inspect each wire (see Table 2). 
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Fig. 18—Wire inspection system. 


Table 2. Specifications 


Object Bonded wires on lead frame IC 

Defects Wires that are broken, too close to- 
gether, or shorted; incorrect wiring; 
deformed balls; and incorrect ball 
placement 

Quality Ball size and shape 

Accuracy +5 um 

Time 0.2 s/wire 


6. Conclusion 

This automated visual inspection system uses 
high-contrast imaging, an accurate measurement 
algorithm, and a versatile wire shape inspection 
algorithm. It can inspect a wire in 0.2 seconds 
and can work in unison with the wire bonder. 
A fully automatic bonding system can be 
constructed by combining this inspection system 
with a wire bonder. 
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This paper discusses a low drive voltage and low chirp modulator integrated DFB laser 
(MI-DFB laser). In this research, the composition of the absorption layer was adjusted exper- 
imentally to obtain a low drive voltage, and the absorption layer thickness was analytically 
optimized. The drive voltage for an extinction ratio of -13 dB is found to be -3 V, which is 
small enough for high-speed IC drivers. A SiN film on the front facet is optimized for an 
antireflective (AR) coat, and gives a low chirp width of 0.016 nm at -15dB (0.004 nm at 
-3 dB) under 10 Gbit/s NRZ pseudo random modulation. The low drive voltage and low 
chirp performance is promising for use in the next generation of multiple-gigabit long-haul 


optical fiber transmission systems. 


1. Introduction 

Current fiber optic transmission systems 
directly modulate the intensity of DFB lasers. 
Directly modulated DFB lasers exhibit large 
wavelength chirp caused by relaxation oscilla- 
tion. This large chirp limits the bit-rate and span 
product because of the dispersion of a fiber. 
But, due to the need to transmit increasing 
amounts of information, multiple-gigabit optical 
fiber transmission systems are in demand. To 
achieve multiple gigabit transmission over large 
distances, low-chirp light sources are essential. 
An MI-DFB laser is a promising device for re- 
ducing the chirp because it has an external 
modulation scheme which gives very low chirp. 

The authors previously reported the MI-DFB 
laser having an absorption layer whose photo- 
luminescence (PL) wavelength was 1.4 pm!?, 
The device had a drive voltage of -5 V and a 
chirp width of 0.01 nm at 3 dB down. This drive 
voltage of -5 V is high for practical use. Al- 
though a chirp width of 0.01 nm is lower than 
that in a directly modulated laser (LD), it is not 
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sufficient for long-haul and multiple-gigabit 
transmissions. Therefore, a reduction in the 
drive voltage and in the chirp are important for 
multiple-gigabit optical fiber transmissions. 

This paper describes the design for a low drive 
voltage and a low chirp width, and then dis- 
cusses the performance of the improved MI-DFB 
laser. 


2. Structure 

A schematic structure of the device is shown 
in Fig. 1. The monolithic device has three re- 
gions on one n-InP substrate. To operate the 
device, a constant current is injected into the 
DFB laser (EA-MOD) region and a modulated 
voltage is applied to the electro-absorption 
modulator (DA-MOD) region. The DFB-LD 
region is 300 wm long and the EA-MOD region is 
200 um long. A buffer region was introduced 
without a contact layer between the DFB-LD 
region and the EA-MOD region for electric 
isolation. A good optical coupling between 
the DFB-LD region and the EA-MOD region was 
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Fig. 1—Schematic structure of the MI-DFB laser. 


AR coat 


achieved by the butt-joined waveguide. The 
stripe was buried with semi-insulating InP for 
low parasitic capacitance!). 

The front facet was coated with a SiN film 
to reduce reflectivity. The reduced reflectivity 
gives a flat frequency response and low chirp 
width without any influence from the relaxation 
oscillation of the laser region”). The back facet 
was coated with two pairs of SiO, film and an 
amorphous Si film for a high reflective (HR) 
coat to achieve high slope efficiency and large 
side-mode suppression. 


3. Design for low operating voltage 

The EA-MOD uses the Frantz-Keldysh ef- 
fect, which is photon absorption assisted by 
tunneling electrons under an electric field®?. 
This effect is enhanced by both small detuning 
between the photon energy and the band gap 
energy, and a large applied field. Taking this into 
consideration, the authors designed the compo- 
sition and the thickness of the absorption layer 
to obtain a low drive voltage. 

Detuning between the light energy and band 
gap energy introduces a trade-off between high 
output power and low drive voltage. This is be- 
cause as detuning decreases, the absorption loss 
under no applied voltage increases. 

Firstly, the composition of the absorption 
layer was chosen considering the output light 
power. The authors experimentally evaluated 
the absorption loss of a solitary EA-MOD whose 
both facets were cleaved. The reflectivity of the 
facet was estimated to be 31 percent. The length 
of the waveguide was 200 um. The absorption 
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Fig. 2—Loss of waveguide under no applied voltage. 


layer had a PL wavelength of 1.4 wm. Transmis- 
sion light power was measured through the EA- 
MOD waveguide changing the incident light 
wavelength with an F-center laser. The absorp- 
tion loss was calculated from this spectrum using 
Hakki and Paoli’s method’. 

These results indicate that the absorption 
loss increases as the detuning between the band 
gap energy of an absorption layer and light ener- 
gy decreases (see Fig. 2). A PL wavelength of 
1.40 um for a light wavelength of 1.55 um 
corresponds to a detuning of 0.08 eV, and is 
expected to have a loss of 30 percent. For this 
result and the previous data of the 1.4-um 
absorption layer), the average effective output 
power from the DFB region is expected to be 
14mW at 150mA. When the detuning is 
0.067 eV, the output power from DFB region is 
halved to 7 mW at the front facet. 

If detuning is smaller than 0.067 eV, the 
output power is too small for practical use. 
Detuning of 0.067 eV is equivalent to a PL 
wavelength of 1.43 wm for light wavelength of 
1.55 wm. Because the detuning for a PL wave- 
length of 1.43 um is smaller than that for a 
1.4 wm absorption layer which has a drive volt- 
age of -5 V, the drive voltage of EA-MOD having 
a 1.43-um absorption layer is expected to be 
lower than -5 V'). Thus the PL wavelength of 
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Fig. 3—Dependence of extinction ratio on the thickness of 
the absorption layser. 


the absorption layer was fixed to 1.43 um®?. 

The authors then calculated the dependence 
of the extinction ratio on the absorption layer 
thickness in an EA-MOD taking account of the 
applied electric field and optical field. The 
device model for this calculation was solitary 
two dimensional EA-MOD. The thickness and 
PL wavelength of the InGaAsP guide layer 
model was assumed to be 0.15 wm and 1.1 wm 
respectively. Concentration of the p-clad layer 
was 5 x 10!7 cm, and one of the n-substrates 
was 2x 1018 cm. The incident light wave- 
length was 1.55 um. Under these conditions, the 
applied electric field was obtained by solving the 
Poisson equation. From the applied electric 
field, the envelope carrier wavefunction was 
calculated. The wavefunction was expressed 
with the Airy function. The absorption co- 
efficient was obtained from the wavefunction 
and the optical field». 

Figure 3 shows: the results of the calcula- 
tions. If the layer is too thin, a small light 
confinement in absorption layer reduces the 
extinction ratio. However, if the layer is too 
thick, the absorption is also reduced because the 
electric field applied to the absorption layer is 
weak. According to the above results, the oper- 
ating voltage is expected to be -3 V when the PL 
wavelength of the absorption layer is 1.43 um 
and the thickness is 0.13 ym. 
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4. Design for low chirp characteristic 

Practical MI-DFB lasers have a finite front 
facet reflectivity. A feedback light does greatly 
affect the lasing condition of a DFB laser. 
Therefore the authors modelled the MI-DFB 
laser to be equivalent to a DFB laser whose 
complex front facet reflectivity is changed by 
applying a voltage to be EA-MOD region. The 
complex front facet reflectivity is represented 
by an absorption coefficient, the reflectivity of 
the real front facet, and a light wave propagation 
constant in the EA-MOD region. 

To analyze the chirp behavior of the DFB 
laser, the authors calculated the dependence of 
the lasing wavelength on DC applied voltage, and 
then applied these results to the chirp behavior 
with a frequency response of the wavelength 
shift. To calculate a lasing wavelength change 
for an extinction ratio of -10 dB, the authors 
considered the change in complex front reflec- 
tivity, the threshold carrier density in the DFB- 
LD region, and kL, where xk is the coupling 
coefficient of the corrugation and L is the length 
of the DFB-LD region. The frequency response is 
calculated by rate equations which included a 
photon lifetime modulation scheme. The a 
parameter of the EA-MOD region and the DFB- 
LD region was 1” and 4 respectively. Heat 
effects were neglected. 

According to the calculated frequency 
response of a lasing wavelength shift, the chirp 
behavior under modulation was almost the same 
as the dependence of the lasing wavelength on 
DC voltage applied to the EA-MOD region, when 
the relaxation frequency of the DFB-LD region 
was sufficiently higher than the modulation fre- 
quency. The calculated chirp behavior agreed 
well with observed temporal behaviors using 
the streak camera system®, Therefore, the cal- 
culation is valid for designing low-chirp devices. 

As the lasing wavelength depends on the 
facet phases, the low-chirp operation was 
estimated using the yield of a wavelength shift 
for -10 dB DC attenuation. The chirp was found 
to decrease as the front facet reflectivity de- 
creased and as KL increased. However a large KL 
induces multi-mode operation due to a strong 
spatial hole burning®? , and there is difficulty in 


117 


T. Watanabe et al.: Low Drive Voltage and Low Chirp Modulator Integrated DFB- - - 


achieving a reflectivity of less than 0.1 percent 
using conventional techniques. Then according 
to the above analysis, a KL of 1.5 and a reflec- 
tivity of 0.1 percent was used because the yield 
of the maximum wavelength shift of less than 
0.01 nm was estimated to be over 50 percent. 


5. Fabrication 

The device was fabricated in four steps, 
three of which were liquid phase epitaxy (LPE) 
and one of which was metalorganic vapor phase 
epitaxy (MOVPE) as shown in Fig. 4". 

First, the DFB-LD region was fabricated. 
The first order corrugation for a wavelength of 
1.55 wm was fabricated partially on the sub- 
strate of the DFB-LD area. Next, a guide layer, a 
thin InP etch-arresting layer, an active layer, and 
a Clad layer were fabricated with LPE. The guide 
layer and the active layer were InGaAsP having 
a PL wavelength of 1.1 wm and 1.54 um respec- 
tively. The clad layer and the etch-arresting layer 
were InP. 

Next the EA-MOD region was fabricated. 


p-InP clad layer 


InP etch-arresting Active layer 


layer 


Guide layer 
Ist order _ n-InP substrate 
corrugatoin 
1) lst LPE growth 
SiO, p-InP clad layer 


Absorption layer 


2) 2nd LPE growth 


Contact layer 


p-InP clad layer 


3) 8rd LPE growth 


S.1.-InP 


4) MOVPE growth 


Fig. 4—Growth procedure. 
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The clad and the active layer on the EA-MOD 
area were etched off using a SiO, mask over the 
DFB-LD area. The guide layer was left selective- 
ly because of the etch-arresting layer. After 
etching, an InGaAsP absorption layer having a PL 
wavelength of 1.43 ym and an InP clad layer 
were refabricated using LPE. A p-clad layer and 
a contact layer were fabricated entirely using 
LPE after removing the SiO, mask over the 
DFB-LD area. 

Next a stripe structure was fabricated with 
selective etching on both the DFB-LD and the 
EA-MOD parts using the same SiO, mask. This 
stripe was buried in semi-insulate InP using 
MOVPE. Next the clad layer in the joint area 
was removed to establish sufficient electrical 
isolation. The resistance was 2 MQ by removing 
the contact layer at the region. 

Finally SiN film was fabricated on the front 
facet using a well-controlled deposition process. 
To achieve a reflectivity of 0.1 percent, the 
authors optimized the thickness of the SiN film 
on the front facet, taking into account again the 
two-dimensional optical field and the refractive 
index of the absorption layer whose PL wave- 
length was 1.43 ym. 


6. Characteristics 

Figure 5 shows the total light output power 
versus the injection current characteristic. The 
threshold current of the device is 24 mA. Out- 
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Fig. 5—Dependence of total light output power on the 
applied voltage. 
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Fig. 6—Extinction ratio of fiber coupled power. 
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Fig. 7—Small signal frequency response. 


put power is 11mW at 100mA. The slope 
efficiency is 0.14mW/mA. The total output 
power from the front facet includes the leak 
power and the guided power. The guided power 
is almost all absorbed under enough applied volt- 
age, about -5 V, and most of the other power is 
the leak power. The coupling efficiency is 
estimated to be 90 percent from the output 
power at no applied voltage and one at -5 V. The 
effective output power is 10 mW at 100 mA. 
These values are comparable to those of the 
previously reported device’. 

Only the effective guided power can coupled 
with fiber. The fiber output power is 1.5 mW. 
The extinction ratio of the fiber coupled output 
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Fig. 8—Light spectrum under 10 Gbit/s NRZ random 
modulation. 


power is shown in Fig. 6. The required voltage 
for the extinction ratio of -13 dB is only -3 ve. 
These results almost agree with the calculated 
extinction ratio of a solitary EA-MOD. This is 
a significant improvement over the previous 
result!?. 

Figure 7 shows the small signal frequency 
response under an applied voltage of -1 V and 
an injection current of 150mA. The cutoff 
frequency is 10.3 GHz, which is almost equal to 
the estimated cutoff frequency of the parasitic 
capacitance of 0.55 pF. 

Figure 8 shows the spectrum of the device 
under 10 Gbit/s NRZ pseudo random modula- 
tion. The measured chirp width at 15 dB down 
is only 0.016 nm (0.004 nm at -3 dB). This is a 
drastic reduction of the chirp of 0.4 to 0.8 nm 
in direct modulation of a DFB-LD. 

The temporary wavelength chirp is also 
important for optical transmission. Figure 9 
shows the temporary image of the wavelength 
chirp under 10 Gbit/s NRZ (1,0) fixed pattern. 
The current injected to the DFB-LD region is 
100 mA. The shift width of the chirp is very 
small. 

Figure 10 shows the eye pattern of the inte- 
grated device under 10 Gbit/s NRZ pseudo 
random modulation. The eye pattern of the 
integrated device has no peaking and no jitter. 
These results indicate that the integrated device 
can be used in high bit-rate long-haul systems. 
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Fig. 9-Temporary image of wavelength chirping at 
10 Gbit/s (1,0) fixed modulation. 


100 ps 


Fig. 10—Eye pattern at 10 Gbit/s NRZ random 
modulation. 


7. The transmission span 

The authors considered a tolerable dispersion 
in transmission for the optimized MI-DFB laser. 
Okiyama et. al. previously reported bit error rate 
(BER) characteristics using our previous MI- 
DFB laser for 10 Gbit/s transmission over 
64.8 km of dispersion shift fiber®’. This MI-DFB 
laser had a large chirp width of 0.01 nm at 3 dB 
down, but there was no power penalty for a 
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dispersion of 63 ps/nm. 

A tolerable dispersion for a power penalty 
of 0.5 dB was estimated to be 230 ps/nm. The 
newly developed MI-DFB laser has a very low 
chirp width of 0.016 nm at -15 dB (0.004 nm at 
-3 dB). The measured chirp width is drastically 
reduced compared to that of the previous MI- 
DFB laser’ Therefore a tolerable dispersion of 
much larger than 230 ps/nm can be expected for 
the latest ultra-low chirp MI-DFB. 


8. Conclusion 

This report discusses a low drive voltage and 
low chirp modulator integrated DFB laser. The 
authors optimized the composition and thick- 
ness of the absorption layer taking into account 
the absorption loss and absorption efficiency. 

As a result, the effective output power was 
10 mW at an injection current of 100 mA, and 
the extinction ratio was -13 dB at -3 V. A drive 
voltage of -3 V is usable for practical systems. 
The AR film on a front facet was optimized 
taking account of the two-dimensional optical 
field and the reflective index of the optimized 
absorption layer. The optimization of AR film 
gives a chirp of 0.016nm at -15 dB under 
10 Gbit/s NRZ pseudo random modulation. 
This low chirp is acceptable for a dispersion of 
larger than 230ps/nm with 10 Gbit/s NRZ 
pseudo random modulation. 

The new MI-DFB laser having low drive volt- 
age and low chirp performance is promising for 
use in the generation of multiple-gigabit, long- 
haul optical fiber transmission systems. 
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